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ABSTRACT 

This  study  is  a  step-by-step  process  for  examining  the 
need  and  probable  effectiveness  of  new  Coast  Guard  SAR 
facilities.  The  vehicle  for  the  study  is  a  proposed  Coast 
Guard  Air  Station  to  be  located  at  Areata,  California,  but 
the  process  is  adaptable  to  any  locale  or  any  type  of  SAR 
facility.  Data  from  historic  SAR  incidents  was  analysed  and 
the  SAR  environment  of  the  area  to  be  served  by  the  air 
station  was  simulated.  The  study  indicated  that  the  air 
station  will  reduce  response  time  to  maritime  SAR  incidents 
15  to  25  minutes  without  considering  possible  delays  caused 
by  low  visibility  at  the  Areata  airport.  Additional  study 
in  specific  areas  was  recommended. 
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!•  INTRODUCTION 

Probably  the  most  well-known  role  of  the  United  States 
Coast  Guard  is  that  of  Search  and  Rescue  (SAR) ,  Originally 
the  Revenue  Marine  Service  considered  SAR  as  a  collateral 
duty,  but  by  1837  its  ships  were  actively  cruising  the 
coasts  of  the  United  States  in  order  to  locate  and  assist 
mariners  in  distress.  Another  organization,  the  U.S. 
Lifesaving  Service,  also  had  a  long  and  notable  record  of 
saving  crewmen  whose  ships  had  gone  aground.  In  1915  these 
two  services  were  merged  to  form  the  U.S.  Coast  Guard,  and 
SAR  was  a  tradition  of  long  standing  that  was  incorporated 
into  the  duties  of  this  new  service. 

This  tradition  is  codified  into  the  laws  of  the  United 
States  in  the  United  States  Code  as  follows: 


1»_  USC  2  "The  Coast  Guard  .  .  .  shall  develop, 
establish,  maintain,  and  operate  .  .  .  rescue  facilities  for 
the   promotion   of   safety  on,  under,  and  over  the  high  seas 

and  waters  subject  to  the  jurisdiction  ox  the  United   States 

ii 

.... 

1_  USC  88  "In  order  to  render  aid  to  distressed 
persons,  vessels,  and  aircraft  on  and  under  the  high  seas 
and  on  and  under  the  waters  over  which  the  United  States  has 
jurisdiction  and  in  order  to  render  aid  to  persons  and 
property  imperiled  by  flood,  the  Coast  Guard  may  perform  any 
and  all  acts  necessary  to  rescue  and  aid  persons  and  protect 
and  save  property."  "The  Coast  Guard  mav  render  aid  to 
persons  and  protect  and  save  property  at  any  time  and  at  any 
place  at  which  Coast  Guard  lacilities  and  personnel  are 
available  and  can  be  effectively  utilized." 

14  USC  14J[  "  The  Coast  Guard  may,  when  so  requested 
by  proper  authority,  utilize  its  personnel  and  facilities  to 
assist  any  Federal  agency,  State,  Territory,  possession.,   or 

Political   subdivision  thereof,  or  the  District  of  Columbia, 
o   perform   any   activity   for   which   such   personnel   and 
facilities  are  especially  qualified.18 


In  order  to  perform  its  Search  and  Rescue  mission,  the 
Coast  Guard  must  allocate  personnel  and  equipment  resources 
to  various  locations  around  the  United  States  to  provide  its 
services   when   needed   by   persons  or  property  in  distress. 


Because  the  Coast  Guard  SAR  mission  is  primarily  concerned 
with  the  saving  of  life  and  property  after  a  distress 
situation  develops,  SAR  is  strongly  time  oriented.  The 
longer  a  person  or  his  property  is  exposed  to  a  hostile 
environment,  the  greater  are  the  chances  that  an  incident 
will  turn  into  a  disaster,  resulting  in  the  loss  of  life 
and/or  property.1  For  this  reason,  the  time  for  a  Coast 
Guard  unit  to  respond  to  a  distress  presents  a  simple 
measure  of  effectiveness  for  its  SAR  activity. 

It  is  conceivable,  of  course,  that  sufficient  resources 
could  be  made  available  in  the  Coast  Guard  budget  to  ensure 
that  a  Coast  Guard  unit  would  arrive  at  the  scene  of  any 
distress  within  any  given  time  increment  anywhere  in  the 
United  States.  As  the  specified  time  increment  decreases 
and  the  distance  from  the  "home"  location  of  the  Coast  Guard 
resources  to  the  position  of  distress  increases,  the  shorter 
response  time  entails  additional  costs.  Eventually,  the 
marginal  improvement  in  response  becomes  prohibitively 
expensive. 

There  is  some  difficulty  in  defining  an  acceptable 
response  time  for  SAR.  One  effort  to  define  such  a  time 
increment  was  given  for  Coast  Guard  aviation  units  as 
follows: 


"SAR  aviation  units  shall,  within  the  Harbor  and 
Coastal  Zone,  be  capable  of  flying  to  the  scene  of  75 
percent  of  the  assistance  cases  within  one-half  hour  and  90 
percent  within  one  hour  and  be  capable  of  recovering 
(rescuing)  persons  in  distress."  2 


The  Harbor  and  Coastal  Zone  was  defined  in  the  Coast  Guard 
Aviation  Plan  as  being  the  zone  including  rivers,  bays,  and 
inlets  out  to  150  miles  seaward. 

The   allocation  of  scarce  SAR  resources  is  an  important 
facet  of  the  Coast  Guard's   SAR   role.    Resources   must   be 


stationed  to  ensure  that  they  can  be  economically  operated 
and  are  suited  to  the  natural  and  SAR  environment  of  their 
assigned  area  while  still  achieving  the  desired  level  of 
response.  At  the  same  time,  care  must  be  exercised  to  avoid 
positioning  too  many  resources  in  any  area.  This  results  in 
"overkill",  inefficiencies,  and  additional  expenses.  In 
order  to  properly  determine  the  optimum  type  and  mix  of 
resources,  their  location,  and  their  manning  levels,  an 
analysis  should  be  made  of  all  factors  bearing  on  their 
activities  in  an  operating  area.  The  analysis  should 
investigate  probable  numbers  and  types  of  distress  incidents 
requiring  assistance  by  Coast  Guard  resources,  the  natural 
environment,  complimentary  search  and  rescue  resources  in 
the  area,  and  any  other  pertinent  factors. 

The  purpose  of  this  thesis  was  to  supply  an  example  of 
such  an  analysis  to  the  needs  for  Coast  Guard 
search-and-rescue  resources  in  an  area  of  Northern 
California.  The  analysis  was  directed  toward  determining  if 
an  additional  SAR  resource  was  reguired  in  the  area  --  an 
air  station  at  the  Arcata-Eureka  Airport.  It  was  written, 
however,  in  such  a  manner  that  it  could  be  applied  to  any 
change  in  Coast  Guard  resources  for  any  area. 

The  current  plan,  as  promulgated  in  Part  Two  of  the 
Coast  Guard  Aviation  Plan  (CG-38G-2) ,  is  to  establish  an  air 
station  at  Areata  with  three  HH-52A  helicopters  in  FY  1978, 
budgetary  liaitations  permitting.  The  rationale  for  the  new 
air  station  is  presented  below: 

"Since  FY  1963  there  has  been  a  substantial  growth  in 
the  number  of  SAR  cases  within  a  radius  of  75  miles  of  the 
Humbolt  Bay  area  of  Northern  California.  Case  totals  have 
increased  24*  in  ten  years  [sic],  from  182  in  FY  1963  to  449 
in  FY  1972.  The  distance  from  the  nearest  om  182  in  FY  1963 
to  449  in  FY  1972.  The  distance  from  the  nearest  Coast 
Guard  air  units  precludes  our  ability  to  provide  timely 
response  with  adequate  recovery  capability.  San  Francisco 
Air  Station  is  210  miles  to  the  south  and  Astoria  is  325 
miles  north  of  Humbolt  Bdy.  Even  with  the  establishment  or 
North  Bend  Air  Station  in  FY  1974  the  nearest  aircrart 
response   (HH-52A's)  will  still  be   1.65  miles  away,  almost  a 


two-hour  flight  from  Coos  Bay. 

The  National  Planning  Association*s  "Marine  Activities 
Forecast"  for  Northern  California  appears  conservative  for 
the  Humbolt  Bay  region;  their  FY  1930  workload  projections 
may  well  be  upon  us  in  tne  early  1970's. 

Public,  private,  and  political  pressure  has  already 
been  exerted  to  bring  Coast  Guard  rescue  aircraft  to  the 
Humbolt  Bay  area*  in  fact  in  1966  Humbolt  County  authorized 
12.8  acres  of  land  at  the  Areata  County  Airport  for  Coast 
Guard  use  on  the  condition  that  a  helicopter  base  be  located 
at  said  airport.  A  reasonable  case  for  establishing  such  a 
unit  can  be  generated  on  the  basis  of  the  increase  in  SAR 
case  load,  aralysis  of  the  nature  of  Coast  Guard  surface 
unit  responses  in  that  area,  and  projected  growth."3 


Two  mission  areas  of  the  Coast  Guard  were  forecasted 
for  the  proposed  station.  They  were  SAR  and  Enforcement  of 
Laws  and  Treaties.*  This  station  is  rather  unique  in  that 
it  is  to  have  such  a  limited  number  of  mission  areas  as  its 
responsibility,  a  fact  which  greatly  simplifies  analysis. 

Because  of  its  emphasis  on  the  Search  and  Rescue  role, 
the  Coast  Guard  has  developed  a  standardized  reporting 
system  to  support  this  program.  Each  time  a  Coast  Guard 
facility  is  used  to  aid  persons  or  property  the  responsible 
Coast  Guard  official  is  required  to  submit  an  Assistance 
Report  (CG  3272) .  The  format  of  this  report  has  changed 
several  times  over  the  years  but,  essentially,  the  same 
information  is  reported  in  each.  A  complete  description  of 
the  property  in  distress  is  given,  including  the  cause  of 
distress.  "Nature  of  Distress"  is  a  term  which  is  used  in 
the  Assistance  Report  to  indicate  the  problem  encountered  by 
the  persons  or  property  that  requires  Coast  Guard 
assistance.  One  example  of  nature  of  distress  is  "vessel 
disabled  and  adrift."  Another  is  "personnel  injury."  When 
the  Coast  Guard  receives  a  report  that  its  assistance  is 
required,  the  nature  of  distress  is  one  important 
consideration  used  to  determine  which  resource  to  dispatch 
to  aid  in  the  situation.  Further  information  includes  a 
brief  description  of  the  weather  encountered  along  with  the 
absolute  position  and  the  position  of  distress   relative   to 
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the  responding  Coast  Guard  facility.  For  the  purpose  of  the 
thesis,  these  items  which  are  essentially  beyond  the  control 
of  man  are  termed  "states  of  nature."  In  addition, 
significant  times  are  recorded,  along  with  the  actions  of 
the  Coast  Guard  facility.  Finally,  the  relative  success  of 
the  actions  cf  the  Coast  Guard  is  reported  in  terms  of  lives 
saved,  helped,  or  lost  and  the  value  of  property  assisted. 
These  data  are  compiled  and  stored  in  Coast  Guard 
Headguarters  for  subseguent  use  in  analysis  and  decision 
making.  A  sample  Assistance  Report  is  included  as  Appendix 
A. 

The  Assistance  Reports  for  Northern  California,  as 
transcribed  to  the  permanent  Coast  Guard  records,  formed  the 
data  base  for  this  thesis.  Originally  the  data  was  in  the 
form  of  cards  obtained  from  the  Twelfth  Coast  Guard 
District.  This  data  was  complete  in  that  it  contained  all 
entries  from  the  Assistance  Reports  for  the  years  1971 
through  1973.  Later,  additional  data  was  received  in  the 
form  of  cards  for  the  years  1970  and  1974  from  Coast  Guard 
Headguarters.  This  data  had  been  transcribed  from  file  tape 
to  cards.  The  data  was  not  complete  in  that  all  the  data 
submitted  on  the  Assistance  Report  was  not  included.  The 
information  concerning  month  of  distress,  position  of 
distress,  and  time  of  day  was  missing  from  this  data.  When 
data  was  missing  from  some  years,  that  fact  was  indicated  in 
the  analysis  chapter  of  the  thesis. 
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II*   AREA  STUDY 

The  term  used  to  describe  the  geographical  region  to  be 
serviced  by  the  proposed  Coast  Guard  facility  was  "area  of 
interest."  This  area  of  interest  should  be  large  enough  to 
include  the  normal  operations  of  the  Coast  Guard  facility, 
taking  into  account  other  Coast  Guard  units  in  the  area, 
geographical  considerations,  and  limitations  imposed  by  the 
type  of  Coast  Guard  resource  being  considered.  For  the 
proposed  Coast  Guard  Air  Station  at  the  Areata- Eureka 
airport,  tne  area  of  interest  was  defined  as  that  region 
between  40-00  and  42-00  North  Latitude  and  between  the 
mountainous  region  approximately  thirty  miles  inland  from 
the  coastline  to  300  miles  offshore.  Latitude  42-00N  is  the 
position  of  the  California-Oregon  border  and  forms  the 
boundary  between  the  Twelfth  Coast  Guard  District  and  tne 
Thirteenth  Coast  Guard  District.  The  region  measures 
approximately  sixty  miles  north  and  south  and  more  than  300 
miles  in  its  east-west  dimensions. 

The  most  significant  maritime  feature  in  the  area  of 
interest  is  Humboit  Bay,  which  is  a  large  sheltered  bay  and 
the  third  largest  natural  harbor  in  California.  It  is  the 
most  important  harbor  on  the  west  coast  between  San 
Francisco  and  the  Columbia  River.  The  cities  of  Eureka  and 
Areata  are  located  near  this  bay.  Eureka  is  the  third  most 
valuable  fishing  port  in  California.  The  major  industry  for 
Areata  is  lumbering  and  lumber  products.  A  third  important 
town  in  the  area  of  interest  is  Crescent  City,  located 
approximately  forty  miles  north  of  Humbolt  Bay.  The  major 
industries  in  Crescent  City  are  lumbering  and  dairying. 

Two  large  river  systems  empty  from  inland  regions  into 
the  ocean  in  the  area  of  interest.  The  more  important  one, 
the  Klamath  River,  empties  approximately  30  miles   north   of 
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Humbolt  Bay.  This  river  is  important  because  nearly 
one-third  of  the  incidents  involving  boaters  within  the  area 
of  interest  occurred  within  five  miles  of  its  mouth.  The 
Eel  River  ends  approximately  five  miles  south  of  Humbolt 
Bay.  Although  few  boats  are  assisted  by  the  Coast  Guard 
near  this  river,  it  has  flooded  during  various  years  and 
Coast  Guard  forces  have  participated  and  assisted  local 
authorities  in  flood  relief  operations. 

A  description  of  the  geography  and  climate  in  the  area 
of  interest  is  provided  by  the  Department  of  Commerce  Local 
Climatologicai  Data  description: 


"There  are  no  hills  in  Eureka  of  any  conseguence.  The 
land  slopes  upward  gently  from  the  Bay  towards  the  Coast 
Range,  which  begins  about  3  miles  east  of  the  station 
(Eureka)  and  reaches  the  top  of  its  first  ridge 
approximately  10  miles  to  the  east.  The  average  elevation 
of  the  ridge  is  about  2,000  feet.  This  ridge  extends  in  a 
semi-circle  from  a  point  20  miles  north  of  Eureka  to  a  point 
25  miles  south. 

The  climate  of  Eureka  being  completely  maritime,  high 
humidity  prevails  the  entire  year,  which  is  divided  into  the 
"rainy"  season  and  the  "dry"  season.  The  rainy  season  begins 
in  October  and  continues  through  April.  About  90  percent  of 
the  year's  precipitation  falls  during  this  period.  The  dry 
season  extends  from  May  through  September  and  is  marked  by 
considerable  fog  or  low  cloudiness.  Usually,  however,  the 
fog  clears  in  the  late  forenoon  with  the  early  afternoons 
generally  sunny. 

Temperatures  are  moderate  the  entire  year.  Although 
the  highest  ever  recorded  was  85  degrees,  and  the  lowest  20 
degrees.  the  usual  range  is  from  a  low  of  about  35  degrees 
to  a  nigh  of  about  75  degrees.  The  daily  range  of 
temperature  averages  from  about  9  degrees  in  the  summer 
months  to  13  degrees  in  the  winter  months,  and  is 
occasionally  not  over  2  to  3  degrees."5 


The  temperature  of  the  waters  offshore  in  the  Eureka 
area  range  between  50  degrees  and  55  degrees  Fahrenheit. 

The  area  of  interest  includes  portions  of  two  counties 
of  California,  Humbolt  County  and  Del  Norte  County.  The 
total  number  of  boats  registered  in  these  two  counties  for 
the  years. 1970  through  1974  were: 
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1L2I  1221  1222  1221  122J1 

2037       2919       4203       4773       5232 

The  low  figures  for  1970  and  1971  were  due  to  the  fact  that 
a  change  in  the  registration  procedure  took  place  between 
those  years  and  subseguent  years.  Typically,  77  to  78 
percent  of  the  registered  boats  are  registered  in  Huabolt 
County.  No  information  was  available  for  the  number  of 
documented  vessels  in  the  area. 

Data  was  not  made  available  which  gave  detailed 
information  concerning  the  lengths  of  the  vessels  in  Humbolt 
and  Del  Norte  counties.  The  usage  of  those  vessels  which 
were  registered  also  was  not  given.  The  state-wide 
distribution  of  lengths  of  registered  boats  was  reported  as 
follows: 6 


JjT/C 

1221 

1221 

1221 

1974 

Less  than 
16  feet 

64.0% 

61.5% 

60.5% 

58.2% 

57.5% 

16  feet  to 
26  feet 

32.0% 

32.5% 

34.0% 

36.2% 

36.9% 

26  feet  to 
40  feet 

4.5% 

4.5% 

5.0% 

5.0% 

4.9% 

40  feet  to 
65  feet 

0.5% 

0.5% 

0.5% 

0.6% 

0.7% 

The  number  of  Search  and  Rescue  incidents  answered  by 
the  Coast  Guard  in  the  area  of  interest  during  the  Fiscal 
Years  1970  through  1974  was  1650,  including  multi-unit  cases 
in  which  more  than  one  Coast  Guard  resource  answered  a  call 
for  distress.  The  data  from  some  of  these  cases  was 
unusable,  resulting  in  the  discarding  of  some  of  the 
incident  reports.  Such  unusable  data  included  reports  for 
which  the  severity  to  personnel  or  property  was  coded  as 
unknown,  data  in  which  the  position  of  distress  was  coded  as 
unknown,  or  data  which   was   incomplete.    Culling   of   such 
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erroneous   data   left   a  total  of  1468  cases  for  five  years. 

The  number  of  cases   occurring   in  any   year  is   given   as 
follows: 

HIS            1221     1111  1111            1111 

327        264        360  195        322 

The  distribution  of  lengths  of  vessels  assisted  by   the 
Coast  Guard  during  the  five  years  were: 

Not  a  Vessel  3.4% 

Less  than  16  feet  32.2% 

16  feet  less  than  26  feet  15.4% 

26  feet  less  than  40  feet  21.5% 

40  feet  less  than  65  feet  20.5% 

65  feet  less  than  100  feet  2.6% 

100  feet  less  than  200  feet  1.1% 

200  feet  less  than  300  feet  0.2% 

Greater  than  300  feet  2.0% 

Unknown  1.1% 
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The  distribution  of  the  usages  of  the  vessels   assisted 
during  those  years  was: 

Recreation  47.'!% 

Fishing  43.5% 

Towing  -   1.0% 

Cargo  1.6% 

Other  1.495 

No  other  code  for  usage  appeared  wore  than  one  percent  of 
the  time.  The  figures  do  not  total  100%  because  in  the 
Assistance  Report  a  usage  code  can  be  assigned  to  cases 
other  than  vessels.  During  the  five  years  there  were 
approximately  13  cases  involving  aircraft,  3  cases  involving 
land  vehicles,  5  cases  involving  land  structures  and  47 
cases  involving  personnel  only  in  the  area  of  interest. 
Charts  representing  the  geographic  distribution  of 
historical  SAR  cases  are  included  in  Appendix  F. 

Most  of  the  Coast  Guard  facilities  now  rendering  search 
and  rescue  services  in  the  area  of  interest  are  under  the 
command  and  control  of  Group  Humbolt  Bay,  which  reports  to 
the  Twelfth  Coast  Guard  District.  The  medium-endurance 
cutter  COMANCHE  (WMEC  202) ,  home-ported  in  Eureka,  is  also 
under  the  direct  control  of  the  District  Headquarters. 
Included  in  the  facilities  of  Group  Humbolt  Bay  are  the 
95-foot  WPB  CAPE  CARTER,,  ported  in  Crescent  City  and  two 
44-foot  MLBs  stationed  in  Eureka.  During  the  salmon  season, 
from  approximately  1  July  to  15  October  of  each  year,  the 
boating  traffic  on  the  Klamath  River  rises  dramatically. 
Swift   river  currents  and  the  interaction  of  the  ocean  swell 
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causes  hazardous  conditions  on  the  river.  During  this 
period,  Station  Klamath  River  is  established  with  a 
complement  of  approximately  six  men  and  trailerable  25-foot 
and  17-foot  boats. 

Coast  Guard  air  support  for  the  area  of  interest  is 
currently  provided  by.  the  Coast  Guard  Air  Station  at  San 
Francisco  International  Airport.  That  location  is 
approximately  210  miles  away  from  Eureka  via  Federal 
Airways.  A  HH-52A  helicopter  requires  approximately  three 
hours  to  reach  the  Arcata-Eureka  airport  from  San  Francisco 
and  requires  refuelling  prior  to  continuing  on  a  SAR  case. 
A  fixed-wing  HC-130  or  HU-16  aircraft  from  San  Francisco  can 
be  over  Eureka  approximately  one  and  a  half  hours  after 
notification. 

Single-engined  HH-52A  helicopters  are  limited  in  their 
operations  and  must  remain  within  25  miles  of  the  shoreline 
unless  they  are  escorted  by  fixed-wing  aircraft,  in  which 
case  they  are  limited  to  100  miles  offshore.  On  several 
occasions,  especially  in  cases  requiring  medical  evacuation 
of  seriously-injured  or  critically-ill  persons,  HH-3F 
helicopters  from  the  Coast  Guard  Air  Station  at  Astoria, 
Oregon  or  San  Diego,  California  were  flown  to  Eureka  by  one 
crew  and  the  medevac  mission  was  performed  by  a  second  crew. 
The  time  required  to  fly  a  helicopter  to  Eureka  is 
approximately  four  hours.  HH-3F  helicopters  can  perform 
missions  up  to  300  miles  offshore.  Assistance  to  personnel 
or  to  property  greater  than  300  miles  offshore  must  be 
provided  by  surface  vessels  or  by  air  dropping  equipment  cr 
people  from  fixed-wing  aircraft.  These  aircraft  can,  of 
course,  assist  vessels  which  are  not  more  than  300  miles 
offshore  in  a  similiar  manner. 

In  addition  to  the  Coast  Guard  units  in  the  Eureka 
area,  there  are  other  resources  which  are   used   for   search 
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and  rescue,  particularly  for  incidents  on  the  land.  The 
Humbolt  County  Sheriff's  Department  reports  that  it  responds 
to  approximately  thirty  inland  SAK  incidents  per  year,  and 
that  this  number  does  not  appear  to  be  growing.  The 
Sheriff's  Department  calls  upon  the  Civil  Air  Patrol,  the 
Ski  Patrol,  the  Marine  Posse,  the  Jeep  Patrol,  and  the 
Mounted  Posse  to  assist  in  inland  SAR  incidents.  It 
requires  approximately  three  hours  between  notification  and 
the  time  a  Civil  Air  Patrol  unit  is  available,  weather 
permitting.7  The  Forest  Service  reportedly  has  use  of  a 
small  helicopter  which  they  operate  in  the  area  during  the 
summer  months.  All  these  resources  have  worked  together 
with  Coast  Guard  units  in  the  past,  both  in  operations  along 
the  coastline  and  overland. 

An  important  aspect  of  the  study  of  the  effectiveness 
of  a  Coast  Guard  air  station  is  the  weather.  Having  an  air 
station  in  an  area  serves  no  purpose  if  the  aircraft 
stationed  there  are  grounded  in  fog.  The  Coast  Guard  Air 
Operations  Manual  (CG-333)  specifies  the  minimum  ceiling  and 
visibility  requirements  for  Coast  Guard  aircraft  departing 
an  airfield.  In  the  case  of  a  single-engined  helicopter  at 
Areata  airport,  the  appropriate  criteria  are  those  for 
airfields  with  a  departure  alternate,  as  defined  in  the  Air 
Operations  Manual,  since  an  airfield  at  Crescent  City  has  a 
precision  instrument  approach  where  the  helicopter  could 
land  if  necessary.  This  airfield  is  approximately  forty 
minutes  flying  time  away  from  Areata.  If  the  weather  for 
this  airfield  is  forecasted  to  be  600  foot  ceiling  with  two 
miles  visibility  for  one  hour  and  forty  minutes  after 
takeoff,  the  Coast  Guard  helicopter  could  depart  Areata  if 
the  visibility  there  is  at  least  one-fourth  mile  or  1600 
feet  visibility  (RVR  16) .  An  exception  is  granted  "When  the 
immediate  urgency  of  the  mission  dictates  .  .  .  Due 
consideration  must  be  given  to  equipment  capabilities, 
alternate    field    availability,    and   the  circumstances 
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constituting  the  argent  requirement. "e 

In  1968  the  National  PJureau  of  Standards  published  a 
report  of  the  low  visibility  conditions  at  the  Areata 
airport  for  the  period  1957  through  1967.9  Here  follows  one 
interpretation  of  that  report  as  it  applies  to  the  proposed 
air  station  at  Areata. 

The  visibility  at  Areata  airport  was  reported  as  being 
below  2100  feet  by  day  or  4300  feet  by  night  (defined  in  the 
report  as  "low  visibility")  as  the  number  of  hours  per 
month.  The  yearly  average  varied  between  a  low  of  26.3  to  a 
high  of  63.3,  with  the  number  of  nighttime  hours 
approximately  double  the  number  of  daytime  hours  in  most 
cases.  The  greatest  number  of  low-visibility  hours  occurred 
in  the  July-January  semi-annual  period.  The  greatest  number 
of  low-visibility  hours  for  any  month  were  September  or 
October,  a  period  in  which  approximately  eighteen  percent  of 
the  historic  SAR  incidents  occurred. 

The  report  contained  the  summary  of  a  three-year  Runway 
Visual  Range  (RVR)  study.  The  times  during  which  RVR  was 
less  than  1000  feet  had  a  mean  of  0.28  hours  and  a  median  of 
0.08  for  day  periods  and  a  median  of  0.15  for  night  periods. 
The  times  ranged  as  high  as  8.3  hours.  In  a  table  of  times 
during  which  the  RVR  was  below  1000  feet,  the  periods  of  low 
visibility  seemed  to  be  centered  between  approximately  one 
hour  before  sunrise  and  noon  and  between  one  hour  before 
sunset  and  midnight.  Approximately  forty  percent  of  the 
historic  SAR  cases  were  during  these  hours  and  they  had 
almost  a  uniform  distribution  of  occurrences  between  the 
hours  of  0800  and  1800  local  time. 

While  the  exception  to  the  takeoff  minimums  in  the  case 
of  an  urgent  mission  permits  a  Coast  Guard  helicopter  to 
take   off  in  low- visibility  conditions,  the  effectiveness  of 
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that  aircraft  is  severely  reduced  in  very  low  visibility. 
Approximately  ten  percent  of  the  historic  SAR  cases  occurred 
in  reported  visibilities  of  less  than  one-half  aile.  The 
effect  of  the  weather  at  Areata  --  both  low  ceilings  and  low 
visibilities  —  upon  operations  of  the  proposed  air  station 
deserves  further  study. 
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Ill*   ANALYSIS  OF  HISTORIC  STATISTICS 

The  Coast  Guard  now  has  facilities  located  so  as  to  be 
able  to  respond  to  incidents  in  virtually  all  portions  of 
the  waters  near  the  seacoast  of  the  United  States.  As  Coast 
Guard  units  respond  to  calls  for  assistance  in  any  given 
area,  the  Assistance  Reports  these  facilities  submit 
represent  a  sample  of  the  conditions  which  will  be 
encountered  by  a  new  Coast  Guard  facility  located  in  that 
area.  One  of  the  first  steps  in  analysing  the  need  for  a  new 
Coast  Guard  facility  in  a  given  area  is  to  examine  the  data 
available  in  the  Coast  Guard  records. 

The  Coast  Guard  SAR  data  readily  lends  itself  to  being 
represented  in  freguency  distributions.  The  frequencies  of 
the  coded  data  are  made  discrete  rather  than  continuous  by 
the  coding  process.  Rather  than  reporting  the  exact 
distance  offshore  that  an  incident  occurs,  for  example,  the 
reporting  facility  fits  the  distance  offshore  into  one  of 
the  ten  codes  (numbers  zero  to  nine)  and  reports  this  coded 
number.  The  process  of  translating  the  distances  offshore, 
which  are  continuous,  into  coded  values  causes  the 
distribution  of  the  data  to  become  discrete.  In  order  to 
facilitate  analysis,  the  discrete  distributions  of  the  data 
can  be  graphically  represented  in  histograms.  Freguency 
distributions  for  weather  conditions,  nature  of  distress, 
month  and  time  of  distress,  and  distance  offshore  of 
distress  for  the  data  from  the  Areata  area  of  interest  are 
given  in  Appendix  B. 

The  characteristics  of  these  frequency  distributions 
are  descriptive  statistics.  A  very  useful  statistic  is  the 
mean,  or  the  average  value,  of  the  observations.  There  are 
several  useful  tests  concerning  the  means  of  the  data  which 
can  be  performed  using  both  the  data  from  the  area  of 
interest  and  the  data  frotr.  the  Coast  Guard  as  a  whole. 
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The  means  of  the  code  values  for  the  data  for  the  five 
years  are  tabulated  below.  The  data  which  are  tabulated  are 
those  elements  of  the  SAB  environment  which  occur  randomly. 
This  data  included  multi-unit  case  data  in  which  a  small 
percentage  of  observations  were  duplicated  because  more  than 
one  unit  was  reporting  the  same  data. 


TABLE  OF  YEARLY  MEANS 
HISTORICAL  DATA 


JS70 

1211 

1211 

1973 

1221 

Nature  of 
Distress 

7.8226 

10.6932 

5.3889 

10.5487 

6.7019 

Distance 
Offshore 

1.3639 

1.5795 

1.2972 

1.7333 

1.3043 

Personnel 
Severity 

1.5260 

1.5455 

1.3083 

1.3538 

1.21U3 

Property 
Severity 

1.5474 

1.4356 

1.3056 

1.3385 

1. 1553 

Sea 
State 

2. 1284 

2.4583 

2.0361 

1.7949 

1.6553 

Wind 
Velocity 

1.4648 

1.8  295 

1.0611 

1.4359 

1.254  7 

Visibility 

4.4924 

4.7992 

4,4444 

4. U 974 

4.6584 

Length 

2.6911 

2.9697 

2. 1944 

2.9744 

2.2547 

Time  of  Day 

12.8838 

12.5985 

12.3472 

13.3436 

13.2516 

Month 

6.3385 

7.5472 

6.8485 

The  means  of  the  data  were  compared  to  determine  if  the 
values  of  the  means  are  statistically  the  same  from  year  to 
year.  By  analysing  the  data,  it  may  be  possible  to  detect 
trends  such  as  cases  occurring  further  offshore  each  year. 
Patterns  such  as  severity  to  property  being  significantly 
greater  than  severity  to  personnel  may  also  be  revealed. 

In  order  to  examine  the  data  for  five  years,  a  one-way 
analysis  of  variance  of  corresponding  means  was  performed. 
The  hypothesis  being  tested  is  that  the  means  in  each  row  of 
data  are  jointly  equal.   The  implicit   assumption   was   made 
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that  the  distribution  of  the  means  of  the  data  was  normal  in 
all  cases.  The  alternate  hypotheses,  in  the  notation 
presented  in  Bolch  and  Huang  [Ref.  1],  are  shown  as: 


H    : 

u.       =    u.       =    u.      =    u.      = 

u. 

0 

12              3             4 

5 

H  :  Not  all  column  means  are  equal 
1  * 


The  appropriate  test  was  Fisher's  F-test  in   which   the 
F-statistic  is  computed  as  a  ratio  as  follows: 

F  =   Metwe§Hz.c2lHEiIi  Variance 
Within- (Column  Variance 

This   test   was   performed   using   the   Biomedical   Computer 

Program   BIMED01V  [Ref.  2]  on  the  IBM  360/67   at   the   Naval 

Postgraduate   School.    The   resulting   F-ratio  was  compared 

with   a   standard   table   of   percentage   points    for   the 

F-statistic  with  four  and  infinite  degrees  of  freedom  at  the 

0.95  confidence  level  to  determine  if  all  the   columns   were 

egual.    If   the   computed   ratio   of  F  was  greater  than  the 

tabulated  critical  value  of  F,  the  null  hypothesis   was  not 

accepted   and   it  was  assumed  that  the  means  of  the  data  for 

the  five  years   were   different   from   year   to   year.    The 

results  of  the  test  are  tabulated  below: 
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ANALYSIS  OF  MEANS 
(DATA  FOR  FIVE  YEARS) 

F-ratio      F-statistic 
Computed        Critical 

Nature  of  Distress  3.66  2.37 

Distance  Offshore  7.79  2.37 

Severity  to  personnel  12.18  2.37 

Severity  to  Property  12.07  2.37 

Sea  State  8.06  2.37 

Wind  Velocity  10. 44  2.37 

Visibility  1.64  2.37 

Length  of  Vessel  14.52  2.37 

Time  of  Notification  2.24  2.37 

Month  (3  years)  19.19  3.00 

In  every  case  except   visibility,   the   test   indicated 

that   the   means   for   the   data   for   the   five   years  were 

different.   In  those  instances  -where  the  calculated   F-ratio 

was  greater  than  the  critical  F-statistic,  it  is  possible  to 

discover  which  year's  data  was  significantly  different   from 

the  rest  by  estimation  of  the  linear  contrasts  of  the  column 

means.   This  test  is  a  pair-wise  comparison   of   all   column 

means   in   which   the   null   hypothesis  H  :  u.   -  u.   =  0  is 

0     i      j 

tested  for  all  possible  (i  j)  pair  combinations  of  the  data. 

A   graphic   depiction   of  the  results  of  these  tests  for  the 
historic  data  is  included  in  Appendix  D. 
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The  result  of  this  analysis  of  variance  can  be 
summarized  by  stating  that  while  conditions  in  an  area  of 
interest  may  be  essentially  the  same  from  year  to  year, 
there  are  usually  some  conditions  that  are  significantly 
different  in  any  one  year.  For  example,  the  value  of 
distance  offshore  that  SAP.  incidents  occur  may  be 
significantly  greater  for  one  year.  There  may  be  a  logical 
explanation  for  this  phenomena,  such  as  the  fish  being 
farther  offshore  that  year,  but  the  analyst  probably  does 
not  have  access  to  this  explanatory  information.  Ke  can 
only  discover  this  fact  statistically  and  plan  for  the 
recurrence  when  allocating  Coast  Guard  resources. 

The  next  step  in  analysing  the  data  for  the  area  of 
interest  is  to  compare  it  with  the  data  of  all  Coast  Guard 
cases  nationwide.  By  comparing  the  means  of  the  values  of 
sea  state,  for  example,  for  the  area  of  interest  with  the 
means  of  all  Coast  Guard  cases  in  any  year  it  is  possible  to 
infer  whether  the  seas  are  rougher  than  generally 
experienced  nationwide.  Similiarly,  by  comparing  the  means 
of  severity  to  persons,  it  can  be  inferred  that  more  rapid 
response  is  necessary  if  the  mean  for  the  area  of  interest- 
is  greater  than  for  the  United  States  as  a  whole. 

Due  to  the  fact  that  the  weather  conditions,  the 
descriptions  of  vessels  in  distress,  and  the  positions  of 
distress  appear  to  vary  in  one  single  area  from  one  year  to 
the  next,  an  attempt  was  made  to  approximate  the  "average" 
conditions  in  the  area  to  use  to  compare  with  the 
nation-wide  data.  The  effects  of  the  variations  could  be 
reduced  somewhat  by  pooling  the  data  and  using  the  mean  for 
the  five  years  in  the  area  of  interest.  This  mean  was 
compared  with  the  nation-wide  mean  for  each  year's  data. 
The  resulting  means  are  shown  in  the  following  table.  These 
means  also  include  multi-unit  case  data  in  which  two  or  more 
units  reported  the  same  data,  but  the  effect  on  the  validity 
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of   the   data   was  judged  to  he  insignificant  because  of  the 
relatively  small  number  of  such  cases. 


TABLE  OF  MEANS  OF  POOLED  DATA 
WITH  NATION-HIDE  DATA 


5YR 
Areata 

1971 

1972 

N§ti^)nwide_ 

1973 

Nature  of 
Distress 

7.8583 

13.9135 

13.1479 

12.9107 

Distance 
Offshore 

1.4221 

1.3629 

1.3649 

1  .3306 

Personnel 
Severity 

1.3849 

1.1685 

1.1677 

1.1579 

Property 
Severity 

1.3542 

1.  1394 

1.1466 

1.1294 

Sea 
State 

2.0170 

1.3104 

1.2811 

1.2978 

Wind 
Velocity 

1.3815 

1.3590 

1.3415 

1  .3545 

Visibility 

4.5729 

4.6738 

4.7227 

4.6905 

A  t-test  was  performed  to  determine  whether  the 
differences  between  column  means  could  be  attributed  to 
chance  or  to  the  fact  that  the  means  were  actually  unequal. 
(The  F-test  previously  used  for  the  analysis  of  means  for 
the  five  years  would  have  also  been  appropriate  in  this 
case,  but  the  t-test  was  used  to  present  an  alternative 
method.)  The  respective  hypotheses  which  are  being  tested 
are: 


H  ;  x    =  x 
0   5yr    71 


H  :  x    =  x 
0    5yr    72 


H  :  x     =  x 
0   5yr    73 


H  :  The  column  means  are  not  equal 
1 
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Since   the   data   represent  large  random  samples  of  size  N 

and  N  ,  the  sampling  distribution  of  the   statistic   x 

2  5yr 

"x     could  be  approximated  closely  by  the  normal.  The  common 
7x 

variance,  sigma  squared,  was  estimated  from  the   populations 

and   pooled.    The  statistician  realizes  that  the  sample  for 

■£     was  not  strictly  independent  from  the   sample   for   x" 
5yr  7x 

for   any   year  since  each  contained  common  observations  from 

the   area   of   interest.    However    the   effect    of   the 

approximately   300  observations  from  the  area  of  interest  on 

the  more  than  70,000  observations   for   all   the   nationwide 

data   for   any   year  was  minimal.  The  t-ratio  was  calculated 

2 
from  the  following   formula,   where   s    was   the   estimated 


variance: 


As  was  the  case  with  the  F-test,  if  the  computed  value 
of  the  t-ratio  was  greater  than  the  critical  t-statistic, 
the  null  hypothesis  was  not  This  was  a  two-tailed  test  in 
which  the  means  of  the  data  for  the  area  of  interest  could 
be  either  greater  than,  equal,  or  less  than  the  means  of  the 
nationwide  data.  The  results  of  the  calculations  for  the 
area  of  interest  compared  with  the  data  for  the  years  1971 
through  1973  are  tabulated  below: 


27 


COMPARISON  OF  MEANS 
POOLED  DATA  WITH  NATION-WIDE  DATA 


1971 
t-ratio 

1972 
t-ratio 

1973 
t-ratio 

t-sta tistic 
critical 

1  1.  13 

9.74 

11.56 

1  .960 

2.03 

1.92 

2.00 

1  .960 

11,39 

11.45 

11.43 

1  .960 

10.74 

10.35 

10.76 

1  .960 

14.77 

15.39 

14.78 

1  .960 

0.56 

0.99 

0.44 

1.960 

1.98 

2.94 

2.20 

1  .960 

Nature  of 
Distress 

Distance 
Offshore 

Severity  to 
Personnel 

Severity  to 
Property 

Sea  State 


Visibility 

From  these  results,  it  was  inferred  that  the  means  of 
the  values  of  nature  of  distress,  the  severity  to  persons 
and  property,  and  the  sea  state  and  visibility  on  scene  were 
significantly  different  from  those  of  the  rest  of  the  United 
States  at  the  95  percent  confidence  level.  The  differences 
between  the  means  of  degrees  of  severity,  for  example, 
indicate  that  the  danger  to  persons  and  property  were 
greater  in  the  Areata  area  of  interest  than  the  United 
States  as  a  whole.  This  hypothesis  is  supported  by  a 
one-tailed  t-test  at  the  95  percent  confidence  level.  The 
means  of  the  values  for  sea  states  and  visibilities  were 
also  greater. 

Having  compared  the  means  of  the  available  year-to-year 
relationships,  the  next  step  was  to  examine  the  data  of  each 
year  with  the  other  data  of  that  year.  This  was  done  to 
determine  if  there  was  any  statistical  association  among 
elements  of  the  data.  Intuitively,  for  instance,  it  could  be 
postulated  that  the  length  of  the  distressed  vessel  and  the 
distance  offshore  were  associated  since  small  vessels  tend 
to  remain  closer  to  shore  while  larger  vessels   with   better 
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navigational  equipment  tend  to  remain  offshore.  Similiarly , 
it  could  be  postulated  that  there  would  be  a  relationship 
between  weather  conditions  and  severity  to  personnel  or 
property.  If  such  relations  were  discovered,  it  might  have 
been  possible  to  hypothesize  which  factors  were  independent 
(not  relying  upon  the  occurrence  of  any  of  the  others)  and 
which  were  dependent  (occurring  as  a  function  of  some  other 
occurrence) .  Such  independent-dependent  relationships  would 
prove  very  useful  in  understanding  the  SAB  environment  in  an 
area.  An  indication  of  the  strength  of  these  linear 
relationships  was  found  by  examining  the  Pearson 
product-moment  coefficient  of  correlation.  As  the 
coefficient  of  correlation,  r,  approached  the  value  1.,  a 
strong  linear  relationship  is  indicated  while  the  indicated 
degree  of  association  decreased  as  r  approached  zero.  A 
positive  value  of  r  indicated  that  the  values  of  one  of  a 
pair  of  variables  increase  as  the  values  of  the  other 
increase  while  a  negative  value  indicated  the  reverse 
relationship. 

The  coefficients  of  correlation  for  data  can  be 
obtained  by  manual  computation  or  from  any  of  several 
computer  statistical  packages.  The  tables  of  coefficients 
of  correlation  presented  in  Appendix  E  were  obtained  from 
the  SNAP/IEDA  Computing  Package  for  the  IBM  360  [  Ref  14]. 
These  values  are  for  the  elements  of  historical  Search  and 
Rescue  incidents  for  the  Areata  area  of  interest  for  the 
years  1970  through  1974. 

Although  useful  for  indicating  relationships  warranting 
further  study,  the  simple  coefficients  of  correlation  were 
not  suitable  for  determining  the  variation  explained  by  the 
relationships.  To  the  unwary,  the  fact  that  the  value  of  r 
for  one  relationship  was  twice  as  large  as  that  for  another 
would  tend  to  lead  one  to  believe  that  the  relationship  was 
"twice  as  strong."   This  was  not  necessarily  the  case.   From 
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the  tables  of  coefficients  of  correlation  it  was  possible  to 
calculate  multiple  correlation  coefficients  for  the 
variables  by  squaring  each  value  of  r.  The  multiple 
correlation  coefficient,  called  Mult  B,  was  the  maximum 
correlation  between  the  dependent  variable  and  all 
independent  variables.   Its  value  is  represented  so: 

2 
Hult  R   =   r   =   Explained  Variation 

Total  "Variation  "" 

From  the  tables  of  coefficients  of  correlation  for  the 
historic  data  it  was  observed  that  the  only  consistent 
binary  relationships  appear  to  be  between  sea  state  and  wind 
velocity  (r  varied  between  0.46  and  0.86),  between  severity 
to  persons  and  severity  to  property  (r  varied  between  0.23 
and  0.63) ,  between  distance  offshore  and  sea  state  (r 
varied  between  0.26  and  0.70),  and  between  distance  offshore 
and  wind  velocity  (r  varied  between  0.15  and  0.78).  The 
fact  that  the  codes  for  incidents  involving  property  are 
centered  at  the  lower  numbers  v.-hile  those  involving 
personnel  only  (personnel  injury  or  involving  divers)  are  at 
the  higher  end  of  the  scale  was  indicated  by  the  positive 
relationship  between  nature  of  distress  and  severity  to 
personnel  and  the  negative  relationship  between  nature  of 
distress  and  severity  to  property.  Since  consistent  binary 
relationships  had  not  been  discovered  through  the 
coefficients  of  correlation,  conclusions  about  which 
variables  were  dependent  and  which  were  independent  could 
not  be  made.  For  the  purpose  of  the  thesis,  the  most 
critical  variables  are  severity  to  personnel  and  to 
property,  because  greater  severity  implies  the  need  for 
faster  response  by  Coast  Guard  units.  For  this  reason,  the 
severity  to  personnel  and  property  variables  were  assigned 
structural  dependence  for  the  purpose  of  further  analysis. 

It   is  intuitively  obvious  that  none  of  the  independent 
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elements  of  a  SAR  incident  such  as  sea  state,  distance 
offshore,  or  nature  of  distress  act  alone  on  the  dependent 
condition  which  is  severity  to  personnel  or  severity  to 
property.  A  straightforward  extension  of  the  association 
between  two  variables  is  the  association  of  several 
independent  variables  and  the  dependent  variable.  This  is 
termed  multiple  correlation,  and  it  is  also  designated  by 
the  term  Mult  E.  As  in  two-way  multiple  correlations,  Mult  R 
is  the  ratio  of  the  explained  variation  to  the  total 
variation.  This  multiple  correlation  coefficient  for  more 
than  two  variables  could  not  be  found  by  simply  adding  the 
correlation  coefficients  found  above  using  only  two 
variables  because  the  interaction  of  all  independent 
variables  acting  simultaneously  must  be  taken  into  account. 
For  this  analysis,  the  values  of  Mult  R  were  determined 
using  the  BIMED02R  program  from  the  Biomedical  Computer 
Package  [Ref.  2].  (They  could  be  found  using  manual 
computations  or  any  of  the  many  statistical  packages 
available.)  The  strength  of  the  relationships  of  the 
variables  nature  of  distress,  distance  offshore,  sea  state, 
wind  velocity,  visibility,  and  length  of  vessel  were 
determined  separately  for  severity  to  property  and  severity 
to  personnel  for  each  of  the  five  years  data.  The  results 
are  tabulated  below: 
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MULTIPLE  CORRELATION  COEFFICIENTS 
Severity  to  Personnel 


J  9  70 

1211 

121Z 

1211 

1211 

Nature  of 
Distress 

0.  12 

0.07 

0.04 

0.22 

0.  17 

Distance 
Offshore 

0.00 

0.00 

0.00 

0.00 

0.02 

Sea  State 

0.04 

0.  11 

0.  15 

0.01 

0.06 

Wind 
Velocity 

0.01 

0.01 

0.03 

0.00 

0.00 

Visibility 

0.00 

0.0  0 

0.00 

0.01 

0.00 

Length 

0^08 

P-il2 

0.2.02 

<L.oo 

o^oo 

Total 

Explained 
Variation  (%) 

0.25 

0.31 

0.24 

0.25 

0.25 

MULTIPLE  CORRELATION  COEFFICIENTS 
Severity  to  Property 


1219. 

1211 

1212 

1£7_3 

J974 

Nature  of 
Distress 

0.  16 

0.02 

0.03 

0.  12 

0.  12 

Distance 
Offshore 

0.02 

0.00 

0.00 

0.02 

0.00 

Sea  State 

0.02 

0.04 

0.02 

0.00 

0.  11 

Wind 
Velocity 

0.00 

0.00 

0.  00 

0.01 

0.00 

Visibility 

0.  00 

0.00 

0,01 

0.00 

0.01 

Length 

<Ll08 

0^.37 

(K00 

0_s0£_ 

CgOO 

Total 
Explained 
Variation  (%) 

0.28 

0.42 

0.07 

0.  15 

0.23 

As  can  be  seen,  the  explained  variation  of  severity  to 
personnel  and  property  was  relatively  low  in  all  cases.  This 
can  be  interpreted  as  indicating  that  severity  to  personnel 
and  property  were  essentially  random.  Regardless  of  what 
the  weather  conditions,  the  distance  offshore,  the  length  of 
boat  in  distress,  or  even  the  nature  of  distress  are,  the 
severity  to  personnel  or  property  was  just  as  likely  to  be 
none  as  it  was  to  be  severe  in  the  historic  data. 
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The  determination  of  the  simple  correlation 
coefficients,  the  multiple  correlation  coefficients  for  two 
variables,  and  the  multiple  correlation  coefficients  for  all 
independent  variables  indicated  that  the  strengths  of 
relationships  between  the  various  factors  in  a  SAR  incident 
such  as  severity,  weather,  nature  of  distress,  and  distance 
offshore  were  relatively  low  (maximum  possible  is  a 
correlation  coefficient  of  1.0).  This  was  interpreted  as 
meaning  that  these  states  of  nature  could  be  simulated  as 
independent  events.  Their  occurrence  in  any  incident  could 
be  generated  as  an  event  independently  from  all  the  rest  of 
the  factors.  (Independent  in  the  sense  that  they  are 
separate  occurrences.)  They  did  not  depend  upon  any  of  the 
other  events  listed.  The  truth  of  this  assumption  was  borne 
out  in  the  validation  of  the  simulation  model  when  it  was 
ascertained  that  the  distributions  of  the  factors  of  a  SAR 
incident  were  the  same  for  simulation-generated  data  and 
historic  data. 

Another  data  analysis  which  might  prove  interesting  in 
the  future  is  the  investigation  of  the  relationships  between 
number  of  boats  registered,  their  lengths,  and  the  number  of 
boats  assisted  by  the  Coast  Guard.  If  correlations  were 
found  between  the  number  and  description  of  boats  registered 
and  documented  in  an  area  and  the  number  of  Coast  Guard 
cases,  it  might  be  possible  to  derive  an  equation  which 
would  predict  the  Coast  Guard  workload  given  the  boat 
registrations  for  the  current  or  some  prior  period.  The 
workload  may  be  a  function  of  lengths  of  vessels, 
propulsion,  usage,  or  age.  Such  a  study  would  require  the 
derivation  of  a  set  of  regression  equations.  To  conduct 
such  an  analysis  would  require  detailed  vessel  registration 
data  in  addition  to  the  Coast  Guard  SAR  data  for  a 
particular  area,  but  the  benefit  may  be  a  consistent, 
reliable  predictor  of  case  loads  and  the  resulting  need  for 
Coast  Guard  facilities. 
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Another  aspect  of  the  historical  data  was  investigated. 
On  the  Assistance  Report,  the  reporting  Coast  Guard  unit  is 
required  to  report  the  severity  to  personnel  and  property. 
If  the  severity  codes  are  higher  for  any  given  area,  it 
could  indicate  that  more  rapid  response  is  necessary.  This 
procedure  was  described  previously.  On  the  Assistance 
Report,  the  severity  to  persons  or  property  is  defined  and 
coded  as  follows10 

0  None  No  personnel/property  involved. 

1  Snail  No  immediate  or  foreseeable 
danger. 

2  Moderate       Some  danger  that  personnel  or 
property  might  be  lost. 

3  Severe         Personnel/property  was  in  danger 
of  being  lost  or  was  actually  lost. 

9     Unknown 

For  reasons  previously  described,  all  cases  with  code 
values  of  nine  were  sorted  out  and  not  used  in  the  data 
analysis. 

In  another  location  on  the  Assistance  Report  form,  the 
reporting  unit  records  what  may  be  described  as  the  result 
of  the  actions  of  the  Coast  Guard.  In  this  section  the 
number  of  lives  saved,  number  of  lives  lost,  number  of 
persons  otherwise  assisted,  and  the  value  of  property 
assisted  are  reported. 

In  an  effort  to  derive  a  cost-benefit  ratio  for  its  SAR 
program,  the  Coast  Guard  has  developed  a  method  in  which 
each  of  these  outcomes  is  assigned  a  dollar  value   and   used 
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in  a  formula  to  determine  an  effectiveness  ratio.  The 
effectiveness  ratio  for  SAR  is  value  of  benefits  divided  by 
the  sum  of  the  value  of  benefits  and  the  value  of  lcsses. 
The  only  losses  reported  are  lives  lost.  The  same  figures 
are  used  to  derive  a  cost-effectiveness  ratio,  which  is  the 
value  of  benefits  divided  by  costs  of  the  SAR  program.  ll 

It  is  apparent  that  the  adjectives  for  describing 
severity  are  open  to  subjective  interpretation.  An  analysis 
was  made  to  determine  if  the  reporting  units  were  consistent 
in  their  reporting  of  degrees  of  severity.  This  is,  in 
effect,  an  evaluation  of  the  quality  of  the  data.  The  only 
manner  in  which  this  could  be  done  was  to  compare  the  number 
of  times  which  severity  to  personnel  was  given  as  moderate 
or  severe  (indicating  that  outcomes  of  loss  or  saving  of 
life  might  be  expected)  with  the  coincidences  of  reports  of 
actual  loss  or  saving  of  life. 

In  the  years  1970  through  1974,  Coast  Guard  resources 
were  involved  in  assisting  2965  persons  in  the  V468  cases 
used  as  the  data  base.  In  these  cases,  38  persons  lost  their 
lives  despite  Coast  Guard  efforts  in  29  cases.  In  56  other 
cases,  113  persons  were  reported  saved,  including  one 
incident  in  which  fourteen  lives  were  reported  saved.  The 
greatest  number  of  lives  lost  in  any  one  case  was  three- 
The  balance,  2814  persons,  were  reported  assisted  in  1200 
cases.  This  data  was  adjusted  for  three  cases  in  which  some 
lives  were  reported  helped  and  some  were  reported  saved,  two 
cases  where  some  lives  were  lost  and  the  rest  of  the  persons 
assisted,  and  one  case  where  one  life  was  saved  and  one 
lost.  Parenthetically,  one  hundred  cases  were  coded  to 
indicate  only  personnel  in  distress  and  fifty-six  were  coded 
as  only  property  in  distress. 
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The  distribution  of  reported  severity  to  personnel  in 
the  five  years  was: 

None        Small      Moderate      Severe 

Number     56  949         305         158 

Percent   3.8%        64.6%       20.8%        10.8% 

The  distribution  of  reported  outcomes  of  Coast  Guard 
assistance  was: 

Lives       Lives      Persons       No 
Lost        Saved     Assisted     Persons 
Cases       Cases       Cases     Involved 

Number        29  56  1200  183 

Percent       2.0$        3.8%       81.7%  12.5% 

No  reconcililiation  was  possible  between  the  numbers  183, 
which  indicated  no  persons  assisted,  and  the  number  56, 
which  indicated  that  no  personnel  were  in  distress.  The  only 
logical  explanation  for  some  of  these  discrepancies  is  that 
a  Coast  Guard  unit  was  dispatched  to  assist  but  aborted  the 
mission  before  reaching  the  scene  of  the  incident. 

Using  this  data,  it  was  possible  to  attempt  to 
determine  how  many  cases  reported  as  having  the  possibility 
of  loss  of  life  were  actually  that  severe—as  reported  by 
the  Coast  Guard  units.  The  total  number  of  cases  where  lives 
were  lost  (29)  plus  the  number  of  cases  where  lives  were 
reported  saved  (56)  divided  by  the  number  of  cases  reported 
as  presenting  the  greatest  severity  to  personnel  (158)  would 
indicate  that  approximately  54  percent  of  those  incidents 
reported  as  being  of  great  severity  to  persons  actually  were 
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so.  In  other  words,  rather  than  10.8  percent  of  all  the 
historical  cases  having  a  severity  state  to  personnel  as 
severe  (code  value  3),  only  about  5.8  percent  of  the  cases 
were  actually  that  severe.  It  was  interesting  to  note  that 
of  all  the  cases  reported  as  moderate  severity  to  personnel, 
in  no  case  was  a  life  actually  reported  lost  or  saved.  It 
appeared  that  persons  filling  out  the  Assistance  Reports 
were  generous  in  assigning  the  degree  of  severity  of  the 
cases  they  reported  and  might  have  been  assigning  a  high 
severity  code  in  more  cases  than  warranted. 

This  disparity  could  be  isolated  only  by  close 
examination  cf  the  data  reported  on  the  Assistance  Reports. 
In  the  blocks  where  severity  to  property  and  value  of 
property  assisted  are  reported,  no  such  analysis  can  be 
made.  Nor  is  there  any  method  of  checking  the  quality  of 
the  figures  given  for  lives  saved.  The  person  filling  out 
the  Assistance  Report  must  enter  the  estimated  value  of  the 
property  assisted  if  no  accurate  value  is  known,  and  there 
is  no  way  to  report  the  value  of  property  lost.  The 
efficiency  ratio  is  thereby  biased  and  there  is  no  way  to 
verify  the  figures  used  to  calculate  it.  This  is  an  area 
where  further  study  might  prove  beneficial  in  order  to 
determine  more  reliable  and  realistic  measures  of 
performance  of  the  SAR  mission. 
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IV*   COMPUTER  SIMULATION  OF  THE  ENVIRONMENT 

After  the  data  for  the  area  of  interest  has  been 
examined,  the  analyst  may  desire  to  develop  a  computer 
simulation  model  of  the  environment  within  that  area. 
Developing  a  logical  portrayal  of  the  SAR  environment  for 
such  a  simulation  is  a  means  of  closely  studying  r.he 
relationships  of  all  the  pertinent  factors  in  that 
environment  and  can  lead  to  a  better  understanding  of  the 
complexities  of  the  environment. 

A  second  use  for  a  simulation  model  to  to  generate  more 
data.  Only  data  for  five  years  was  available  for  the  Areata 
area  of  interest,  and  this  was  deemed  insufficient  for  the 
purposes  of  analyzing  the  need  for  a  new  air  station. 

Another  application  for  such  a  computer  simulation 
model  is  to  manipulate  it  in  order  to  test  the  Coast  Guard 
facilities  in  an  area  against  differing  environments.  In 
the  Areata  area  of  interest,  for  example^  there  have  been 
very  few  incidents  reported  inland.  This  is  because  or! y 
water-borne  units  were  stationed  in  the  area  and,  except  for 
a  few  isolated  cases,  the  Coast  Guard  has  no  means  of 
responding  to  inland  SAR  incidents.  Analyzing  the  past,  data, 
one  could  conclude  there  were  no  incidents  requiring  Coast 
Guard  assistance  which  were  not  on  the  water.  It  is  not 
true  that  what  was  not  recorded  never  happened,  of  course, 
and  the  means  of  inserting  such  land  cases  as  aircraft 
crashes,  lost  hunters,  and  automobile  crashes  is  presented 
through  a  computer  simulation  model.  In  another  example,  if 
some  relationship  between  number  of  distress  incidents  and 
the  number  of  boats  registered  in  an  area  were  found,  future 
case  loads  could  be  simulated  by  causing  more  incidents  to 
occur  in  the  computer  simulation  model. 
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The  purpose  of  the  model  for  the  area  of  interest,  as 
written  for  this  thesis,  was  to  approximate  the  environment 
in  which  the  Coast  Guard  forces  in  that  area  will  operate. 
The  parameters  of  the  model  were  based  upon  historical  data 
which  had  been  gathered  during  the  course  of  the  normal 
Search  And  Rescue  operations  of  the  Coast  Guard sin  that 
area.  Once  the  model  was  functioning  properly  and  its 
output  was  validated  against  the  historic  data,  it  could  be 
used  to  generate  more  data  to  use  investigate  the  need  for 
the  Areata  Air  Station. 

The  historic  data,  as  originally  received,  had  to  be 
sorted  in  order  to  obtain  only  those  cases  which  had 
occurred  uithin  the  area  of  interest.  The  information  which 
was  needed  to  model  the  environment  in  the  simulation  was 
that  pertaining  to  the  weather  reported  during  the  distress 
incidents,  the  dates  and  times  of  the  incidents,  the  types 
of  distress,  and  the  positions  of  the  distress.  These  items 
are  the  elements  of  a  search  and  rescue  incident  which  occur 
randomly.  All  further  information  reported  on  the 
Assistance  Report,  such  as  distance  to  scene  of  the 
distress,  action  taken,  and  time  the  assisting  resource  gets 
underway  or  arrives  on  scene  are  raore-or-less  a  result  of 
the  random  factors.  As  indicated  previously,  the  random 
elements  were  also  termed  "states  of  nature"  in  this  thesis. 

The  data  which  was  made  available  tor  the  model  of  the 
Areata  area  was  for  the  years  1S70  through  1974.  Data  for 
the  position  of  distress  and  the  date  for  the  incidents 
occurring  in  the  years  1970  and  1974  was  missing  and  could 
not  be  recovered  in  time  to  be  used  in  the  thesis.  The 
historic  data  was  sorted  and  examined  statistically  to  make 
it  usable  for  the  modelling  of  the  environment  of  the 
proposed  Areata  Air  Station.  The  product  of  this  analysis 
was  a  set  of  cumulative  distributions  representing  the 
frequencies   of  the  states  of  nature  which  had  been  observed 


39 


by  the  Coast  Guard  units  during  search  and  rescue  operations 
for  the  years  1970  through  1974.  When  a  model  of  the 
environment  was  constructed,  these  historic  cumulative 
distributions  were  relied  upon  to  ensure  that  the  model 
generated  states  of  nature  which  corresponded  to  the 
historic  facts.  Since  it  was  proved  in  the  previous  chapter 
that  the  elements  of  the  SAH  environment  occurred  randomly, 
the  model  was  designed  to  generate  various  elements  in  any 
distress  incident  randomly.  Within  the  computer  model,  a 
random  number  was  generated  and  compared  to  the  historic 
cumulative  distribution  of  the  states  of  nature.  Most 
computer  facilities  include  a  program  in  their  software 
package  which  will  generate  uniformly  distributed  random 
numbers.  At  the  Naval  Postgraduate  School,  the  program 
which  was  used  was  called  RANDOH  [Ref.  16].  This  program 
has  been  statistically  proved  as  being  a  very  good 
random-number  generator  and  for  the  purposes  of  this  thesis, 
the  numbers  obtained  from  that  program  can  be  assumed  to  be 
truly  "random."  This  assumption  was  supported  in  tests  to 
be  described  later. 

In  order  to  derive  the  code  value  for  any  element  of 
the  environment  in  the  simulation,  the  computer  program 
first  called  for  a  random  number.  This  random  number  was 
used  in  combination  with  an  array  stored  in  the  computer 
which  contained  the  cumulative  historic  distributions  of  all 
the  permissible  codes  for  that  element.  In  effect,  the 
array  is  a  table.  The  random  number  is  compared  with  all 
values  in  the  table  in  order  to  determine  which  code  value 
corresponds  to  the  random  number.  When  the  program  finds 
the  range  of  values  which  contains  the  value  of  the  random 
number,  it  determines  what  the  code  value  is  and  returns 
that  value  to  the  simulation  program  to  represent  the  state 
of  nature. 

The   program  for  the  computer  simulation  was  written  in 
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the  FORTRAN  IV  language.  Its  scope  was  limited  to  model 
only  the  environment  and  to  generate  SAR  cases  within  the 
specified  area  of  interest.  This  meant  that  actions  of  the 
Coast  Guard  units,  such  as  towing  disabled  vessels  to  port 
or  medevacing  patients  to  hospitals,  were  not  modelled 
because  these  were  outside  the  scope  of  this  analysis.  If 
the  simulation  generated  a  case  in  which  a  boat  was  lost  or 
disoriented  at  a  specified  position,  the  time  spent  by  a 
Coast  Guard  unit  in  searching  for  that  boat  was  not 
modelled,  although  this  and  other  factors  could  be  added  to 
the  computer  program. 

Another  factor  in  the  operations  of  the  Coast  Guard 
which  was  not  modelled  was  the  readiness  status  of  the  SAR 
facilities.  It  was  assumed  that  if  an  incident  occurred,  a 
Coast  Guard  unit  could  get  underway  immediately  to  respond 
without  recalling  the  crew  or  getting  ready  for  sea.  This 
assumption  means  that  when  analysing  response  times  for  the 
Coast  Guard  units,  the  calculated  time  to  respond  is  the 
minimum  time  to  arrive  on  scene.  The  computer  program  used 
to  simulate  the  environment  was  included  in  a  later  section 
of  this  thesis. 

The  first  step  of  the  computer  model  was  to  establish 
the  dimensions  of  the  arrays  and  read  into  these  arrays  the 
historic  distributions.  Storage  arrays  were  designated 
which  corresponded  to  the  sea  state,  wind  velocity, 
visibility,  latitude,  longitude,  minutes  of  latitude  arid 
longitude,  nature  of  distress,  severity  to  personnel  and 
property,  and  time  of  distress.  The  data  which  was  read  into 
these  arrays  was  used  at  various  places  in  the  computer 
program  to  generate  a  distress  incident. 

The  main  computer  program  was  used  to  advance  the 
calendar  through  the  model  year  and  to  call  the  subroutines 
that  were  required.   For  the  purposes  of  the  simulation,  the 
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Gregorian  year  was  "standardized"  in  order  to  eliminate  the 
confusion  and  difficulties  associated  with  different  days  of 
the  week  falling  on  different  dates  in  different  years. 
Experience  has  shown  that  the  incidence  of  SAR  cases  is 
highly  cyclical  and  dependent  upon  the  season,  the  month, 
and  the  day  of  the  week.  Each  model  year  consisted  of  336 
days,  or  12  months  of  28  days  each.  The  first  day  of  each 
month  was  a  Sunday,  and  the  same  day  of  the  week  fell  on  any 
given  date  regardless  of  the  month  or  the  year.  The 
historic  day-of-week  data  of  SAR  incidents  was  analysed  for 
the  years  1971,  1972,  and  1973  and  the  experience  of  these 
three  years  was  converted  to  the  model  year.  In  order  to 
determine  the  average  number  of  incidents  which  occurred  on 
the  first  Sunday  of  February,  for  instance,  the  number  of 
cases  reported  in  the  Areata  area  for  the  first  Sunday  was 
averaged  and  assigned  to  the  storage  location  for  February 
the  first  in  the  model  year.  For  the  months  which  have  more 
than  28  days,  some  of  the  historic  cases  were  lost  because 
there  was  no  corresponding  model  date  in  which  to  place 
them.  The  number  of  times  this  happened  during  the  three 
historic  years  was  so  small  that  the  effect  of  losing  these 
incidents  was  considered  insignificant.  In  order  to 
increase  the  total  number  of  cases  in  the  area  of  interest, 
the  user  of  the  computer  simulation  could  increase  the  daily 
average.  In  a  like  manner,  the  workload  for  any  given  day 
could  be  adjusted. 

In  order  to  determine  the  number  of  incidents  which 
were  to  occur  on  a  model  day,  a  subroutine  was  called  by  the 
main  program.  The  number  of  SAR  incidents  which  have 
occurred  on  any  day  in  the  past  fulfilled  the  classical 
Poisson  conditions;  the  number  of  opportunities  for  an 
incident  was  infinite,  the  probability  of  ail  incident 
occurring  for  any  given  opportunity  approached  zero,  yet 
incidents  happened.  The  average  number  which  happened  on  any 
day   is   known.    The   subroutine   used  the  current  date  and 
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month  to  query  the  array  in  which  the  average  number  of 
incidents  for  each  day  was  stored.  It  used  a  Poisson 
distribution  to  determine  the  number  of  cases  which  were  to 
occur  during  the  current  year  on  that  dare  and  returned  that 
number  to  the  main  program.  If  no  cases  were  to  occur  that 
day,  the  main  program  skipped  any  further  steps,  incremented 
the  date  by  one  day,  and  began  the  iteration  again., 

If  a  SAR  incident  was  to  occur  on  any  day,  the  main 
program  called  a  subroutine  which  determined  the  hours  *;hich 
the  incident  or  incidents  were  to  happen.  If  more  than  one 
incident  was  scheduled  to  occur,  the  subroutine  generated  a 
corresponding  number  of  times  of  day,  arranged  them  into 
ascending  order,  and  returned  the  first  to  the  main  program 
to  be  used  as  the  time  of  the  first  incident.  The  main 
program  controlled  the  simulation  and  generated  the  required 
number  of  cases  prescribed  in  the  earlier  subroutine.  As 
each  incident  was  generated.,  the  main  program  referred  to 
the  time  subroutine  to  obtain  the  time  of  the  current 
incident.  Each  time  it  was  called,  the  subroutine  returned 
the   hours  and  minutes  (clock  time)  of  the  current  incident. 

Having  determined  that  an  incident  was  to  occur  on  any 
day  and  the  time  of  that  incident,  the  main  program  next 
called  a  subroutine  which  determined  nature  of  distress 
representing  the  incident  and  the  severity  involved.  The 
Coast  Guard  has  responded  to  virtually  every  type  of 
distress  that  can  befall  man  or  his  possessions,  but  in  the 
Areata  area  relatively  few  of  them  occurred  during  the  years 
1970  through  1974.  It  was  assumed  that  this  sample  was 
sufficiently  representative  so  that  the  model  would  not  have 
to  generate  types  of  distress  which  did  not  occur  during 
these  years.  Each  of  the  types  of  distress  which  were 
reported  during  these  years  also  had  a  degree  of  severity  to 
either  persons  and/or  property  associated  with  them.  Both 
the    distributions    of   the   types   of   distress   and   the 


43 


distributions  of  the  degrees  of  severity  to  persons  and 
property  were  stored  in  an  array  in  the  computer  model.  The 
subroutine  which  determined  the  nature  of  distress  generated 
a  random  number,  compared  it  with  the  cumulative 
distributions  in  the  nature-of-distress  array,  and  derived 
the  nature  of  distress  for  the  incident  it  was  in  the 
process  of  creating.  Once  the  nature  of  distress  was 
determined,  the  subroutine  could  then  determine  the  degree 
of  severity  to  personnel  or  property  in  a  similiar  manner. 

The  main  program  next  determined  the  position  of  the 
incident.  The  historic  distribution  of  SAR  incidents  was 
stored  in  an  array  in  the  computer  model,  but  since  the 
overwhelming  number  of  cases  had  historically  been  handled 
by  coast  Guard  surface  units,  this  distribution  was  not 
truly  representative  of  the  positions  of  incidents  which 
would  be  handled  by  the  Areata  air  station.  The  area  of  the 
Pacific  Ocean  to  be  serviced  by  the  proposed  air  station  was 
between  40-00N  and  42-00N  Latitude  and  roughly  124-10W  and 
125^00H.  During  the  years  1971  through  1973  approximately 
94.4  percent  of  the  SAR  incidents  occurred  within  this  area. 
In  the  model  this  area  was  reduced  to  a  grid,  each  element 
of  which  represented  five  minutes  of  latitude  (five  miles) 
and  five  minutes  of  longitude  (approximately  3.7  miles).  In 
the  model,  each  grid  was  assigned  a  number  representing  the 
historic  cumulative  distribution  of  cases  which  had  occurred 
within  its  confines  in  three  years.  The  subroutine  which 
determined  the  position  of  distress  generated  a  random 
number  and  compared  it  with  the  cumulative  numbers  stored  in 
the  grid  array  in  order  to  determine  the  latitude  and 
longitude  of  the  distress  incident  being  generated  by  the 
main  program.  If  the  random  number  *as  greater  than  94.1, 
it  meant  that  the  incident  occurred  outside  the  confines  of 
the  288  grids,  or,  in  other  words,  that  it  had  either 
occurred  over  the  land  (less  than  124-1 Oa)  or  farther  nest 
than  125-00H.   If  this   was   the   case   and   the   nature   of 
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distress  indicated  that  the  incident  had  occurred  on  the 
water  (such  as  vessel  disabled,  or  fire  on  a  vessel) ,  the 
subroutine  generated  a  Latitude  for  the  distress  and  a 
distance  offshore  greater  than  125-00W.  An  incident  in  such 
a  location  normally  could  not  be  handled  by  the 
single-engined  HH-52  helicopter  which  would  be  stationed  at 
the  proposed  air  station.  If  the  seriousness  of  the 
situation  warranted,  a  long-range  HH-3F  helicopter  would  be 
flown  to  Areata  from  another  air  station  to  perform  the 
mission.  If  the  nature  of  distress  indicated  that  the 
incident  occurred  over  the  land  (such  as  personnel  injury  or 
land  vehicle  crash)  the  subroutine  generated  a  bearing  and 
range  from  the  proposed  air  station  to  the  scene  of  distress 
to  represent  the  position  of  the  incident.  The  bearing  and 
range  were  generated  randomly  because  of  lack  of  any 
historic  data  to  indicate  the  distribution  of  types  and 
positions  of  land  cases. 

The  logical  result  of  knowing  the  position  of  the 
distress  incident  is  the  determination  of  the  distance  to 
that  position  from  the  various  Coast  Guard  facilities  in  the 
area  of  interest*  This  computation,  as  well  as  that  which 
determined  the  time  for  a  Coast  Guard  assisting  resource  to 
arrive  on  scene  of  the  distress,  was  performed  by  another 
subroutine  which  was  called  by  the  main  program  after  the 
position  of  the  incident  had  been  determined.  This 
subroutine  was  used  to  determine  the  minuraum  response  time 
for  a  Coast  Guard  unit  to  arrive  on  scene  of  the  generated 
incident. 

The  final  random  state  of  nature  in  the  environment  of 
a  SAR  incident  is  the  weather.  The  weather  conditions 
encountered  by  the  Coast  Guard  facilities  during  the 
prosecution  of  their  search-and- rescue  operations  were  in 
the  five  years  data..  For  each  incident,  the  wind  velocity, 
the  sea  state,   and   the   visibility   was   reported   on   the 
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Assistance  Report.  The  main  program  called  a  subroutine 
which  generated  these  weather  conditions  according  to  the 
historic   distributions  stored  in  an  array  within,  the  node!. 

Although  pertinent  to  the  operations  of  the  Coast 
Guard,  the  length  of  the  distressed  vessel  was  not  generated 
because  of  the  implied  necessity  of  analysing  position  and 
weather  and  their  combined  effects  on  length.  It  would  be 
necessary  to  determine  if  the  historic  cases  involving 
smaller  vessels  tend  to  center  in  any  specific  area  and 
generate  lengths  that  corresponded  to  these  restrictions. 
For  the  purpose  of  this  thesis,  it  was  assumed  that  any  of 
the  Coast  Guard  Cutters  in  the  Areata  area  of  interest  was 
capable  of  responding  to  any  of  the  incidents  generated. 
The  purpose  of  the  thesis  was  to  investigate  the  need  for 
rapid  response  and  not  for  determining  need  for  towing 
capability. 

Insofar  as  possible,  internal  consistency  was  ensured 
in  the  simulation  model.  If  the  nature  of  distress 
generated  was  one  which  normally  occurs  only  over  land,  a 
position  inland  was  generated  and  the  sea  state  code  was 
made  to  correspond  to  an  inland  case.  If  the  nature  of 
distress  was  vessel  aground,  the  simulation  of  position  was 
constrained  to  one  adjacent  to  the  coastline.  Enforcing 
such  internal  consistencies  would  tend  to  cause  small 
distortions  in  the  distributions  of  the  environmental 
elements,  but  the  degree  of  this  distortion  could  not  be 
discerned. 

While  it  was  possible  to  provide  for  internal 
consistencies  within  any  one  case.,  the  model  is  raemoryless 
and  consistencies  from  case  to  case  could  not  be  ensured. 
This  meant  that  if  one  case  were  generated  in  any  position 
with  wind  and  seas  calm,  another  case  could  be  generated  in 
exactly   the   same   position   one  hour  later  with  conditions 
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such  as  would  be  experienced  in  a  full  gale.  Neither  the 
weather  from  the  previous  case  nor  the  fact  that  a  Coast 
Guard  unit  was  very  close  to  the  scene  of  distress  was 
"remembered"  in  the  model.  This  was  considered  to  be  no 
great  problem  since,  again,  the  purpose  of  the  simulation 
was  not  to  model  the  actions  of  the  Coast  Guard  units  but 
instead  to  model  the  environment.  The  results  of  modelling 
the  environment  were  to  be  used  to  generate  statistics  and 
not  to  recreate  reality.  Output  from  the  simulation  program 
was   included  in  the  computer  output  section  of  this  thesis. 

In  1971  the  National  Bureau  of  Standards,  in  a  joint 
effort  with  the  Coast  Guard,  developed  a  Search  and  Rescue 
Simulation  Model  (SARSIM) to  assist  Coast  Guard  management  in 
planning  for  its  SAR  mission  [Ref.7].  The  model  is  highly 
detailed  and  quite  wide  in  scope.  The  model  developed  for 
this  thesis  is  in  no  way  meant  to  replace  that  resource.,  but 
SARSIM  was  not  available  to  the  student  for  use.  As 
previously  stated,  one  of  the  purposes  of  developing  a 
computer  simulation  model  is  to  gain  greater  insight  into 
the  SAR  environment  for  an  area.  This  insight  might  not  be 
gained  by  using  SARSIM.  It  should  also  be  noted  that  the 
run  time  of  one  year's  simulation  on  the  IBM  360/67  at  the 
Naval  Postgraduate  School  is  approximately  30  seconds,  which 
translates  into  an  estimated  computing  cost  of  about  three 
dollars. 
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v#   VALIDATION  OF  THE  COMPUTES  WOP EL 

The  computer  model  of  the  environment  was  not  useful 
until  it  could  be  proved  that  it  was  generating  conditions 
that  represented  the  actual  environment.  In  order  to 
validate  the  model  it  was  necessary  to  check  its  internal 
consistency  and  its  agreement  with  the  historic  environment 
of  the  Areata  area.  The  first  test  was  relatively  simple; 
the  computer  output  was  examined  to  ensure  chat  such 
inconsistencies  as  divers  in  distress  inland  or  vessels 
aground  ten  niles  offshore  were  not  being  generated. 

Testing  the  consistency  of  the  output  of  the  simulation 
of  the  environment  from  the  model  with  the  actual  historic 
environmental  conditions  required  using  the  same  tests  that 
were  used  in  the  analysis  of  the  original  data.  The  program 
for  the  computer  simulation  was  written  so  that  the 
distribution  of  events  in  the  model  could  never  exceed  the 
historic  distributions,  but  it  was  necessary  to  test  whether 
the  distributions  of  the  output  conformed  to  these  historic 
distributions. 

There  were  two  ways  in  which  an  inconsistency  could 
occur.  First,  the  process  of  enforcing  internal  consistency 
within  a  case  could  alter  the  resulting  distributions  of 
frequencies.  The  nature  of  distress,  for  instance,  was 
generated  randomly  and  the  on-sceno  Heather  was  also 
generated  randomly.  In  each  case,  the  resulting  state  of 
nature  was  derived  by  comparing  the  generated  random  number 
with  a  table  of  cumulative  distributions  of  historic  data. 
The  distribution  of  simulated  events  within  the  model  could 
never  exceed  the  historic  distribution  if  this  process  were 
unaltered,  since  both  states  were  derived  independently. 
When  internal  consistency  was  enforced,  however,  the 
independent  derivation  of  states  was  lost.  Randomly 
generating   a   nature   of   distress  which  must  have  occurred 
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inland  forced  a.  sea  state  code  which  indicated  a  land 
incident.  It  also  forced  an  inland  position  of  distress.  A 
large  number  of  such  events  could  alter  the  resulting 
distribution  of  the  dependent  states  such  as  the  position 
and  sea  state  code. 

The  second  possibility  for  error  in  the  output  of  the 
computer  model  was  an  improperly-operating  random  number 
generator.  As  indicated,  the  states  of  nature  depended  upon 
the  generation  of  a  random  number  which  was  compared  with  a 
table  of  distributions  to  derive  the  correct  simulated 
state.  If  the  random  number  generater  consistently  produced 
a  high  or  low  number,  the  resulting  distribution  of  states 
of  nature  would  be  significantly  biased.  The  only  means 
available  to  ensure  that  these  errors  had  not  occurred  was 
to  statistically  compare  the  data  generated  by  the 
simulation  model  with  the  historic  data  which  had  served  as 
the  basis  for  the  model. 

In  order  to  validate  the  computer  model,  three  runs 
representing  three  years  were  made  using  different  seed 
values  for  the  random  number  generator.  The  random  number 
generator  uses  the  seed  value  to  compute  subsequent  numbers. 
When  the  random  number  generator  is  given  the  same  seed 
value,  it  will  produce  an  identical  stream  of  random 
numbers,  Changing  the  seed  value  causes  the  random  number 
generator  to  compute  an  entirely  new  set  of  numbers.  As  was 
done  with  the  historic  data,  the  descriptive  statistics  of 
the  generated  data  were  obtained:.  Histograms  comparing  the 
pooled  historic  data  and  generated  data  were  included  in 
Appendix  C.  The  means  of  the  generated  elements  of  the 
environment  are  given  in  the  following  table.  The  means  of 
the  five-year  pooled  data  are  given  because  these  pooled 
distributions  were  used  in  the  model  to  derive  the  various 
states  of  nature. 
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TABLE  OF  MEANS 
SIMULATION  GENERATED  DATA 


Ban  1 

Run  2 

Run  3 

5YR 
Pooled 

Nature  of 
Distress 

8.3447 

11  .5581 

7.7428 

7.8  583 

Personnel 
Severity 

1.3826 

1.5356 

1.3333 

1.3849 

Property 
Severity 

1 .4280 

1.3521 

1.4203 

1.3542 

Sea  State 

2.0341 

2.1236 

2. 1051 

2.0170 

Hind 
Velocity 

1 .1667 

1.  1985 

1.2319 

1.3815 

Visibility 

4.4811 

4c3970 

4.5652 

4.572  9 

Time  of  Day 

12.7273 

12.67  79 

12. 1848 

12.6679 

Month 

7. 1705 

7.0150 

6.9855 

7.0  342 

The  value  for  distance  offshore  is  not  given  because 
while  historically  the  greatest  number  of  cases  (73  percent) 
occur  within  three  miles  of  the  coastline,  this  percentage 
could  not  be  duplicated  in  the  computer  model  because 
irregularities  in  the  coastline  were  not  duplicated. 
Furthermore,  the  grid  size  for  the  generated  position  of 
distress  was  greater  than  three  miles  on  any  side.  As  a 
result,  the  smallest  distance  offshore  which  could  be 
generated  was  more  than  three  miles.  Another  data  element 
not  listed  is  length  of  distressed  vessel,  which  was  not 
generated  in  the  simulation. 

As  was  done  with  the  historic  data,  a  one-way  analysis 
of  variance  was  performed  comparing  each  of  the  runs  with 
the  historic  data  for  each  of  the  five  years.  The 
hypothesis  being  tested  was: 
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H  :  u.      =  u.    =  u.    =  u.     =  u.    =  u. 
0     1,2,3      70      71      72       73      74 


H  :  The  means  are  not  equal 
1 


These  tests  were  performed  for  each  element  in  the 
table  above  using  the  BIMED01V  program  [Ref  2]  which  gave  an 
F-ratio  as  described  in  the  data  analysis  section.  If  the 
F-ratio  was  greater  than  the  critical  F-statistic  at  least 
one  of  the  column  means  was  significantly  different  at  the 
95  percent  confidence  level  than  the  rest  of  the  means.  The 
results  of  the  test  are  tabulated  below: 


ANALYSIS  OF  MEANS 
SIMULATION-GENERATED  DATA 


P- statistic 

Critical 


Nature  of        4.08       2.94       2.99         2.21 
Distress 

Personnel        9.49       11.28       10.02         2.21 
Severity 

Property         9.82        9.38       9.85         2.21 
Severity 

Sea  State        6.34        6.68       6.58         2.21 

Wind  10.32       10.10        9.73         2.21 

Velocity 

Visibility        1.48        1.73        1.39         2.21 

Time  of  Day       4.85       4.68       8.49         2.60 
Month  12.43       12.23       12.38  2.60 

The  result  of  the  tests  led  to  the  conclusion  that  ail 
the  means  of  the  simulated  data  differed  significantly  from 
the  means  of  the  pooled  data.  In  all  the  cases  where  the 
F-ratio  exceeded  the  critical  F-statisticf  a  pair-wise 
comparison  of  column  means  by  estimation  of  linear  contrasts 
was  performed.  This  was  done  in  order  to  determine  which  of 
the  column  means  were  different,  resulting  in  the   rejection 
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of  the  hypothesis  of  equality.   The  graphic  depiction  of  the 
results  of  these  tests  are  included  in  Appendix  D. 

A  summary  of  the  results  of  these  tests  is  that,  except 
in  a  few  cases,  the  mean  of  the  generated  data  did  not 
differ  significantly  from  the  means  of  the  historic  data. 
Generally,  any  time  there  was  a  significant  difference 
between  simulated-data  means  and  historical-data  means,  the 
historical-data  mean  had  been  found  to  be  significantly 
different  from  all  the  historical  means  in  the  tests 
performed  in  the  original  data  analysis.  This  result  tended 
to  validate  the  presumption  that  the  output  of  the  model 
conformed  to  the  actual  environment. 

In  order  to  test  the  effect  of  forcing  internal 
consistency  within  the  model  and  the  operation  of  the  random 
number  generator,  the  data  for  each  ran  was  compared  with 
the  pooled  data  from  which  it  was  supposedly  derived.  For 
this  comparison  the  t-test  was  performed.  Again,  the 
hypothesis  that  the  row  mean  for  the  generated  data  was  not 
significantly  different  from  the  mean  of  the  five-year 
pooled  data  was  tested  at  the  95  percent  confidence  level. 
A  table  of  the  results  of  that  test  is  presented  below: 
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ANALYSIS  CF  MEANS 
GENERATED  DATA  VERSUS  POOLED  DATA 


t-ratio 
Run  J_ 

t-ratio 
Run  2 

0.35 

2.  13 

0.04 

2.86 

1.36 

0.04 

0.  13 

0„90 

3.27 

2.85 

0.78 

1.42 

0.  19 

0.03 

Nature  of        0.35        2.13        0.09        2.015 
Distress 

Personnel        0.04        2.86        1.11        2.015 
Severity 

Property  1.36        0.04        1.28        2.015 

Severity 


t-ratio 
Run  3 

0. 

09 

1. 

11 

1. 

28 

0. 

7  4 

2. 

14 

0. 

07 

1. 

77 

t--sta  tistic 
Critical 


Sea  State        0.13        0„90       0.74        2.015 

Hind  3.27        2.85       2.14        2.015 

Velocity 

Visibility        0.78        1.42       0.07        2.015 

Time  of  Day       0.19        0.03        1.77        2.3  53 

Month  0.81        0.16        0.30        2.353 

In  most  cases  the  differences  between  the  simulated 
data  and  the  pooled  five-year  actual  data  are  not 
significant  at  the  95  percent  confidence  level.  In  some 
instances  the  t-ratio  was  greater  than  the  critical 
t-statistic,  which  indicated  that  the  differences  were 
significant.  This  can  be  expected  to  occur  one  time  in 
twenty  by  the  construction  of  the  test.  As  long  as  it  does 
not  recur  consistently,  it  was  not  considered  significant f 
since  this  would  occur  with  any  data. 

In  the  case  of  wind  velocity,  however,  the  differences 
were  consistently  greater  than  expected.  In  each  of  the 
three  simulation  runs,  the  mean  of  the  wind  velocities  was 
significantly  lower  than  the  mean  of  the  five-year  pooled 
data.  The  explanation  for  these  differences  was  found  when 
examining  the  distributions  of  the  historic  data.  The 
permissible  code  values  for  reporting  the  wind  range  from 
the  value  zero  to  nine.  The  value  nine  is  used  to  indicate 
unknown  conditions.  This  value  appeared  three  percent  of 
the  rime  in  the  historic  data,   but   was   not   permitted   to 
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appear  in  the  simulated  data.  When  the  nines  were  removed 
from  the  historic  data  the  new  mean  became  1.1571  instead  of 
1.3815,  and  when  this  mean  was  tested  against  the  simulation 
runs  there  was  no  significant  difference  between  them.  It 
appeared  that  the  inclusion  of  the  value  nine  weighted  the 
mean  sufficiently  in  the  historic  data  to  cause  a 
significant  difference  when  compared  with  a  vector  of  wind 
velocity  values  where  the  value  nine  was  not  present.  Even 
if  all  the  historic  observations  had  been  a 
higher-than-normal  code  value,  such  as  six,  the  differences 
between  historic  and  simulated  data  would  not  have  been 
significant.  The  code  value  of  nine  is  permissible  for  the 
other  weather  conditions,  but  the  number  of  times  it 
occurred  was  not  sufficient  to  cause  a  significant 
difference  between  the  historic  data  means  and  the  those  of 
the  simulated  data. 

fin  important  comparison  of  data  was  that  pertaining  to 
nature  of  distress.  The  code  values  of  nature  of  distress 
range  between  00  and  99,  but  not  all  numbers  are  used  to 
describe  a  nature  of  distress.  In  the  case  of  the  Areata 
area  of  interest,  several  of  the  permissible  values  were 
never  used  in  the  data  for  five  years.  It  was  important 
that  the  distributions  of  the  simulated  nature  of  distress 
reflect  the  actual  distributions,  and  this  coincidence  had 
to  be  tested.  The  distribution  of  data  for  nature  of 
distress  is  given  in  the  table  below: 
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Simulated 
Distribution 

4, 

.2 

58, 

.0 

4. 

.8 

3. 

,6 

3, 

,6 

1. 

,0 

7. 

,0 

4. 

5 

0. 

1 

2. 

6 

0. 

4 

DISTRIBUTION  OF  NATURE  OF  DISTRESS 
(Figures  in  Percent) 

Nature  of  Historical 

Distress  Distribution 

Other  Vessel  3.5 
Conditions 

Vessel  Disabled  56.3 
or  Adrift 

Vessel  Disabled  6. 1 
and  Anchored 

Vessel  Aground  3.6 

Vessel  Capsized  4. 4 
or  Sinking 

Vessel  Fire  1.6 
or  Explosion 

Vessel  Flooding  6.1 
or  Sinking 

Vessel  Overdue/Missing  5..  7 

Vessel  Collision  0.4 

Vessel  Disoriented-Lost  3.2 

Vessel  Endangered  0.5 
by  Weather 

Aircraft  Distress  0.8                0.9 

Land  Vehicle  Distress  0,2                 0.0 

Land  Structure  0.3                 0.4 

Diver  in  Distress  0.4                 0.7 

Personnel  Medevac  0*7                 0.9 

Person  Drowning  0.4                 0.6 

Person  in  Water  0.8                 1.3 

personnel  Sickness  2.6                2.5 

Personnel  Injury  1,2                 1.6 

Person  Lost  on  Land  0.3                 0,3 

Other  Miscellaneous  0.3                 0.3 
Conditions 

Flare  Sighting  0.6                 0.7 


The   appropriate   statistic   used  to  determine  if  these 

distributions  were  equal  was  the  chi--square   goodness-of-f it 

2 
test.   The  chi   statistic  is  given  as  follows: 
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.  2   _  2 

Chi   =  \  j[observed  freguency  -  expected  frequency! 
4-—         ex  peel  eel  irequency' 


In  this  equation,  expected  frequency  was  the  historic- 
frequency  and  observed  frequency  was  the  simulated  frequency 
for  eacn  of  the  nature-of-distress  codes.  The  value  of  the 
ratio  was  found  for  each  of  the  entries  in  the 
nature-of-distriss  distribution  table  and  summed  to  find  the 
chi-squared  statistic.  This  statistic  was  computed  as  2.68; 
the  critical  chi-square  statistic  with  22  degrees  of  freedom 
at  the  95  percent  confidence  level  is  12.34.  It  was 
concluded  that  there  was  no  significant  difference  between 
the  two  distributions  of  nature  of  distress.  Sioiiiar 
chi-squared  tests  were  performed  on  all  elements  of  the  data 
to  compare  the  goodness  of  fit  between  historic  data  and 
simulated  data  distributions  with  the  same  results. 

The  conclusion  reached  as  a  result  of  all  these  tests 
comparing  the  historic  data  and  the  data  generated  by  the 
computer  model  of  the  environment  was  that  the  simulation 
data  was  not  significantly  different  from  the  historic  data. 
It  was  concluded  that  the  data  generated  by  the  simulation 
could  be  used  to  augment  the  historic  data  and,  within 
limits,  could  be  used  to  represent  the  environment.  in  the 
area  of  interest. 
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VIe  ANALYSIS  OF  ALTERNATIVES 

After  historic  data  had  been  analysed  and  more  data  had 
been  generated  using  the  computer  model,  the  data  which  was 
available  was  used  to  determine  the  effectiveness  of 
alternate  methods  of  providing  SAR  coverage  for  the  Areata 
area  of  interest.  One  of  the  alternatives  was  the  proposed 
air  station  at  Areata.  A  computer  program  was  written  which 
performed  this  step. 

The  computer  program  was  written  to  compute  times 
required  for  various  Coast  Guard  units  to  arrive  on  the 
scene  of  distress  from  their  "home"  positions.  This  time  is 
termed  response  time.  The  response  times  were  calculated 
for  surface  units  at  Humbolt  Bay  and  Crescent  City,  for 
HH-52A  helicopters  from  San  Francisco,  for  HH-52A  response 
from  the  proposed  air  station  at  Areata,  and  for  HH-3F 
helicopter  response  from  San  Francisco.  The  first  three 
facilities  are  currently  serving  the  area  of  interest,  and 
the  latter  two  were  alternate  options  for  SAR  air  coverage 
in  the  Areata  area.  The  computer  program  calculated  only 
transit  times  for  Coast  Guard  units  from  their  base  to  the 
scene  of  distress  a+  assumed  speeds.  The  response  times  did 
not  include  total  time  between  alerting  the  SAR  facility  to 
the  time  to  either  launch  the  unit  or  get  it  underway.  The 
response  time  did  include  some  "dead  time"  in  calculating 
response  of  aircraft  transiting  between  air  stations  other 
than  Areata  and  the  Areata  airport.  These  delays  were  added 
to  account  for  reasonable  ferry  times,  for  refueling  times, 
and  for  other  logistics  delays. 

The  speed  of  advance  for  surface  vessels  was  set  at 
eighteen  knots  in  ail  cases.  It  was  assumed  that  the  unit 
would  always  travel  at  this  speed  regardless  of  weather.  It 
was   also  assumed  that  anv  surface  unit  could  respond  to  any 
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SAR  incident  anywhere  in  the  area  of  interest.  In 
actuality,  of  course,  the  choice  of  the  appropriate  surface 
unit  would  be  a  function  of  weather,  distance  to  scene, 
distance  offshore,  nature  of  distress,  and  alert  status  of 
the  unit  at  the  time  of  distress. 

The  transit  time  for  an  HH-52A  helicopter  from  San 
Francisco  Air  Station  to  Areata  was  assumed  to  be  2.7  hours 
for  the  210  miles.  Once  the  helicopter  arrived  in  Areata, 
it  was  assumed  that  one  hour  would  be  required  to  refuel  the 
helicopter  before  it  could  be  launched  for  the  SAR  mission. 
The  transit  time  for  HH-3F  helicopters  from  San  Francisco  to 
Areata  was  assumed  to  be  1.75  hours.  If  the  incident 
occurred  less  than  100  miles  from  the  coastline,  the 
helicopter  was  assumed  to  be  capable  of  flying  directly  to 
scene  from  San  Francisco  without  refuelling.  The  speed  over 
ground  for  HH-52A  helicopters  was  assumed  to  be  85  knots 
despite  wind  conditions.  The  speed  used  for  HH-3F 
helicopters  was  120  knots. 

HH-52A  helicopters  from  Areata  and  San  Francisco  were 
assumed  to  be  permitted  to  fly  25  miles  offshore  unescorted. 
If  the  straight- line  distance  from  the  Areata  airport  to  the 
scene  of  distress  was  greater  than  25  miles  but  the  scene 
could  be  reached  by  flying  parallel  to  the  coast  to  the 
correct  latitude  and  then  perpendicular  to  the  coastline  to 
the  scene  of  distress,  this  two-leg  distance  was  computed. 
For  distress  positions  between  25  and  100  niles  offshore,  an 
HH-52A  helicopter  could  conduct  the  mission  only  if  escorted 
by  a  fixed-wing  aircraft.  It  was  assumed  that  1.5  hours 
were  required  to  fly  a  fixed-wing  aircraft  from  San 
Francisco  to  Areata  before  the  helicopter  could  begin  the 
mission. 

For  positions  of  distress  greater  than  100  miles  and 
less  than   300   miles   offshore,   an   HH-3F   helicopter   was 
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required  to  effect  the  rescue.  it  was  assumed  that  six 
hours  were  reguired  to  fly  the  helicopter  to  Areata,  refuel 
it,  man  it  with  a  fresh  crew,  and  launch  it  for  the  mission. 
This  time  was  not  required  if  the  helicopter  were  stationed 
in  San  Francisco. 

The  computer  program  which  calculated  response  times 
for  Coast  Guard  units  used  the  positions  of  distress  which 
were  generated  in  the  computer  simulation,  A  modification 
of  this  program  calculated  response  times  using  data  from 
three  years  historic  incidents.  As  was  expected,  the 
average  response  time  varied  from  year  to  year  because  the 
positions  of  distress  varied.  Besides  using  the  three  years 
historic  data,  seven  years  simulated  data  was  used  to 
determine  response  times,  yielding  the  equivalent  of  ten 
years  activity.  The  assumption  was  made  that  the 
distribution  of  the  positions  of  distress  would  not  change. 
It  was  also  assumed  that  the  distributions  of  the 
percentages  of  the  severity  to  personnel  and  property  codes 
did  not  change.  A  sample  output  from  the  computer  program 
was  included  in  the  Computer  Output  section. 

Without  an  air  station  at  Areata,  California,  the 
calculated  average  response  time  for  all  over-vate- 
incidents  was  0.66  hours,  or  about  40  minutes.  With  an  air 
station  at  Areata,  the  average  response  time  was  0.56  hours, 
or  about  34  minutes.  With  HH-3F  helicopters  at  San 
Francisco,  the  average  response  time  for  the  over-water 
cases  was  the  same  as  that  with  no  air  station  at  Areata,  or 
0.66  hours.  This  was  so  because  in  almost  every  case  a 
surface  unit  arrived  on  scene  in  less  time  than  it  took  for 
the  HH-3F  to  fly  from  San  Francisco  to  Areata.  The  few- 
cases  in  which  this  was  not  true  were  those  cases  occurring 
greater  than  approximately  forty  miles  offshore,  where  the 
response  time  for  air  and  surface  facilities  was  equal. 
Historically,  SAR  incidents  have  occurred  more  than  5)0  miles 
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offshore   in   the   area  of  interest  less  than  ten  percent  of 
the  time. 

The  effectiveness  of  the  proposed  air  station  can  not 
be  measured  as  a  reduction  in  average  response  time  to  all 
incidents,  since  a  helicopter  would  be  used  only  in  moderate 
or  severe  cases.  Moderate  and  severe  cases  historically 
occurred  about  thirty  percent  of  the  time.  The  average 
response  time  to  moderate-severity  cases  without  Areata  Air 
Station  was  0.7  hour  (about  42  minutes).  The  range  for  that 
average  was  from  0.48  hour  to  0.91  hour  in  different  years. 
The  average  response  time  with  Areata  Air  Station 
established  was  0.42  hour  (about  25  minutes).  The  range  of 
yearly  means  was  from  0.30  hour  to  0,55  hour.  It  was 
assumed  that  a  helicopter  would  not  launch  from  San 
Francisco  on  a  moderate-severity  case. 

The  cases  with  the  greatest  severity  to  persons  or 
property  reported  averaged  about  ten  percent  in  the  historic 
data.  The  average  response  time  to  these  cases  without  an 
air  station  at  Areata  (but  with  helicopter  assistance  from 
San  Francisco)  was  0.75  hour  (about  45  minutes).  The  yearly 
averages  ranged  between  0.58  hour  and  one  hour*  The  average 
response  time  with  an  air  station  at  Areata  was  0.42  hour 
(about  25  minutes)..  This  yearly  average  ranged  between  0.25 
hour  and  0.74  hour.  The  average  response  time  with  HH-3F 
helicopters  at  San  Francisco  was  just  slightly  less  than  the 
average  with  BH-52A  helicopters  at  San  Francisco. 

The  average  response  time  to  both  moderate  or  severe 
cases  without  Areata  Air  Station  and  with  either  HH-52AS  or 
HH-3FS  at  San  Francisco  was  0.72  hour  (about  43  minutes) . 
With  an  air  station  at  Areata,  the  average  response  time  was 
0.41  hour  (about  25  minutes). 

In  addition  to  making  these  absolute  determinations  of 
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effectiveness,  it  was  necessary  to  compare  the  effect  of  the 
additional  air  station  relative  to  the  "acceptable  SAR 
response  time"  as  set  forth  in  the  Aviation  Plan.  The 
tables  below  show  the  percentage  of  times  the  response  times 
of  the  Coast  Guard  units  were  greater  than  the  Aviation  Plan 
standard.  The  tables  include  both  historic  data 
computations  and  simulated  data  computations. 


PERCENTAGE  OF  RESPONSE  TIMES 
GREATER  THAN   ONE- HALF  HOUR 

1971     1972     1973     Run1     Run2     Run3 


Without 
Areata 
Air  Sta. 

0.50 

0.63 

0.57 

0.58 

0.61 

0.60 

With 
Areata 
Air  Sta. 

0.  17 

0.56 

0.45 

0.36 

0.40 

0.41 

Acceptable 

0.25 

0.2  5 

0*25 

0.2  5 

0.25 

0.25 

PERCENTAGE  OF  RESPONSE  TIMES 

GREATER  THAN   ONE  HOUR 

1221  1212.  1973    RunJ    Run2    Bun3 

Areata1        0.12     0.12     0.30     0.16     0,18    0.20 
Air  Sta. 

With 

Areata         0.08     0.10     0,22     0.08    0.11     0.12 

Air  Sta. 

Acceptable     0.10     0.10     0.10     0.10     0.10     0.10 

The  "acceptable  response  time"  is  not  mere  than  2  5%  of 
cases  requiring  one-half  hour  and  not  more  than  10% 
requiring  one  hour.  The  addition  of  a  Coast  Guard  air 
station  at  Areata  does  not  permit  the  Coast  Guard  to  meet 
the  established  standard  in  all  incidents,  but  it  does 
permit  them  tc  reach  the  standard  a  greater  number  of  times. 

A  much  more  favorable  case  was  made  for  an  air  station 
at  Areata  if  the  incidents  occurring  over  land  were 
considered.   In  these  incidents,   the   transit   time   for   a 
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helicopter  from  San  Francisco  to  Areata  was  saved:  1.75 
hours  for  an  HH-3F  and  2.7  hours  for  an  HH-52A.  Undoubtedly, 
the  number  of  Coast  Guard  responses  to  land  incidents  would 
increase  if  an  air  station  was  established  at  Areata,  but 
there  was  no  means  of  analyzing  the  effectiveness  of  the  air 
station  relative  to  the  SAfi  facilities  presently  serving  the 
area.  Whether  this  supplement  to  present  over-land 
resources  justifies  the  establishment  of  a  new  air  station 
or  not  is  a  policy  decision. 
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VII®   CONCLUSIONS 

Since  there  was  no  established  relationship  between 
lives  lost  or  saved  and  response  time  when  this  thesis  was 
written,  the  results  of  this  study  must  be  stated  merely  as 
information  for  the  decision  maker.  An  air  station  at 
Areata,  California  will  probably  reduce  Coast  Guard  response 
time  to  SAR  incidents  involving  moderate  and  severe  danger 
to  persons  and  property  in  the  area  of  interest  an  average 
of  fifteen  to  twenty-five  minutes.  It  will  reduce  the 
percentage  of  times  which  Coast  Guard  units,  other  than 
fixed-wing  aircraft,  cannot  arrive  on  scene  of  a  distress 
which  is  up  to  150  miles  offshore  within  one-half  hour  from 
approximately  sixty  percent  to  approximately  forty  percent 
of  the  time.  Its  presence  will  reduce  the  percentage  of 
times  that  response  time  in  this  zone  is  greater  than  one 
hour  from  approximately  twenty  percent  to  approximately  ten 
percent.  These  numbers  are  based  on  case  loads  of 
approximately  250  to  350  per  year  in  the  Areata  area.  Of 
these  cases,  approximately  75  to  105  cases  are  of  moderate 
or  severe  severity.  These  figures  do  not  account  for  the 
possible  delay  in  response  by  aircraft  caused  by  low- 
ceilings  or  low  visibilities  at  the  Areata  Airport. 
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APPENDIX  B 
Histograms  of  Five-Year  Data 
in 
Areata  Area  of  Interest 

The  following  seven  histograms  are  the  distributions  of 
the  code  values  reported  in  the  historic  data  for  the  Areata 
area  of  interest.  Each  graph  is,  in  effect,  five  separate 
histograms  superimposed  upon  each  other.  Rather  than  showing 
the  distributions  of  the  code  values  for  any  specific  year, 
the  graphs  are  intended  to  show  the  general  distributions 
experienced  in  the  area  of  interest  as  reported  in  the 
Assistance  Reports.  These  distributions  were  pooled,  or 
averaged,  to  give  the  cumulative  distributions  used  in  tne 
data  analysis  aDd  the  computer  simulation  model.  A  table  of 
code  values  depicted  in  this  appendix  can  be  found  in 
Appendix  G.  These  histograms  were  drawn  using  the  CALCOMP 
Model  765  and  the  plotting  package  [Ref.  15]  at  the  Naval 
Postgraduate  School. 
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APPENDIX  C 
Histograms  of  Pooled  Data 

and 
Simulation-Generated  Data 

The  following  six  histograms  were  included  to  shou  the 
pooled  distribution  of  code  values  for  the  five  years 
historic  data  in  the  Coast  Guard  records  and  the 
distribution  of  the  pooled  data  from  three  simulation  runs. 
The  hypothesis  that  these  distributions  were  equal  was 
accepted  in  an  analysis-of-means  test  in  the  chapter  of  the 
thesis  dealing  with  validation  of  the  computer  model.  A 
table  of  code  values  depicted  on  these  histograms  can  be 
found  in  Appendix  G. 


73 


MONTH 

CODE  FREQUENCV  DISTRIBUTIONS 

HJSTGRICRL  DRTR  i9?0  -  19?4 
"IMUlHTIQN  GENERATED  DfiTfi 


c 


♦  33 


h~ 

tf 

z: 

II 

C_D 

^1 

cx 

W 

UJ 

Q_ 

fc 

iQ 


1        2 


_> 


4        5        6        ?        3 
MONTH 


10      11     12 


"\ 


7U 


Sf.UERIT^   TO   PERSONNEL 

CODE_FREQUENCV  DISTRIBUTIONS 
HISlQRICRL  DRTR  L9?0  -  1974 
SIMUlRIION  GENERATED  DRTR 


5a 


SO 


..  ia 


f— 


o 


2Z 


10 


0 


1 


SEUERITV  CODE 


IS 


SEUERITV   TQ   PROPERTY 

CODE  FREQUENCY  DISTRIBUTIONS 
HISTORICRL   DRTR   .1.970   -   1374 
SIMULRTION   GENERRTFR   DRTR 


so 


<  0 


M 


h  -• 

ZI 

11 

«J_J 

_ 

LJ 

CC 

u 

LU 

a 

CI. 

20 


0 


1 


2 


SEUERITV   CODE 


76 


50 


WIND  UELOCITV 

CODE  FPPfQUEHCV  DISTRIBUTIONS 
HISTQRICRL  DRTR  1970   -   1974 
SIMUlRTIQN  Gf.MLRRTffi  DRTR 


,  io 


UJ     - 


JQ 


4 


Q 


1 


n 
*_ 


3  4  5 

HIND  STRTf"  CODE 


7 


5 


77 


SER  STATE 
code  frequency  distributions 
historicrl  grip.  \9?q  -   1974 

MUlRTION  GENERRTFD  DRTfi 


h-    2 

i. 

I 

+ 

t 

CC    u 
CL.     i/i 

i 

i- 

CT> 

w                           * 

■ 

CE 

_J 

i 

_ J 

1 

1 

— -=d 

0 


L 


3  4    5     6 

SIR  STRTE".  CODE 


? 


a 


7a 


UISIBILITV 

CQDE  FREQUENCE  DISTRIBUTIONS 
HISTORICAL  DfiTfi  .L97Q  -  1974 
"IMULfiTION  GENERATED  DfiTfi 


C 


t   « 


10 


^ 


Q 


<■*> 


3  4  5  6 

UISIBILITV  CODE 


? 


a 


79 


APPENDIX  D 

Results  of  Analysis-of-Means  Tests 

(Estimation  of  Linear  Contrasts) 

The  tables  in  this  appendix  are  a  graphic 
representation  of  the  result  of  mathematical  tests  performed 
on  the  historic  data  and  the  simulation  data.  The 
estimation-of-linear-contrasts  test  was  performed  on  each 
data  element  to  determine  if  the  means  of  the  code  values 
for  each  year  were  equal.  The  test  was  performed  in  the 
analysis-of-data  chapter  to  determine  if  the  means  of  the 
code  values  varied  from  year  to  year.  It  was  performed  in 
the  validation-of-computer-raodel  chapter  to  determine  if  the 
means  of  the  code  values  generated  by  the  computer  model 
were  different  from  the  historic  data.  These  tests  are 
described  in  Bolch  and  Huang  [Ref.  2]. 
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APPENDIX  E 

Correlation  Matrices 
(Measures  of  Relationships  Between  Two  Variables) 

The  tables  in  this  appendix  show  the  linear 
relationships  between  two  data  elements  for  the  historic 
data  from  the  years  1970  through  1974  for  the  Areata  area  of 
interest.  If  one  data  element,  such  as  the  code  value  for 
sea  state,  always  increased  whenever  another  element,  such 
as  wind  velocity,  increased,  the  correlation  would  approach 
the  value  1.0.  If  the  value  of  the  correlation  coefficient: 
is  negative,  it  indicates  a  negative  relationship. 
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APPENDIX  F 
Charts  of  Positions  of  SAB  Incidents 
(1971  through  1973) 

The  following  computer-generated  charts  are  scaled 
approximations  of  positions  of  SAR  incidents  in  the  Areata 
area  of  interest  for  the  years  1971  through  1973.  The 
charts  are  included  only  to  show  concentrations  of  positions 
and  the  relative  positions  of  salient  locations  within  the 
area  of  interest.  Although  not  apparent  on  the  charts,  more 
than  one-third  of  the  distresses  for  each  of  the  years 
occurred  near  the  mouth  of  the  Klamath  River,  just  south  of 
Crescent  City. 

These  charts  were  drawn  to  scale  using  the  CALCOMP 
Model  765.  The  software  for  drawing  the  charts  was  made 
available  through  the  Plotting  Package  available  at  the 
Naval  Postgraduate  School  [Ref.  15  J.  In  order  to  draw  the 
grid  and  the  coastline,  all  points  were  stored  in  an  array 
in  the  program  which  controlled  the  pen  to  draw  the 
outlines.  All  data  points  and  all  labelling  were  done  in 
the  same  manner.  The  program  is  highly  adaptive;  the  author 
drew  similiar  charts  for  the  activities  of  the  Coast  Guard 
in  the  San  Piego  and  the  San  Francisco  area. 
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APPENDIX  G 
Selected  Code  Values  for  Assistance  Reports 

This  appendix  is  a  table  of  code  values  for  selected 
elements  of  the  SAR  environment  which  are  reported  on  the 
Coast  Guard  Assistance  Report  (CG  3272) .  The  elements  which 
were  selected  were  those  which  were  analyzed  in  the  data 
analysis  chapter  and  those  which  were  simulated.  This  table 
was  included  to  illustrate  to  the  reader  who  is  not  familiar 
with  the  Coast  Guard  SAR  Assistance  Report  system  the  method 
in  which  the  conditions  are  reported.  It  should  be  noted 
that  the  code  values  do  not  represent  egual  increaents  on 
the  continuum  of  possible  conditions.  The  spread  of 
distance  offshore,  for  example,  is  three  miles  for  code 
value  one  and  150  miles  for  code  value  seven. 

NATURE  OF  DISTRESS 


Vessel  Conditions 


00  -  Other  vessel  condition: 

01  -  Disabled,  adrift 

02  -  Disaoled,  anchored 

03  -  Aground 

04  -  Capsized 

05  -  Eire/Explosion 

06  -  Flooding/Sinking 

07  -  Overdue/Missing 

08  -  Collision 

09  -  Unfamiliar  with  area/ 

Disoriented 

11  -  Endangered  by  weather 

12  -  Endangered  by  ice 

Aircraft  Conditions 

20  -  Other  aircraft  conds. 

21  -  Ditch/forced  landing 

22  -  Crash 

23  -  Low  on  fuel 

24  -  Bail  out 

25  -  Fire/Explosion 

27  -  Overdue/Missing 

28  -  Mechanical  casualty 

29  -  Unfamiliar  with  area/ 

Disoriented 


Lcin_cl   or    Offshore   Structure 
Conditions 

50  -  Other, 

51  -  Flooding 

55  »■  Fire,  Explosion 
57  -  Hissing 

Diver  Conditions 

— — — — — ~  -r 

70  -  Other  diver  conditions 

71  -  Stranded 

72  -  Personnel  in  water 

73  -  Bends 

74  -  Air  Embolism 

75  -  Emphysema 

76  -  Equipment  failure 

77  -  Overdue/Missing 

78  -  Predator  attack 

Personnel^  Other  than  5iy.§£ 


00  -  Other  personnel  conds. 

81  -  Drowning 

82  -  Personnel  in  water 

84  -  Sickness 

85  -  Injury 
87  -  Hissing 
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Land  Vehicle  Conditions 


Miscellaneous  Conditions 


40  -  Ot 

42  -  Cr 

43  -  Ov 


her  land  vehicle  conds  90  ~  Other  conditions 
ash  (accident)  98  -  Flare  sighting 

erdue/Missing  99  -  Unknown 


0 
1 
2 

3 
4 


0 
1 

2 
3 
9 


0 
1 

2 
3 
9 


0 
1 
2 
3 
4 


0 
1 

2 

3 
4 


M 

Lan 
0  - 
3.  1 
10, 
25. 


STANCE  OFFSHORE 

d 

3  miles 

-  10  miles 
1-25  miles 
1  -  50  miles 


5-50.1-  100  miles 

6  -  100.  1  -  150  miles 

7  -  150.  1  -  300  miles 

8  -  Greater  than  300  miles 

9  -  Unknown 


SEVERITY  OF  DISTRESS  -  PERSONNEL 


None  -  No  personnel  were  involved. 

Small  -  No  immediate  or  foreseeable  danger  to 

personnel. 

Moderate  -  Some  danger  than  personnel  might  be  lost. 

Severe  -  Personnel  were  in  danger  of  loss  or  were  lost, 

Unknown. 

SEVERITY  OF  DISTRESS  -  PROPERTY 

None  -  No  property  was  involved. 

Small  -  No  immediate  danger  to  property. 

Moderate  -  Some  danger  that  property  might  be  lost. 

Severe  ■>•  Property  was  in  danger  of  loss  or  was  lost. 

Unknown. 


SEA  STATE 


Calm 

1  -  2  feet 
3-4  feet 
5-6  feet 


4-7-10  feet 

5-11-20  feet 

6  -  Greater  than  20  feet, 

8  -  Land 

9  -  Unknown 


WIND 


Calm 
0.  1  - 
10.  1 
20.  1 
30.  1 


10  knots 

-  20  knots 

-  30  knots 
-40  knots 


0  -  1/4 
Greater 
Greater 
1.1-3 
3.  1  -  5 


mile 

than  1/4 
than  1/2 
miles 
miles 


1/2  mile 
1  mile 


5 
6 
7 

8 

9 

- 

40.  1  -  50  knots 
50.1  -  60  knots 
6  0.1  -  7  0  knots 
Over  7  0  knots 
Unknown 

ILITY 

e 

5  - 

6  -• 

7  - 

9  - 

5.1  -  10  miles 
10.1-  15  miles 
15.  1  -  20  miles 
Unknown 
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SAMPLE  COMPUTER  OUTPUT 

The  following  is  a  sample  output  from  the  simulation 
program.  It  represents  the  cases  generated  for  the  month  of 
March.  In  addition  to  the  printed  output  shown,  the  program 
also  punched  a  card  with  the  elements  of  data  punched  in  the 
same  format  as  was  found  in  the  historic  data  which  was 
received  on  cards.  The  cards  which  were  punched  were 
retained  and  used  in  the  data  analysis  and  validation,  of  the 
simulation  model. 
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A  COMPUTER  PROGRAM  TO  SIMULATE 
THE  SEARCH  AND  RESCUE  ENVIRONMENT 

OF  THE 
PROPOSED  COAST  GUARD  AIR  STATION 

AT 
ARCATA,  CALIFORNIA 
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DIMENSI CN 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 

ILAND 


LLIND(283) ,XCUM(288) 
WXNM729),  I  STATE  (729) 
HYST(23) ,HCUM(23, 16) 
XT  I  ME ( 24) 
AVG( 12,28) 


IS  AN 

HAPPENi 

IFLAG  IS  AN 

LATLON 


INDICATOR  TO  SHOW  IF  A  CASE 
ID  ON  LAND. 
INDICATOR  USED  IN  SUBROUTINE 
TO  SHOW  THAT  THE  POSITION  OF  A 


SAR  INCIDENT  WAS  FARTHER  WEST  THAN 

12  5-OOW. 
IAIR   IS  AN  INDICATOR  TO  SHOW  IF  A  HELO 

WOULD  BE  LAUNCHED  ON  THIS  CASE  DUE 

TO  THE  SEVERITY. 
ITOTCS  IS  A  COUNTER  FOR  THE  TOTAL  CASESo 
IHELCS  IS  A  COUNTER  FOR  HELICOPTER  CASES 

ILAND  =  0 
IFLAG  =  0 
ITOTCS  =  0 
IHELCS  =  0 
IAIR  =  0 


THIS  SECTION 
FOR  THE  DATA 
SUBROUTINES. 


ESTABLISHES  COMMON 
WHICH  WILL  BE  USED 


STORAGE 

IN  SEVERAL 


C 
C 


CCMMON 
COMMON 
CCMMON 
CCMMON 
CCMMON 
CCMMON 
CCMMON 
CCMMON 
CCMMON 
CCMMON 
CCMMON 
CCMMON 
CCMMON 
CCMMON 
CCMMON 
CCMMON 
CCMMON 
CCMMON 
COMMON 
COMMON 
CCMMON 
CCMMON 
CCMMON 
CCMMON 
CCMMON 
CCMMON 


/COM1 
/COM2 
/COM3 

/C0M4 
/C0M5 
/C0M6 
/COM7 
/COMS 
/COM9 
/COM1 
/COM! 
/COM1 
/COMl 
/COMl 
/COM  I 
/COMl 
/COMl 
/COMl 
/COMl 
/COM  2 
/COM2 
/COM2 
/COM2 
/COM? 

/com 

/COM2 


/XCU 
/ISE 
/  1ST 

/WXN 
/  riY 
/HCU 
/LAT 
/MIN 
/LMI 
O/XT 
1/NC 
2/AV 
3/DI 
4/IL 
5/1  H 
6/  IA 
7/IS 
8/IW 
9/1  V 
O/IP 
1/ME 
2/ NT 
3/NU 
4/1P 
5/LL 
6/  D 


M 

ED 

ATE 

R 

ST 

M 

S 

NS 

I  ME 

ASES 

G 

ST 

AND 

ELCS 

I  R 

EA 

I  N  D 

IZ 

ER 

W 

YM 

DES 

ROP 

IND 

DDIST 
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CALL    OVFLOW 
C 
C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c  c 

C  THE  HISTORICAL  CUMULATIVE  DISTRIBUTIONS  C 

C  FOR  THE  LATITUDE  AND  LONGITUDE  OF  THE  SAR  C 

C  INCIDENTS  FOR  THE  AREA  OF  CONSIDERATION  C 

C  FOR  THE  ARCATA  COAST  GUARD  AIR  STATION  C 

C  (LAT  40-00N  TO  LAT  41-00N,L0NG  124-00W  TO  C 

C  LONG  125-OOW)  ARE  READ  OFF  CARDS  AND  STORED  C 

C  IN  AN  ARRAY  CALLED  XCUM  WITH  A  SIZE  OF  C 

C  288  (12  X  24)  LOCATIONS.   THIS  MEANS  THAT  C 

C  EACH  OF  THE  288  SQUARES  REPRESENTS  A  C 

C  AREA  FIVE  DEGREES  IN  LATITUDE  BY  FIVE  C 

C  DEGREES  IN  LONGITUDE.   THE  CUMULATIVE  C 

C  NUMBERS  IN  THE  STORAGE  LOCATIONS  REPRESENT  C 

C  THE  NUMBER  OF  CASES  WHICH  OCCURRED  IN  THE  C 

C  YEARS  1971  THROUGH  1973  IN  THAT  AREA.  C 
C 

READ  (5,100)  (XCUM(  I)  ,1  =  1,288) 

100  FORMAT  ( 12 (F5.1)  ) 

DO  10  N  =  1,288 

LLIND(N)  =  N 

10  CONTINUE 
C 

c  c 

C  AN  ARRAY  CALLED  ISTATE  IS  FILLED  WITH  C 

C  NUMBERS  WHICH  IDENTIFY  EACH  ARRAY  LOCATION.  C 

C  THIS  ARRAY  WILL  BE  USED  IN  SUBROUTINE  C 

C  NATURE  TO  CALCULATE  WEATHER  CONDITIONS.  C 

c  c 

C  THE  HISTORICAL  DISTRIBUTION  OF  WEATHER  CON-  C 

C  DITIONS(WIND  VELOCITY,  SEA  STATE,  AND  VI  SI-  C 

C  3ILITYHS  READ  OFF  CARDS  AND  STORED  IN  AN  C 

C  ARRAY  CALLED  WXNR  OF  SIZE729  (9X9X9).  C 

C  THE  DISTRIBUTION  IS  BASED  ON  THE  OBSERVED  C 

C  WEATHER  AT  THE  SCENE  OF  SAR  INCIDENTS  IN  C 

C  THE  AREA  OF  STUDY  FOR  THE  YEARS  1970  TO  C 

C  1974.  C 
C 

DC  20  I  =  1,729 

ISTATE( I)  =  I 

20  CONTINUE 

RE£D(5,200)  (WXNR(I)fl  =  1,729) 

200  FORMAT  (9(F7.3) ) 
C 

C  THE  "NATURE  OF  DISTRESS"  HISTORICAL  CUMU-  C 

C  LATIVE  DISTRIBUTION  IS  READ  OFF  CARDS  ALONG  C 

C  WITH  THE  CUMULATIVE  DISTRIBUTIONS  OF  DEGREE  C 

C  OF  SEVERITY.   THE  DATA  IS  BASED  UPON  FIVE  C 

C  FIVE  YEARS  SAR  ACTIVITY.   THE  NATURE  OF  C 

C  DISTRESS  ARRAY — HYST--IS  OF  SIZE  23.  C 

C  BECAUSE  SOME  OF  THE  POSSIBLE  TYPES  GF  C 

C  DISTRESS  DID  NOT  OCCUR  IN  THE  EUREKA  AREA  C 

C  IN  THE  FIVE  YEARS  HISTORICAL  DATA,  THEY  C 

C  ARE  PRESUMED  NEVER  TO  OCCUR  IN  THE  RUNNING  C 

C  OF  THE  MODEL.  EACH  TYPE  OF  DISTRESS  HAS  A  C 

C  4X4  MATRIX  REPRESENTING  DEGREE  OF  SEVER-  C 

C  ITY  FOR  THAT  TYPE  OF  DISTRESS.   THE  SIZE  C 

C  OF  THE  STORAGE  REPRESENTING  DEGREE  OF  SEVER-  C 

C  IS  23  X  16.  C 

c  o 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

READ  (5,300)  ( HY ST ( K )  , K  =  l , 2 3 ) 
300  FORMAT  (8F10.2) 

DO  30  K  -  lf23 

READ( 5,40  0)  (HCUM(K,J)  ,  J  =  l,  16) 
400  FCRMAT( 16(F5.2) ) 
30  CONTINUE 
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c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


c 


c 


C  THE  FIRST  STEP  OF  THE  SAR  SIMULATION  IS  TO  C 

C  GENERATE  THE  TIME  OF  THE  DISTRESS.   BY  DATA  C 

C  ANALYSIS  IT  WAS  DETERMINED  THAT  THE  PRO-  C 

C  BABiLITY  OF  A  SAK  INCIDENT  OCCURRANCE  C 

c  varied  hour  by  hour,  day  by  day,  and  month  c 

C  BY  MONTH.   INI  ORDER  TO  SIMULATE  THE  SAR  C 

C  OCCUKANCES  WHILE  KEEPING  THE  MODEL  WITHIN  C 

C  REASONABLE  BOUNDS  THE  FISCAL  YEAR  DATA  WAS  C 

C  ADAPTED  TO  A  "MODEL  YEAR".  THE  MODEL  YEAR  C 

C  CONSISTS  OF  336  DAYS.  EACH  OF  TWELVE  MONTHS  C 

C  HAS  FOUR  WEEKS  EACH  OF  WHICH  HAS  SEVEN  DAYS.  C 

C  THE  FIRST  OF  EACH  MONTH  IS  A  SUNDAY,  AND    .  C 

C  EACH  SEVENTH  DAY  IS  A  SUNDAY.   EACH  MONTH  C 

C  HAS  28  DAYS.  C 

C  C 

CONTINUE 

C  IN  ORDER  TO  DETERMINE  THE  DISTRIBUTIONS  CF  C 

C  SAR  INCIDENTS  FOR  THE  MODEL  YEAR,  THE  SAR  C 

C  DATA  FOR  FOUR  YEARS  WAS  ANALYZED.   THE  WEEK  C 

C  OF  EACH  CALENDAR  MONTH  WAS  REARRANGED  SO  C 

C  THAT  IT  CORRESPONDED  TO  THE  MODEL  YEAR  AND  C 

C  THE  SAR  OCCURANCES  HAPPENING,  FOR  INSTANCE  C 

C  ON  THE  SECOND  MONDAY  OF  JUNE  FOR  EACH  OF  C 

C  THE  FOUR  YEARS  WAS  AVERAGED  IN  ORDER  TO  C 

C  OBTAIN  THE  AVERAGE  (EXPECTED)  NUMBER  OF  C 

C  SAR  OCCURANCES  FOR  THAT  DAY.   THESE  FIGURES  C 

C  ARE  READ  INTO  AN  ARRAY  CALLED  AVG(M,N).  C 

C  ,(Mn  IS  THE  MONTH  NUMBER  AND  "N"  IS  THE  DAY  C 

C  NUMBER.  C 

c  c 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C 

DC  40  M  =  1,12 

READ  (5,500)  (AVG(M,N),N  =  1,28) 
500  FORMAT  U6F5  . 1 ,  /  12P5.  1 ) 
40  CONTINUE 
C 
C 

C  THE  HISTORICAL  CUMULATIVE  DISTRIBUTION  OF  C 

C  TIMES  OF  DAY  FOR  SAR  INCIDENTS  IS  READ  INTO  C 

C  AN  ARRAY  CALLED  ITIME. 


C 


READ    (5,600) (XTIMEl I) ,1=1,24) 
600    FORMAT     U6F5.1  ,/ 8F5.1  ) 


C 

C+t*    *&*    J*  >*••  +}"   ■*'■'  **-  *f*   ^'■"  >}*■   ***  "**^  ••'-'  » ■■*■  *J*  ***■  JU  v>  O^  *■*.  %!*-  i.ir   **»  »V  *?+  *■'»  *V  *f*   ****  *J+   *■***  '*'■'  ***■  -»*•*  n'*  >"»•  *$*   *•&?  *t+   -*+   *'i*  *■'-'  *t+   *A*  *Af    >J*  J*  ***  *A*   »' -  -J*-  O*  -A*   •_!■»  **■•  **.*  >^  xl*/*"* 
rr  *«-  -v  rj-  t-  *v  *v  ^  ^r*  *r  -.-  t*  -;•»  *»*■  v  *r-  t-  *r»  t-  -v  a*  -y*  -v  "r  -c-  ***■  *r  *r-  *r  *y  *e  *r  *r-  *r  *i*  -v-  -v*  -r  *c  -r  nr  *r  -r  T*  ■¥>  ***•  *v  *e  t  nr  t*  *c  -)*•  V  *v  V  *r  -v*  l^, 

C 

C  EXECUTION  OF  SIMULATION  BEGINS 

C 

*i"  -?  -v-  *i*  -i'  V-  «y*  5?  *i*  ^  ^  *v-  <?  *?  ^  ■*?  'r  '!*  ^r  -v  '.*  -*•  -£  *v  *•*  *?  ^  ^  -•"  •v:  'i*  3<^  •  r  n5  A  -  •**  -^*  «v  *r  '<c  *r  "r  V  •r-  -S*  V  *»*  *?  -*^  *«*  't»  *»■*  ^*  -?  -v  -r>  '«»  *><•  (_, 

C 

c 

C  THE  SIMULATION  IS  INSTRUCTED  TO   RUN  FOR  C 

C  TWELVE  MONTHS  OF  TWENTY-EIGHT  DAYS  EACH.  C 

C 

DC  70  M  =  1,12 

DO  70  N  =  1,28 
C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c  c 

C  AS  THE  MODEL  PROGRESSES  THROUGH  THE  YEAR,  C 

C  THE  CALENDAR  IS  ADVANCED.   A  SUBROUTINE  C 

C  IS  CALLED  TO  PRINT  OUT  THE  DATE,  THE  DAY  OF  C 

C  WEEK,  AND  THE  MONTH.  C 

C  c 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

CALL  DATE  (M,N) 
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c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c  c 

C  SUBROUTINE  ICASES  IS  CALLED  TO  DETERMINE  C 

C  THE  NUMBER  OF  CASES  ( NCASES)  WHICH  WILL  C 

C  OCCUR  DURING  THE  CURRENT  MODEL  DAY.  C 

C  IN  ORDER  TO  DETERMINE  THE  NUMBER  OF  CASES  C 

C  WHICH  WILL  OCCUR,  THE  SUBROUTINE  IS  PASSED  C 

C  THE  AVERAGE  NUMBER  OF  CASES  WHICH  OCCURRED  C 

C  ON  THIS  MODEL  DAY  FROM  THE  HISTORICAL  C 

C  CUMULATION  ARRAY  AVG(M,N).  C 

C  c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

CALL  ICASES  (M,N) 

IF  (NCASES.LT.l )  GO  TO  70 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  c 

C  THE  NEXT  TASK  IS  TO  DETERMINE  WHAT  TIMES  C 

C  OF  THE  DAY  THAT  THE  CASES  WILL  OCCUR.   FROM  C 

C  HISTORICAL  DATA  IT  CAN  BE  DETERMINED  WHAT  C 

C  THE  DISTRIBUTIONS  OF  SAR  OCCURANCES  ARE.  C 

C  GIVEN  THE   NUMtfER  OF  CASES  THAT  WILL  OCCUR  C 

C  ON  ANY  DAY  (NCASES)  THE  NEXT  STEP  IS  TO  C 

C  SIMULATE  AT  WHAT  TIMES  THEY  WILL  OCCUR.  C 

C  SUBROUTINE  NT  I  ME  IS  CALLED  TO  DETERMINE  C 

C  WHAT  TIMES  OF  DAY  THE  SAR  INCIDENTS  WILL  C 

C  OCCUR.  C 

C  C 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

DO  71  NCSE  =  1, NCASES 

CALL  NT  I  ME  (NCSE) 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  c 

C  SUBROUTINE  NATURE  IS  CALLED  IN  ORDER  TO  C 

C  DETERMINE  THE  NATURE  OF  DISTRESS.  C 

C  SUBROUTINE  LATLON  IS  CALLED  TO  GENERATE  C 

C  THE  LATITUDE  AND  LONGITUDE  OF  THE  SAR  C 

C  INCIDENT,   I  FLAG  IS  A  MARKER  WHICH  C 

C  IS  USED  TO  INDICATE  THAT  NATURE  HAS  C 

C  GENERATED  A  CASE  WHICH  OCCURS  INLAND.  C 

C  SUBROUTINE  DISTNC  COMPUTES  AND  PRINTS  THE  C 

C  DISTANCE  FROM  COAST  GUARD  STATIONS  C 

C  TO  THE  SCENE  OF  DISTRESS.  C 

C  SUBROUTINE  WEATHR  GENERATES  THE  WEATHER  AT  THE     C 

C  SCENE  OF  DISTRESS.  C 

C  c 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

CALL  NATURE 

CALL  LATLON  ( IFLAG; 

IF  (IFLAG.  EQ.l)  GO  TO  50 

CALL  DISTNC 
50  CALL  WEATHR 

IFLAG  =  0 

ILAND  =  0 

I  AIR  =  0 

ITOTCS  =  ITOTCS  *  1 
C 
C         DATA  CARD  IS  PUNCHED. 


C 


WRITE  ( 7,700)  MtNiNTYMt NODES, MEW, L AT  *  MI NS , LMINS , I  PER, 
1  I.  PROP,  I  SEA,  I  WIND,  I  VIZ  yDDDI  ST  „.,,»,,,, 

7  00  FORMAT! 1 6X, I  2, 3X, 212, 6X, 12,11, 212, 3X, I Z , IX, 211, lX,3il, 
14X,F4.0) 
71  CONTINUE 
70  CONTINUE 
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c 
c 
c 
c 
c 


THE  SUMMARY 
IS  PRINTED. 


OF  THE  OPERATIONS  OF  THE  PERIOD 


WRITE  (6,800)  ITOTCS,  IHELCS 
800  FORMAT  (///«  THE  TOTAL  CASES 


1'PERIOD 
STOP 
END 


GENERATED  FOR 


IS 


13 


TOTAL  CASES  WITH  HELO 


THIS 
IS  « 


c 

c 


14 


C  *  *  £  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  £  *  *  *  *  *  *  *  *  *  sic  #  *  *  #  *  *  #  *  *  if  if  if  if.  *  *  if  if  #  *  #  4  *  #  *  *  *  *  *  r 

C  r  C 

SUBROUTINE    DATE     (M,N) 
C  C 

C* *********************************** *********  3Ct*##£##;;c*:>;<*3Cc*Q 

c 
c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c  c 

C  THIS  SUBROUTINE  RECEIVES  THE  VALUES  "M"  C 

C  AND  ««N"  FROM  THE  MAIN  PROGRAM.   "M"  IS  THE  C 

C  MONTH  NUMBER  (1  TO  12)  AND  "N"  IS  THE  DAY  C 

C  NUMBER  (L    TO  23).  IT  TRANSFORMS  THESE  NUMBERS  C 

C  AND  PRINTS  THE  CORRECT  DATE, DAY,  AND  MONTH.  C 

C  C 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

REALMS  DAY(7)/' SUNDAY '  ,  ' MONDAY '  , ■ TUESDAY • ,' WEDESDAY'  , 
1 'THURSDAY* , 'FRIDAY' , • SATURDAY'/ 
C 

REAL* 8  MONTH ( 1 2 ) /»  JANUARY •  ,  ' FE BRUARY «  ,' MARCH •  ,' APR  I L« , 
1 'MAY' ,  «  JUNE'  , 'JULY' , « AUGUST' ,' SE PT. ',' OCTOBER • , 
2sNOVEMBER« , '  DECEMBER' / 
C 

WRITE  (6,100)  N, MONTH (7) 
100  FORMAT  (//«  TODAY  IS  THE  ',12,'  DAY  OF  • , A8 ) 

IF(N.EQ.l.G:H.N.EQ.3.0R„N.EQ.lL>.0R.N.EQ.22)  GO  TO  11 
IF(N. EQ. 2. OK. N.EQ. 9. OR. N.EQ. 16. OR. N.EQ. 23)  GO  TO  12 
IF(N.EQ. 3. OR. N.EQ. 10. OR. N.EQ. 17. OR. N.EQ. 24)  GO  TO  13 
IFUM.EQ, 4. OR. N.EQ. 11. OR. N.EQ. 18. OR. N.EQ. 25)  GO  TO  14 
IF(N.EQ.5.0R.N.EQ.12.CR.N.EQ.19.0R.N.EQ.26)  GO  TO  15 
I F( N.EQ. 6. OR. N.EQ. 13.  OR.  N .  E  Q  0  20  .  OR  ,N.  EQ  .  2  7)  GO  TO  16 
IF(N.EQ.7.0K.N.EQ.l4.0R.N.EQ.21.GR.N.EQ.28)  GO  TO  17 

11  WRITE  (6,19)  DAY(l) 
GO  TO  9  9 

12  WRITE  (6,19)  DAY(2) 
GO  TO  99 

13  WRITE  (6,19)  DAY (3) 
GO  TO  99 

14  WRITE  (6,19)  DAY(4) 
GC  TO  99 

15  WRITE  (  6  $  19)  DAY (5) 
GO  TO  99 

16  WRITE  (6,19)  DAY(6) 
GO  TO  99 

17  WRITE  (6,19)  DAY (7) 
GO  TO  99 

19  FORMAT  ( '      TODAY  IS  ' ,A8) 
99  RETURN 
END 


Q  ***  if.if.if.  ***  *  if.******  *  ***  *  *************  *********  **************£ 
C  C 

SUBROUTINE    ICASES    <M,N) 
C  C 

c************ ********* *********  *************** *************c 

c 
c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c  c 

C  THIS  SUBROUTINE  RECEIVES  THE  AVERAGE  NUMBER  C 

C  OF  CASES  WHICH  HAVE  OCCURRED  ON  THIS  DAY  C 

C  AS  DAVG.   FROM  THIS  NUMBER  IT  GENERATES  C 

C  THE  NUMBER  OF  CASES  WHICH  WILL  OCCUR  IN  THE  C 

C  SIMULATION  MODEL.   IT  RETURNS  THIS  NUMBER  C 

C  TO  THE  MAIN  PROGRAM  AS  NCASES.  C 

C  C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

C 

c 

CCMMON  /CGM2/ISEED 
COMMON  /COM11 /NCASES 
CCMMON  /COM12/AVG 
-  DIMENSION  AVG(12t28) 
C 

DAVG  =  AVG(M,N) 
A   =  EXP(-DAVG) 
S  =  1.0 
NCASES  =  0 
10  CALL  RANDOM  (ISEED,U,1) 
S  =  S  *  U 

IF  (S.LT.A)  GO  TO  99 
NCASES  =  NCASES  +  1 
GO  TO  10 
99  WRITE  (6*100}  NCASES 
100  FORMAT  (•      THE  NUMBER  OF  CASES  FOR  THIS  DAY  IS   «, 
112) 
RETURN 
END 
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c 
c 

c  c 

Subroutine  ntime   (ncse) 
c  c 

C 

c 
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


c 

c 

subroutine  ntime  determines  the  times  that 

c 
c 

c 

THE  SAR  INCIDENTS  HAPPEN  ON  ANY  GIVEN  DAY. 

c 

c 

STORED  IN  A  COMMON  STORAGE  LOCATION  IN  THE 

c 

c 

THE  SUBROUTINE  IS  THE  NUMBER  OF  CASES  OR 

c 

c 

INCIDENTS  WHICH  WILL  OCCUR.   THE  SUBROUTINE 

c 

c 

FIRST  DETERMINES  IF  THERE  ARE  MORE  THAN  ONE 

c 

c 

INCIDENTS  FOR  THE  DAY.   IF  NOT,  THE  SUB- 

c 

c 

ROUTINE  GENERATES  ONE  RANDOM  NUMBER,  AND 

c 

c 

COMPARES  THIS  NUMBER  AGAINST  THE  CUMULATIVE 

c 

c 

DISTRIBUTION  STORED  IN  XTIME  TO  DETERMINE 

c 

c 

WHAT  TIME  OF  DAY  THE  CASE  WILL  HAPPEN.  IT 

c 

c 

THEN  GENERATES  A  NUMBER  AND  CONVERTS  THIS 

c 

c 

TO  THE  MINUTE  VALUE.   THE  SUBROUTINE  NEXT 

c 

c 

PRINTS  THE  CASE  NUMBER  FOR  THE  DAY,  THE 

c 

c 

HOUR  AND  THE  MINUTE  OF  THE  CASE, 
CONTINUE 

c 

c 

IF  MORE  THAN  ONE  CASE  WILL  HAPPEN  ON  ANY 

c 

c 

DAY,  THE  SUBROUTINE  GENERATES  AN  EQUIVALENT 

c 

c 

NUMBER  OF  TIMES,  SORTS  THEM  IN  ORDER  OF 

c 

c 

ASCENDING  SEQUENCE,  AND  STORES  THEM  IN  AN 

c 

c 

ARRAY  CALLED  NEXTIM.   AS  THE  MAIN  PROGRAM 

c 

c 

PROCEEDS  THROUGH  ITS  SEQUENCE,  IT  DETERMINES 

c 

c 

HOW  MANY  CASES  WILL  OCCUR  (NCASES),  AND 

c 

c 

THEN  PASSES  TO  THE  SUBROUTINE  THE  CASE 

c 

c 

NUMBER  THAT  IT  IS  CURRENTLY  GENERATING 

c 

c 

(NCSE),   THE  SUBROUTINE  TAKES  THIS  NUMBER, 

c 

c 

DETERMINES  THE  HOUR  OF  OCCURRANCE  FOR  THAT  CASE, 

c 

c 

GENERATES  A  RANDOM  MINUTE  OF  OCCURANCE  WITHIN 

c 

c 

THAT  HOUR,  AND  WRITES  THIS  INFORMATION  ON 

c 

c 

ON  THE  PRINTOUT  BEFORE  RETURNING  TO  THE 

c 

c 

c 

CCX( 

c 

MAIN  PROGRAM. 

c 
c 
c 

XCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

DIMENSION  XTIME (24 ) 

DIMENSION  T(  20)  ,TI.ME(  20)  , NEXTIMt  20i 

CCMMON  /C0M2/ISEED 

CCMMON  /C0M10/XTIME 

COMMON  /COMLI/NCASES 

CCMMON  /COM22/NTYM 

IF(NCASES.GT.l)  GO  TO  30 

CALL  RANDOM  (ISEED,TT,1) 

TYME  =  100.  *  TT 

DO  3.0  I  =  1,24 

IF  (TYME.LE.XTIMEl I) )  GO  TO  20 
10  CONTINUE 
20  NEXTIM ( 1)  =  I  -  1 

GO  TO  71 

30  IF  (NCSE.GT.l )  GO  TO  71 

CALL  RANDOM  ( I SE ED , T , NCASES ) 

DO  60  J  =  I, NCASES 

TIME(J)  =  100.  *  T(J) 

DO  40  I  =  1,24 

IF  (TIME( J) .LE.XTIME( I) )  GO  TO  50 
40  CONTINUE 
5  0  NEXTIM  (J)  =  I  -  1 


60 


I  END  =  NCASES  -  1 
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70 
71 


100 


DO  70  M  =  1, I  END 
JEND  =  NCASES  -  M 


DC  70  N 

IF 

I  TEMP 

NEXTIMt N+l 

NEXTJM(N) 


=  ItJEND 

1M<N).LE .NEXT  I M( N+l) ) 

NEXTiM(N+l) 

)  =  NEXTJM(N) 

=  ITEMP 


GO  TO  70 


CONTINUE 
CCNT1NUE 
L  =  NCSE 
CALL  RANDOM  { 
XMINS  =  60.  * 
IF(NEXTIM(L) . 
WRITE  (o,100) 
FORMAT  (/■ 
1F4.0) 
NTYM  =  NEXTIM(L) 
RETURN 
END 


SEEO,TMtl) 

TM 
EQ*24)  XMINS  =  0 
Lf NEXTIM(L) tXMINS 
CASE  ',12,'  TIME  OF  INCIDENT  IS  «tI2 
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C***  **  ***************  ******  ****  ******  ******  ******  ****  ******q 

c  c 

SUBROUTINE    NATURE 

c  c 

C ************************** ********* * * *  *  * * * * *******  v  *  *  * * *  * * C 

c 
c 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c  c 

C  THE  PURPOSE  OF  THIS  SUBROUTINE  IS  TO  GEN-  C 

C  ERATE  THE  NATURE  OF  DISTRESS  REPRESENTED  BY  C 

C  THE  SAR  INCIDENT  GENERATED  IN  THE  i-i A I N  PRO-  C 

C  GRAM.   EACH  SAR  INCIDENT  OCCURS  MORE  OR  C 

C  LESS  RANDOMLY  OURING  ANY  YEAR.   THE  CUMU-  C 

C  LATIVE  DISTRIBUTIONS  OF  THE  SAR  OCCURRANCES  C 

C  FOR  THE  YEARS  1970  TO  1974  WERE  TABULATED  C 

C  AND  READ  INTO  THE  MAIN  PROGRAM  AS  (HYST).  C 

C  ALONG  WITH  THESE  PERCENTAGES,  THE  DEGREE  C 

C  OF  SEVERITY  CUMULATIVE  PERCENTAGES  FOR  EACH  C 

C  OF  THE  TYPES  OF  DISTRESS  WERE  ALSO  TABU-  C 

C  LATED  AND  READ  INTO  A  STORAGE  ARRAY  (HCUK).  C 

C  THIS  SUBROUTINE  CAN  THEREFORE  GENERATE  A  C 

C  RANDOM  NUMBER  AND  USE  THIS  NUMBER  TO  DETER-  C 

C  MINE  THE  NATURE  OF  THE  SAR  INCIDENT.  C 

C  KNOWING  THE  NATURE  Or  DISTRESS,  THE  PRO-  C 

C  GRAM  GENERATES  ANOTHER  RANDOM  NUMBER  AND  C 

C  USES  THE  SECOND  NUMBER  TO  DETERMINE  THE  C 

C  SEVERITY  OF  THE  DISTRESS.  THE  NATURE  OF  C 

C  DISTRESS  AND  SEVERITY  ARE  RETURNED  TO  THE  C 

C  MAIN  PROGRAM  FOR  FURTHER  USE.  C 

C  C 

c  c 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

C 

COMMON  /CQM2/ I  SEED 

COMMON  / COM3 /HYST 

COMMON  /CQM6/HCUM 

CCMMON  /C0M14/ILAND 

COMMON  /C0M15/IHELCS 

CCMMON  /CCM16/I AIR 

CCMMON  /C0M20/IPER 

CCMMON  /COM25/NODES 

COMMON  /C0M24/IPR0P 
C 

DIMENSION  HYST(23) , IHURTC 16) ,HCUM(23.16) 
C 

CALL  RANDOM  (ISEED,ATI) 
C 

C         COMPARE  THE  RANDOM  NUMBER  AGAINST  THE 
C         CUMULATIVE  DISTRIBUTION  STORED  IN  ARRAY  HYST. 

AINDX  =  100.0  *  A 

DO  10  K  =  1,22 

IF  (AINDX. LE. HYST(K) )  GO  TO  15 
10  CONTINUE 
C         DEPENDING  ON  THE  VALUE  OF  K,  BRANCH  PROGRAM 
C         TO  STATEMENT  REPRESENTING  NATURE  OFD 
C         TO  STATEMENT  REPRESENTING  NATURE  OF  DISTRESS. 
C 
C 

15  IF(K.EQ.l)  GO  TO  199 

IFtK.EU.2)  GO  TO  299 

IFU.EQ.3)  GO  TO  399 

IF(K.EQ*4)  GO  TO  499 

IF(K.EQ.5)  GO  TO  599 

lF(K.Ew.6)  GO  TO  6  99 

IF(K.EQ.7)  GO  TO  799 

IF(K.EQ.8)  GO  TO  899 

IF(K.EQe9)  GO  TO  999 

IF(K.EO.IO)  GO  TO  1099 

IF(K.EQ.ll)  GO  TO  1199 

IF(K.EQ.12)  GO  TO  1299 
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c 
c 
c 
c 
c 
c 


IF(K 
IF(K 
IF(K 
I  F(K 
IF(K 
IF(K 
IF(K 
IF(K 
IF(K 
IF(K 
IF(K 


.EQ 
.EQ, 
.EQ  , 
•  EQ . 

.EQ, 

.  E  Q  , 

.EQ, 

.EQ, 

.  EQ 

.EQ, 

.EQ, 


13)  GO  TO 

14)  GO  TO 
15 J  GO  TO 
16JG0  TO 
17)  GO  TO 

GO 
GO 
GO 
GO 
GO 
GO 


18) 
19) 
20) 
21) 
22  ) 
23) 


TO 
TO 
TO 
TO 
TO 
TO 


1399 
1499 
1599 
1699 
1799 
1399 
1999 
2099 
2199 
2299 
2399 


THIS  SECTION  PRINTS  THE  NATURE  OF  DISTRESS 
AS  WELL  AS  PERFORMING  HOUSEKEEPING  FUNCTIONS 
WHICH  ENSURE  INTERNAL  CONSISTENCY  WITHIN  THE 
MODEL. 


199 

100 


299 
200 


399 

300 


499 

400 


599 
500 


699 

600 


799 

700 


899 

800 


999 

90C 


1099 
1000 


1199 
1100 


WRITE' 6 
FORMAT 

1 'VESSEL 
NODES  = 
GO  TO  9 
WRITE  ( 
FORMAT 

1' AND  DI 
NODES  = 
GO  TO  9 
WRITE  ( 
FORMAT 

1«  AND  AN 
NODES  = 
GO  TO  9 
WRITE  ( 
FORMAT 
I  LAND  = 
NODES  = 
GO  TO  9 
WRITE  ( 
FORMAT 
NODES  = 
GO  TO  9 
WRITE  ( 
FORMAT 

1'  EXPLOS 
NODES  = 
GO  TO  9 
WRITE  ( 
FORMAT 

l'OR  SIN 
NODES  = 
GO  TO  9 
WRITE  ( 
FORMAT 

1«  MISS  IN 
NODES  = 
GO  TO  9 
WRITE  ( 
FORMAT 

l'COLLIS 
NODES  = 
GO  TO  9 
WRITE (6 
FORMAT 

1 'DISORI 
NODES  = 
GO  TO  9 
WRITE  ( 
FORMAT 

1 f  ENOANG 
NODES  = 
GO  TO  9 


,100) 

{•  NATURE 

CONDITIO 

0 
999 
6,200) 
(i  JH5 

SA8L.ED 
I 

999 

6,300) 
(«  THE 
CHORED." ) 

2 
999 

6,400) 
(  «  THE 

2 

999 

6,5  00) 
(•  THE 

4 
999 

6,600) 
('  THE 
ION  ON 

5 
999 
6*700) 
(«  THE 
KING.' 

6 
999 

6,800) 
<•  THE 
G  VESS 

7 
999 

6,900  ) 
('  THE 
ION,5  ) 

8 
999 
,1000) 
('  THE 
ENTATED 

9 
999 

6, 1100) 
('  THE  NA 
ERED  BY  W 

11 
999 


OF  DISTRESS  IS  MISCELLANEOUS  », 
N.  '  ) 


NATURE  OF  DISTRESS  IS  VESSEL  ADRIFT  ', 
»  ) 


NATURE  OF  DISTRESS  IS  VESSEL  DISABLED  », 


NATURE  OF  DISTRESS  IS  VESSEL  AGROUND') 


NATURE  OF  DISTRESS  IS  VESSEL  CAPSIZED.5) 


NA 
VE 


TURE  OF  DISTRESS  IS  FIRE  OR  ', 
SSEL.«  ) 


NATURE  OF  DISTRESS  IS  VESSEL  FLOODING  S 


I 


NA 
EL. 


TURE  OF  DISTRESS  IS  OVERDUE  OR  '? 
•  ) 


NATURE  OF  DISTRESS  IS  VESSEL  «, 


NA 
V 


TURE  OF  DISTRESS  IS  LOST  OR  »  , 
ESSEL. « ) 


TURE  OF  DISTRESS  IS  VESSEL  '  , 
LATHER.' ) 
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1299 
1200 


1399 

1300 


). 


1499 
1400 


1 


1599 
1500 


1 


1699 
1600 


1799 
1700 


18  99 
18  00 


1999 
1900 


2099 
2  0  00 


2199 
2100 


2299 
2200 


2399 
2300 


WRITE  ( 
FORMAT 
•CRASHE 
HAND  = 
NODES  = 
GO  TO  9 
WRITE  ( 
FORMAT 
•CRASHE 
HAND  = 
NODES  = 
GO  TO  9 
WRITE  ( 
FORMAT 
' STRUCT 
I  LAND  = 
NUDES  = 
GO  TO  9 
WRITE  ( 
FORMAT 
•BENDS 
NODES  = 
GO  TO  9 
WRITE  ( 
FORMAT 
•MEDEVA 
NODES  = 
GO  TO  9 
WRITE  ( 
FORMAT 
NODES  = 
GO  TO  9 
WRITE  ( 
FORMAT 
NODES  = 
GO  TO  9 
WRITE  ( 
FORMAT 
! SICKNE 
NODES  = 
CALL  RA 
IF  (Z*L 
GO  TO  9 
WRITE  { 
FOPMAT( 
NODES  = 
CALL  RA 

IF    (Z.L 

GO  TO  9 
WRITE  ( 
FORMAT 
•  LAND. ' 
I  LAND  = 
NODES  = 
GO  TO  9 
WRITE  ( 
FORMAT 
.'  MISCEL 
NODES  = 
GO  TO  9 
WRITE  ( 
FORMA T( 
NODES  = 


6, 1200) 
('  THE  NA 
D  OR  LOST 

I 

20 
999 

6, 1300) 
(•  THE  NA 
D  OR  LOST 

1 

40 
999 

6, 1400) 
(  '  THE 
URE.« ) 

1 

50 
999 

6, 1500) 
(«  THE  NA 

OR  LOST. 

70 
999 

fa, 1600) 
(•  THE 
C«  ) 

80 
999 

6,  1700) 
(  •  THE 

81 
999 

6,  1800) 
{•  THE 

82 
999 

6,1900) 
(»  THE 
SS.«  ) 

84 
NDOM  (ISE 
T.0.5)  IL 
999 

6,2000) 
•  THE  NAT 

85 
NDOM  (ISE 
5)  IL 


TURE  OF  DISTRESS 
.'  ) 


TURE  OF  DISTRESS 
.»  ) 


IS  AIRCRAFT  ' , 


IS  LAND  VEHICLE  ' , 


NATURE  OF  DISTRESS  INVOLVES  LAND  «  , 


TURE  OF  DISTRESS 
') 


IS  DIVER  WITH  •  , 


NATURE  OF  DISTRESS  IS  PERSONNEL  «, 


NATURE  OF  DISTRESS  IS  PERSON  DROWNING.') 


NATURE  OF  DISTRESS  IS  PERSON  IN  WATER.8) 


NATURE  OF  DISTRESS  IS  PERSONNEL  ', 


ED,Z,1 ) 
AND  -  1 


URE  OF  DISTRESS  IS  PERSONNEL  INJURY, «) 

EDi  Z, 1 ) 

AND  =  1 


NATURE  OF  DISTRESS  IS  PERSON  LOST  ON  «, 


T.O 

999 

6,2100) 
(!  THE 
) 

1 

87 
999 

6,2200) 
(•  THE 
LAN ECUS 

90 
999  . 
6,2300) 
1  THE  NATURE  OF  DISTRESS  IS  FLARE  SIGHTING.') 

98 


NA 
C 


TURE  OF  DISTRESS 
ONDITION.1 ) 


IS  OTHER  ■  , 
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c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c  c 

C  THIS    PORTION    OP    THE    SUBROUTINE    COMPUTES    THE  C 

C  SEVERITY    OF    THE    SAR    INCIDENT.       THE    FOLLOWING  C 

C  LABELS    ARE    USED.     .     .  C 

C  c 

C  I  PER   =   THE  SEVERITY  TO  PERSONNEL  (0  TO  3)  C 

C  I  PROP   =   THE  SEVERITY  IT  PROPERTY  (0  TO  3)  C 

C  IHURT(I)   =   THE  DESIGNATOR  OF  THE  GRID  C 

C  USED  TO  COMPUTE  I  PER  AND  I  PROP.  C 

c  c 

C         CODE  0   =   NO  SEVERITY  TO  PERSONS  OR  PROPEkTY      C 
C         CODE  1   =   SMALL  SEVERITY  (LITTLE  DANGER)  C 

C  CODE  2   =   MODERATE  DANGER.  MIGHT  USE  HELO.        C 

C         CODE  3   =   GREAT  DANGER*   DEFINITELY  USE  HELO*     C 
C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 
c 

9999  DO  20  I  =  1,16 

I  HURT (I)  =  I 
20  CONTINUE 

CALL  RANDOM (  ISEED,P,1) 

HURT  =  100.  *  P 

DO  30  J  =  1,16 

IF  (HURT.LE.HCUM(  K,  J)  )  GO  TO  35 
30  CONTINUE 
35  I  PER  =  (I  HURT! J I  -  1)  -  ((IHURT(J)  -  l)/4  *  4) 

I  FROP  =  (IHURTU  )  -  l)/4 

WRITE  (6,3G00)  I  PER,  I  PROP 
3000  FORMAT  (!  SEVERITY  TO  PERSONNEL  IS  B , 1 2 , 5X , 
Is  SEVERITY  TO  PROPERTY  IS  ',12) 

IFt ( IPER.LT.2) .AND. (  IPR0P.LT.2) }     GO  TO  50 

WRITE  (6,4000) 
4000  FORMAT? »  HELICOPTER  WOULD  BE  USED  ON  THIS  CASE.') 

IHELCS  =  IHELCS  +  1 

I  AIR  =  1 
50  CONTINUE 

RETURN 

END 
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Q********!!;*^^^*^)):^^:^*  $  $  $  $  $  $  *  *  *  *  *  *  #  *  *  #  *  *  #  *  V  *  *  *  *  *  *  *  ♦  *  *  *  *  *  #  C 

c  c 

SUBROUTINE    LATLOiM    (1FLAG) 
C  C 

Q*#  *#*  *  $$##*$>{($  $$£$£$*  $*$$$$ £$**$$  $$$  **  &*$$**  $  *  *  £  *  *  ***$  #  *  *  *£ 

c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c  c 

c  c 

C  THE  LATLON  SU3RGUTINE  RECEIVES  THE  SEED  C 

C  VALUE  FROM  THE  MAIN  PROGRAM  TO  GENERATE  A  C 

C  RANDOM  NUMBER.  IT  USES  THIS  VALUE  TO  COM-  C 

C  PUTE  THE  LATITUDE  AND  LONGITUDE  OF  THE  SAR  C 

C  INCIDENT.   THE  SUBROUTINE  IS  USABLE  FOR  C 

C  LATITUDES  41  AND  42N  AND  LONGITUDES  124-00  C 

C  TO  125-OOW.   FOR  THE  FEW  SAR  INCIDENTS  THAT  C 

C  OCCUR  OUTSIDE  THIS  BOUNDARY,  A  SEPARATE  C 

C  PROCEDURE  IS  USED,   THE  LATITUDE  AND  LGN-  C 

C  GITUDE  ARE  COMPUTED  USING  DIOPHANTINE  EQUA-  C 

C  TIONS.   IN  THE  SUBROUTINE,  THE  FOLLOWING  C 

C  VARIABLES  ARE  USED.  .  .  C 

c  c 

C  LAT   =   THE  LATITUDE  OF  THE  SAR  INCIDENT*  C 

C  LONG   =   THE  LONGITUDE  OF  THE  SAR  INCIDENT.  C 

C  MTN   =   THE  MINUTES  OF  LATITUDE.  C 

C  LMINS   =   THE  MINUTES  OF  LONGITUDE.  C 

C  C 

c  c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

CCMMON  /C0M1/XCUM 

CCMMON  /C0M2/ISEED 

CCMMON  /C0M7/LAT 

COMMON  /C0M8/MINS 

CCMMON  /C0M9/LMINS 

CCMMON  /C0M13/DIST 

COMMON  /CCM14/IL AND 

CCMMON  /C0M20/  I  PER 

COMMON  /C0M21/  MEW 

COMMON  /C0M2  4/  I  PROP 

CCMMON  /COM2  5/LLIND 

CCMMON  /COM? 6/  ODD  I  ST 
C 

DIMENSION  LLIND(2  88J , XCUM<2  88) 
C 

C         IF  THE  INDICATOR  I  LAND  IS  EQUAL  TO  ONE,  THE 
C         SAR  INCIDENT  OCCURS  INLAND  AND  THE  PROCEDURE 
C         FOR  GENERATING  THE  LATITUDE  AND  LONGITUDE 
C         OF  THE  DISTRESS  IN  AN  OFFSHORE  POSITION  IS 
C         BYPASSED. 
C 

DDDIST  =  0. 

IF  (ILAND.EQ.l)  GO  TO  99 
C 
C 

CALL  RANDOM  i I SE ED , XSQR , 1 J 

SQR  =  100.  *  XSQR  ♦  1. 
C 

C  THE  RANDOM  NUMBER  IS  COMPARED  WITH  THE 

C         ARRAY  CONTAINING  THE  HISTORICAL  DISTRIBUTIONS 
C         OF  SAR  INCIDENTS. 


C 


DO  50  N  =  1,288 
IF (SQR. LE.XCUM(N J )  GO  TO  60 
50  CCNTINUE 


C 

C  IF  THE  RANDOM  NUMBER  IS  GREATER  THAN  THE 

C  DISTRIBUTIONS  IN  THE  ARRAY,  THE  INCIDENT 

C  OCCURRED  EITHER  INLAND  OR  FARTHER  WEST  THAN 

C  125-OOW.   IN  THAT  CASE,  POSITION  IS  COMPUTED 

C  IN  A  DIFFERENT  MANNER. 

C 


123 


IF(SQR.GT.XCUM(N) ) GO  TO  70 
C 

C         LATITUDE  AND  LONGITUDE  ARE  COMPUTED  USING 
C  THE  VALUE  OF  THE  SQUARE  IN  THE  ARRAY  WHICH 

C         CONTAINED  THE  RANDOM  NUMBER, 


C 


60  LONG  =  124 

I  NO  =  LLIND(N) 


LMINS  =  ((UNO  -  i>/24)  +  1)  *  5 
C 

C         IF  THE  INDICATOR  ILAND  HAS  BEEN  SET  EQUAL  TO 

C  TO  2,  IT  INDICATES  THAT  THE  SAR  INCIDENT 

C  MUST  OCCUR  ADJACENT  TO  THE  SHRRE,  SUCH  AS 

C         OCCURS  WHEN  A  VESSEL  GOES  AGROUND.   IN  THAT 

C         CASE,  THE  MINUTES  OF  LONGITUDE  ARE  CONSTRAINED 

C         SO  THE  POSITION  WILL  BE  NEXT  TO  THE  SHORE. 

C 

IF  ( ILAND. EQ. 2)  LMINS  =  10 

LAT  =  40 

ILAT  =   ((IND-1)  -   (UNO  -  U/24)*24)  +  I 

IF( ILAT .GT.12)  LAT  -  41 

IFCILAT. GT.12)  ILAT  =  ILAT  -  12 

MINS  =  5  *  ILAT 

WRITE  (6,100)  LAT, MINS, LONG, LMINS 
100  FORMAT  (■  SAR  INCIDENT  OCCURS  LATITUDE  »,2I3, 
1»  MINS  NORTH,  LONGITUDE  '^IS,8  MINS  WEST.1) 

GC  TO  999 
C 
C 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


C 
C 
C 

IF  THE  INDICATOR  SQR  IS  GREATER  THAN  THE 

C 

c 
c 

C 

CUMULATIVE  PROBABILITIES  STORED  IN  XCUM 

c 

C 

(THE  MAXIMUM  OF  WHICH  IS  94.4)  IT    MEANS 

c 

c 

THAT  THE  SAR  INCIDENT  OCCURRED  OUTSIDE  THE 

c 

c 

LONGITUDE  PARALLELS  OF  124-OQW  AND  L24-59W 

c 

c 

THIS  OCCURANCE  IS  SIMULATED  IN  A  RANDOM 

c 

c 

MANNER  BECAUSE  OF  THE  ABSENCE  OF  ADEQUATE 

c 

c 

HISTORICAL  DATA  SHOWING  THE  CUMULATIVE 

c 

c 

OCCURaNCES.   THE  LATITUDE  IS  RANDOMLY 

c 

c 

GENERATED  ASSUMING  EQUAL  PROBABILITIES 

c 

c 

CAN  CCCUR  FOR  ANY  LATITUDE  BETWEEN  40-00N 

c 

c 

AND  41-OON.   THE  DISTANCE  FROM  EUREKA  IS 

c 

c 

ALSO  RANDOMLY  GENERATED  AND  THE  DISTANCE 

c 

c 
c 

TO  SCENE  SUBROUTINE  (DISTNC)  IS  BYPASSED. 

c 
c 

c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

c 
cc 

70  CALL  RANDOM*  ISEED, XL, 1) 

XLATN  =  12.  *  XL  +  1. 

LLATN  =  XLATN 

IF(LLATN.GT.6)  GO  TO  75 

LAT  =  40 

MINS  =  5  *  LLATN 

GO  TO  80 
75  LLATN  =  12  -  LLATN 

LAT  =  41 

MINS  =  5  *  LLATN 
80  CALL  RANDOM  USEED,D,1) 

DIST  =  300.  *  D 
C 

C         THE  DISTANCE  OFFSHORE  MUST  BE  GREATER  THAN 
C  50  MILES  Ok  IT  WOULD  HAVE  OCCURRED  WITHIN 

C         THE  AREA  COVERED  BY  THE  GRID  XCUM. 


C 

C 

C 

c 


IF  (DIST-LT.50,)  GO  TO  80 


124 


WRITE  (6,200)  LAT,MINS,DIST 
200  FORMAT*/'  SAR  INCIDENT  OCCURS  LATITUDE  't     213, 

1«  KINS  NORTH.  DISTANCE  TO  SCENE  IS  • ,F5 . 0» ■ MILES .* ) 
LMINS  =  0 
IFLAG  =  1 

IF  ((DIST. GE. 50. J . AND.CDIST.LT. 100.) )  MEW  =  5 
IF  ( (DIST. GE. 100.) .AND. (DIST, LT. 150.) )  MEW  =6 
IF   CDIST.CE .150. )   MEW  =  7 
DDDIST  =  DIST 

IF  (  C  IPER.LT.3)  ,  AND.  (  I  PROP.LT.  :•»  i  )  GO  TO  90 
C 

C  IF  THE  DISTANCE  OFFSHORE  IS  GREATER  THAN 

C         THE  CAPABILITY  OF  THE  HH-52  HELICOPTER* 
C         AND  THE  PROGRAM  HAS  DECIDED  THAT  AN  AIR- 
C         CRAFT  IS  REQUIRED  FOR  THE  CASE  (SEVERITY 
C         CODE  GREATER  THAN  2),  THE  PROGRAM  COMPUTES 
C  TIME  FOR  A  HELICOPTER  FROM  SAN  DIEGO  TO 

C         ARRIVE  ON  SCENE. 
C 

WRITE  C6,300) 
300  FORMAT  <<  DISTANCE  TO  SCENE  IS  BEYOND  HH-521  , 
1 •  CAPABILITY.' ) 
ETIME  =  6.  +■  DIST  /  120. 
WRITE  (6,400)  ETIME 
400  FORMAT  («  TIME  REQUIRED  FOR  HH-3F  HELICOPTER6, 
1*  TO  ARRIVE  IS  9,F3.i.  ,'  HOURS-1) 
GO  TO  999 
C 

C         TIME  FOR  A  LARGE  UNIT  FROM  HUMBOLT  BAY  IS 
C         COMPUTED. 
C 

90  ETIME  =  DIST  /  18. 
WRITE  (6,500)  ETIME 
500  FORMAT  C  TIME  FOR  A  SURFACE  UNIT  FROM  ', 
1« HUMBOLT  BAY  TO  ARRIVE  ON  SCENE  IS  ctF4.1f 
2«  HOURS.') 
GO  TO  999 
C 
C 

C         IN  ORDER  TO  GENERATE  A  POSITION  OF  DISTRESS 
C  INLAND  A  BEARING  TO  THE  SCENE  FROM  THE  AIR 

C  STATION  IS  GENERATED  AND  THEN  A  RANDOM 

C         RANGE  IS  GENERATED 
C 
C 

99  CALL  RANDOM  (ISEED,B,1) 
BRNG  =  180.  *  B 
IBRNG  =  BRNG 
CALL  RANDOM  (ISEED,R,1) 
RANGE  =  100.  *  R 
DIST  =  RANGE 
IRNG  =  DIST 

WRITE  (6,600)  IBRNG, IRNG 
600  FORMAT  ('  INLAND  SAR  INCIDENT  OCCURS  BEARING  • , 

113,'  AT  A  DISTANCE  OF  ',13,*     MILES  FROM  AIR  STATION.') 
IFLAG  =  1 
MEW  =  0 
LAT  =  0 
M1NS  =  0 
LMINS  =  0 
DDDIST  =  DIST 
ETIME  =  DDDIST  /  85. 
WRITE  (6,700)  ETIME 
700  FORMAT  (•  TIME  FOR  HELICOPTER  FROM  ARCATA  TO 

l'ON  SCENE  IS  ',F4.1,'   HOURS.') 
999  RETURN 
END 
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C  C 

SUBROUTINE    DISTNC 

C  C 

C  *  *  *  *  *  *  *  *  *  *  *  #  *  *  *  *  *  *  3?  *  *  #  *  *  *  *  *  *  *  *  ♦  *  *  *  *  *  *  *  *  *  *  *  ♦  *SSt####**#*##*#*C 

c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c  c 

c  c 

C         THIS  SUBROUTINE  COMPUTES  THE  DISTANCE  FROM  C 

C          THE  PROPOSED  AIR  STATION  TO  THE  POSITION  C 

C          OF  THE  SAR  INCIDENT  GENERATED  BY  SUBROUTINE  C 

C          LATLON.   THE  FOLLOWING  LABELS  ARE  USED.  .  .  C 

C  C 

C         XNS  =  THE  DIFFERENCE  IN  DEGREES  OF  LATI-  C 

C               TUDE  BETWEEN  THE  AIR  STATION  AND  THE  SAR  C 

C               INCIDENT.   EACH  MINUTE  OF  LATITUDE  IS  C 

C              EQUIVALENT  TO  ONE  NAUTICAL  MILE.  C 

C          XEW  9    THE  DISTANCE  TO  THE  SAR  INCIDENT  C 

C               ALONG  A  PARALLEL  OF  LATITUDE.   THE  C 
C               DEPARTURE,  OR  THE  DISTANCE  BETWEEN  DEGREES    C 

C               OF  LONGITUDE,  IS  TAKEN  FROM  THE  LATITUDE  C 

C              OF  +1  DEGREES  NORTH  AND  IS  45. 2d  MILES.  C 

C  C 

C         IT  IS  ASSUMED  THAT  IF  THE  DISTANCE  OFF  -  C 

C          SHORE  (XEW)  IS  GREATER  THAN  25  MILES  THE  C 

C          HELICOPTER  WOULD  FLY  PARALLEL  TO  THE  COAST  C 

C         AND  THEN  PROCEED  DIRECTLY  OFFSHORE  WHEN  C 

C          AT  THE  COREECT  LATITUDE.   THIS  ALTERS  THE  C 

C         METHOD  OF  COMPUTING  THE  DISTANCE  TO  THE  C 

C         SCENE  OF  THE  SAR  INCIDENT.  C 

c  c 

C         THIS  ALGORITHM  IS  RESTRICTED  TO  LATITUDES  C 

C         40  AND  41N  AND  LONGITUDE  124W,  C 
CONTINUE 

c  c 

C                       THE    POSITION    OF    EUREKA    AIR    STATION    IS:  C 

C                                                  41-00N          124-05W  C 

C  c 

C         THE  POSITION  OF  GROUP  HUMBGLT  BAY  MOORING  IS*.  C 

C                    40-47N    124-10W  C 

C  C 

C         THE  POSITION  OF  CRESCENT  CITY  MOORING  IS:  C 

C                    4I-45N    124-10W  C 

C  c 

c  c 

c  c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

CCMMON  /C0M7/LAT 

COMMON  /C0M8/MINS 

CCMMON  /C0M9/LMINS 

CCMMON  /C0M13/DIST 

CCMMON  /C0M14/ILAND 

CCMMON  /C0M16/IAIR 

CCMMON  /C0M21/MEW 

MM  INS  =  0 
C 

IF  (IAIR.EQ.O)  GO  TO  30 
C 
C         COMPUTE  DISTANCE  FOR  HELICOPTER  FROM  AIR  STA 


C 


IF  (LAT.EG.40)  XNS  =  60.  -  MINS 

XMINS  =  LMINS 

XEW  =  45.28  *  A6S(XMINS  -  10. 01/60. 0 

IF  (XEW. GT. 25.0)  GO  TO  10 

DIST  =  ((XNS  **  2.)  *-  (XEW  **  2.})  **  0.5 

GC  TO  20 
10  DIST  =  XNS  ♦■  XEW 
20  WRITE  ( 6, 100)  DIST 
100  FORMAT (■  DISTANCE  FROM  AIR  STA  TO  POSITION', 
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c 
c 
c 
c 

c 
c 

c 

c 


c 
c 
c 
c 


1"     OF     SAR     INCIDENT     IS       ',F8.2,     '     MILES.') 
ETIME    =    DIST/S5.0 
WRITE     (6,150)     ETIME 
150    FORMAT     (•     ENROUTE    TIME    FOR    HELICOPTER    TO    «, 
l'SCENE    CF    DISTRESS    IS    »,F3.1,     '     HOURS.') 

DECIDE    ON    CLOSEST    SURFACE    UNIT. 

30    IF    ( (LAT.EQ.41) .AND. (MiNS.GT.3G)  )    GO    TO    50 

COMPUTE  DISTANCE  TO  SCENE  FOR  SURFACE  UNIT 
FROM  HUMBOLT  BAY. 

XNS  =  0. 

IF  (  (LAT.EQ.40)  .AND.  (  MINS. GT. 47)  )  XNS  =  MINS  ■ 

IF  (XiMS.GT.O)  GO  TO  40 

XNS  =  47.  -  MINS 
40  XMINS  =  LMINS 

XEW  =  45.28  *  (XMINS  -  10.J/60.0 

DIST  =  ((XNS  **  2.)  *  (XEW  **  2.))  **  0.5 

WRITE  (6,200)  DIST 
200  FORMAT  ('  THE  DISTANCE  TO  SCENE  FOR  A  ', 
1' SURFACE  UNIT  FROM  HUMBOLT  BAY  IS  B,F5.1) 

ST  I  ME  =  DIST/I80O 

WRITE  (6,250)  ST  I  ME 
250  FORMAT  (•  THE  TIME  REQUIRED  FOR  A  SURFACE  «, 
1*UNIT  TO  ARRIVE  ON  SCENE  IS  «, F3.lt1  HOURS.1) 

GO  TO  60 


COMPUTE  DISTANCE  JO    SCENE  FROM  CRESCENT  CITY. 


47 


50  I 
X 
G 

55  X 

56  X 
X 


300 


D 
W 
F 
1« 
S 

w 

60  W 
400  F 
80  C 
I 
I 
1 
I 
I 
I 
I 
R 
E 


F  (MIN 

NS  =  M 

0  TO  5 

NS  =  4 

MINS  = 

EW  =  4 

1ST  = 

RITE 

ORMA 

CITY 

TIME 

RITE 

RITE 

CRMAT( 

ONTINU 


(  (XE 
(  (XE 
((XE 
(  (XE 
(  (XE 
((XE 
( X  E  W 

:ETURN 

ND 


Sc 

IN 

6 

5. 

L 
5. 
(( 
6, 
(' 
S 

D 

e, 

6 
1 

f- 

W 

w 

w 
W 
W 
W 


LT 


45)  GO  TO  55 
45. 


0.5 


-  MINS 
MINS 

2  8  *  ABS( XMINS  -  10.0) /60.0 
XNS  **  2. )  +     (XEW  **  2. ) )  *: 
300)  DISV 

THE  DISTANCE  TO  SCENE  FROM  CRESCENT  «, 
*  ,F5.1 5 •  MILES.'  ) 
I  ST/ 18.0 
250}  STIME 
400)XEW 
THE  DISTANCE  OFFSHORE  IS  »,F5.1,«  MILES8) 

GE.O.)  .AND.(XEW.LE.3.M  MEW  =  1 
GT.3. ) .AND.tXEW.LE.10.) )  MEW  =  2 


GT.10. ) .AND. (XEW.LE.2  5. ) )  MEW  =  3 
GT.25o  )  .AND.  (  XEW.LE.50.  )  )  MEW  =  4 
GT.50.) .AND. (XEW.LE.100.) )  MEW  =  5 
GT.100.  )  .AND. (XEW. LE. 150. ))  MEW  =  6 
iT.150.  )  MEW  =  7 
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q ****** ******  *******************************#**************£ 

C  C 

SLBROUTINE  WEATHR 
C  C 

C *********************************************** ***********£ 

c 
c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c  c 

c  c 

c       this  subroutine  computes  the  weather  con-       c 

C  DITIONS  AT  THE  SCENE  OF  THE  SAR  INCIDENT  C 

C  WHICH  WAS  GENERATED  IN  THE  MAIN  PROGRAM.  C 

C  THE  CUMULATIVE  DISTRIBUTIONS  OF  THE  WEA-  C 

C  THER  CONDITIONS  ARE  BASED  UPON  THE  HIS-  C 

C  TORICAL  WEATHER  AT  ALL  SCENES  OF  SAR  C 

C  INCIDENTS  FOR  THE  YEARS  1970  TO  1.974.  C 

c  c 

C  THE  VARIABLE  LABELS  USED  IN  THIS  SUBROUTINE       C 

C         ARE  AS  FOLLOWS.  .  .  C 

C  c 

C  ISTATE   =   THE  CELL  NUMBER  FOR  THE  WEATHER  C 

C                    CONDITIONS.   IT  IS  DERIVED  BY  C 

C                     GENERATING  A  RANDOM  NUMBER  AND  C 

C                     COMPARING  THIS  NUMBER  WITH  THE  C 

C                     HISTORICAL  CUMULATIVE  PERCENTAGES  C 

C                     FOR  THE  WEATHER  CONDITIONS.  C 

C  IWIND   =   THE  WIND  STATE  CODE  (0  TO  6).  C 

C  ISEA   =   THE  SEA  STATE  CODE  (0  TO  8).  C 

C  I  VIZ   =   THE  VISIBILITY  STATE  CODE  {0  TO  3).  C 

C  C 

C  C 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

DIMENSION  ISTATE  (729),  WXNR(729) 

c 

CCMMON  /C0M2/ISEED 

CCMMON  /C0M3/I STATE 

COMMON  /C.0M4/WXNR 

CCMMON    /COM] 4/ I  LAND 

CCMMON    /C0M17/    ISEA 

CCMMON    /C0M18/IWIND 

COMMON    /CQM19/IVIZ 
C 

5    CALL    RANDOM! ISEED,W,1 ) 

WX    =    100.    *    W 

DO    20    I     =    1,729 

IF(WX.LE.WXNR( li  )     GO    TO    2  5 
20    CONTINUE 
C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
C  c 

C         THE  GRID  FOR  THE  WEATHER  CONDITIONS  AT  THE         C 
C         THE  SCENE  OF  THE  SAR  INCIDENT  HAS  THE  C 

C         FOLLOWING   DIMENSIONS.  .  .  C 

C  c 

C         SEA  STATE  (ISEA)   X-AXIS  WITH  DIMENSION  9          C 

C                            (0  TO  8)  C 

C         WIND  (IWIND)  Y-AXIS  WITH  DIMENSION  9  (0  TO  3)      C 

C          VISIBILITY  (I  VI*)  Z-AXIS  WITH  DIMENSION  9          C 
CONTINUE 
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c  c 

C  USING  DIUPHANTINE  EQUATIONS,  WE  CAN  DERIVE  C 

C  THE  CODE  FUR  THE  WIND,  SEA,  AND  VISIBILITY  C 

C  KNOWING  THE  CELL  NUMBER  (ISTATEJ  REPRESENTING  C 

C  THE  JUNCTION  OF  THE  THREE  CODES.   THIS  C 

C  CELL  NUMBER  IS  GENERATED  RANDOMLY  WHEN  "W"  C 

C  IS  GENERATED  ABOVE.   THE  GENERAL  EQUATION  C 

C  FOR  THIS  CELL  NUMBER  GIVEN  THE  VALUES  OF  C 

C  THE  THREE  AXIS  IS  "L   -   NM(K-l)  +  N(J-l)  C 

C  +  I  WHERE  IN  THIS  EQUATION  L  IS  ISTATE(I),  C 

C  N  =  9  (THE  NUMBER  OF  ROWS),  M  =  9  (THE  C 

C  NUMBER  OF  COLUMNS)  AND  K  IS  THE  UZ"  AXIS  C 

C  VARIABLE,  J  IS  THE  "XH  AXIS  VARIABLE,  AND  C 

C  I  IS  THE  "Y"  AXIS  VARIABLE.   THE  BASIC  C 

C  EQUATIONS  TO  DETERMINE  EACH  OF  THESE  VALUES  C 

C  ARE  AS  FOLLOWS.  .  .  C 

CONTINUE 

C  C 

C  J   =   (L-NM(K-l)-l  )/N  *■  1  C 

C  K   =   (L-l)NM  +  i  C 

C  I   =   (L-NM(K-I)-l J  -  ( ( L-NM(K-l)-i)/N*N)+l  C 

C  C 

C  THESE  EQUATIONS  MUST  BE  ALTERED  IN  THE  C 

C  MODEL  FOR  THE  AIR  STATION  SIMULATION  BE-  C 

C  CAUSE  SOME  OF  THE  ALLOWABLE  CODES  ARE  C 

C  ZERO  AND  THE  BASIC  EgUATIONS  GIVE  VALUES  C 

C  GREATER  THAN  ONE.   IN  ORDER  TO  TAKE  INTO  C 

C  ACCOUNT  THE  FACT  THAT  ZEROES  ARE  PRESENT,  C 

C  THE  »+l"  FACTOR  IS  DROPPED  FROM  EACH  EQUA-  C 

C  TION  AND  THE  •'  K- 1 "  FACTOR  FOR  "J"  AND  C 

C  "I"  EQUATIONS  IS  CHANGED  TO  SIMPLY  "K".  C 

C  c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

25    NM    =    81 

N    =    9 

c 

I  VIZ    =    (  I  STATE! I  )    -    1J/NM 

Ilr.IND    =    (ISTATE(I)    -    NM*(  IVIZ)-1)/N 

ISEA=(  ISTATEt  I)-NM*(  IVIZJ-1  )-(  (  i  ISTATE(  I  )-NM*  (  I V  I  Z  }-  1. } 
1/N)*N) 

IF  (ISEA.EQ.8)  GO  TO  5 

IF<  ILAND.EQ.l)  ISEA  =  8 
C 

WRITE(6,2G0)  IV  I Z , I  WIND, I  SEA 
200  FORMAT  (•  THE  VISIBILITY  ON  SCENE  IS  CODE  ',12, 
1«  THE  WIND  ON  SCENE  IS  CODE  «,I2r/ 
2'  THE  SEA  STATE  ON  SCENE  IS  CODE  «,I2) 

RETURN 

END 
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c 
c 
c 

C  A  PROGRAM  TO  COMPUTE  RESPONSE  TIMES 

C 

C  USING 

C 

C  COMPUTER- SIMULATION  DATA 

C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccc.cccc 
c  c 

c  c 

C  SFRAN  AIR  STATION  37-37N  122-20W  C 

C  ARC ATA  AIR  STATION  41-OON  124-05W  C 

C  HUMbOLT  BAY  MOORINGS  40-47N  124-10W  C 

C  CRESCENT  CITY  MOORINGS  41-45N  124-10W  C 

c  c 

c  c 

CONTINUE 

c  c 

C  THIS  PROGRAM  COMPUTES  THE  RESPONSE  TIMES  FOR  C 

C  COAST  GUARD  UNITS  TO  ARRIVE  ON  SCENE  FOR  C 

C  SEARCH  AND  RESCUE  INCIDENTS-   THE  POSITION  C 

C  OF  THE  INCIDENTS  AND  THE  SEVERITY  IS  READ  C 

C  FROM  PUNCHED  CARDS.   IF  THE  INCIDENT  OCCURS  C 

C  ON  LAND,  THE  PROGRAM  BRANCHES  AND  DOES  NOT  C 

C  CONSIDER  RESPONSE  BY  SURFACE  UNITS.   THE  C 

C  PROGRAM  DETERMINES  RESPONSE  TIMES  FROM  THE  C 

C  PROPOSED  COAST  GUARD  AIR  STATION  AT  ARCATA  AS  C 

C  WELL  AS  THE  PRESENT  AIRBORNE  CAPABILITIES  OF  C 

C  THE  COAST  GUARD  IN  SAN  FRANCISCO.   AS  AN  C 

C  ALTERNATE  POSSIBILITY,  THE  RESPONSE  TIMES  C 

C  FOR  HH-3F  HELICOPTERS  LOCATED  IN  SAN  FRANCISCO  C 

C  ARE  ALSO  CALCULATED.  C 

CONTINUE 

C  THE  DECISION  CRITERIA  USED  IS  THAT  A  HELI-  C 

C  COPTER  FROM  ARCATA  WOULD  RESPOND  TO  A  DISTRESS  C 

C  INCIDENT  IF  THE  REPORTED  SEVERITY  WAS  MODERATE  C 

C  OR  SEVERE  (CODE  2  OR  3i,   A  HELICOPTER  C 

C  WOULD  RESPOND  FROM  SAN  FRANCISCO  ONLY  IF  THE  C 

C  SEVERITY  WAS  CODE  3.   T;ir-  PROGRAM  Di  FERMINES  C 

C  RESPONSE  TIME  FOR  EACH  FACILITY  ACCORDING  TO  C 

C  THESE  CRITERIA  AND  PRINTS  THEM.   IT  THEN  C 

C  DETERMINES  MINIMUM  RESPONSE  TIME  FOR  A  DISTRESS  C 

C  WITH  AND  WITHOUT  THE  PRESENCE  OF  THE  C 

C  PROPOSED  COAST  GUARD  AIR  STATION  AT  ARCATA.  C 

C  c 

c  c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

WRITE  (6,400) 
400  FORMAT  (///•  COAST  GUARD  ', 

1»  RESPONSE  TIMES  —  SIMULATED  DATA' ///I 
WRITE  (6,300) 
300  FORMAT  (<  CASE       HUMBOLT  BAY       CRESCT  CITY1, 
1«         ARCATA  AIR     SFRAN  AIR     MINIMUM     OLD', 
2'  MINIMUM     LAND'//) 
C 

C         TOTALS  ARE  SET  TO  ZERO. 
C 

TCT  =  0. 
TCTNEW  =  0. 
TCTH3  =  0. 
T0THH3  =  0. 
NRLAND  =  0 
NRHH3  =  0 
NRFAR  =  0 
TCT2  =  0. 
TGT3  =  0. 
TLTNU2  =  0. 
T0TNU3  =  0. 
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c 

c 


c 
c 
c 
c 
c 

c 
c 
c 


c 
c 


c 
c 
c 
c 
c 


c 
c 
c 
c 
c 


c 
c 
c 
c 


20 


XMCD  =  0. 
XSEV  =  0. 


0VER1  =  0. 
0VER5  =  0. 
GREAT5  =  0. 
GREAT1  =  0. 
M0RE5  =  0 
MCRE1  =  0 
CCAST  =  0, 


COMPUTATION  BEGINS. 


DO  990  J  =  1, NUMB 


TOTA 

LS 

LAND  = 

0 

HDIST  = 

0. 

HTIME  = 

0. 

CDIST  = 

0. 

CTIME  = 

0. 

ADIST  = 

0. 

ATI  ME  = 

0. 

SFDIST 

=  0 

SFTIME 

=  0 

PTIME  = 

0. 

CEW  =  0 

e 

A  El,  =  0 

6 

SFEW  = 

0. 

HEW  =  0 

• 

UTIME  = 

0. 

ARE  SET  TO  ZERO  FOR  EACH  CASE 


10  READ  (5,100)  PLAT, PMI  N  S,  PL  M  INS  , IPEP, IPROP, ISEA, DDDIST 
100  FORMAT  (34X,2F2.0,3X,F2.0, 1X,2I1 ,1X,1 l,5X*F3cO) 

CHECK  FOR  LAND  CASE  BY  EXAMINING  ISEA  VALUE. 
IF  LAND  CASE,  BRANCH  TO  AIR  STATION  COMPU- 
TATIONS.  INCREMENT  LAND  CASE  COUNTER. 

IF  (ISEA.EQ-3)  NRLAND  =  NRLAND  +  1 
IF  ( ISEA.EQ.8)  GO  TO  40 

IF  DISTANCE  OFFSHORE  IS  GREATER  THAN  100  MILES, 
INCREMENT  COUNTER  BY  ONE.   THIS  CASE  WOULD  BE 
HANDLED  BY  A  HH-3F. 

IF  (DDDIST. GT. 100.)  NRHH3  =  NRHH3  -5-  1 
HEW  =  DDDIST 


COMPUTE  DISTANCE  TO  SCENE  FROM  HUMBOLT  BAY  MOORINGS 


(PLAT. EQ. 41 . ) 
(PLAT. EC. 41 . ) 
(PMINSeGT.47. I 
( PMINS.LE.47. ) 
(PLMINS.GT.IO. 
(PLMINS .LE. 10. 
HDIST  =  SQRT(HNS 
IF  (DDDIST. GT. 0. ) 
HTIME  =  HDIST  /  18. 


IF 
IF 
IF 
IF 
IF 
IF 


HNS  = 
GO  TO 

HNS  - 


PMINS  *  13. 
20 

PMINS  -  47., 
HNS    47.  -  PMINS 
HEW  =  ( PLMINS  -  10.) 
HEW  =  (10.  -  PLMINS) 
*  2  +-  HEW  **  2) 
HDIST  =  SQRKHNS  **  2 


7  547 
7  54  7 


f  DDDIST  **  2) 


131 


c 
c 
c 


c 
c 
c 
c 
c 
c 
c 
c 


c 
c 
c 
c 
c 
c 
c 
c 


c 
c 
c 
c 
c 
c 


c 
c 
c 


c 
c 
c 
c 
c 
c 
c 
c 


COMPUTE  DISTANCE  TO  SCENE  FROM  CRESCENT  CITY 

IF  ((  PLAT. EQ. 41 .)  .AND. (PMINS. GE. 45.)  )  CNS  = 
1  (PMINS  -  45,) 

IF  ((  PLAT.EQ.41  .) .AND. (PMINS. LT.  45.) )  CNS  = 
1  (PMINS  -  45.) 

IF  (PLAT. EQ. 40.)  CNS  =  45.  +  (60.  -  PMINS) 

CEW  =  HEW 

CDIST  =  SORT  (CNS  *#  2  +  CEW  #*  2) 

IF  (DDDIST sGT.O. )  CDIST  =  SQRT(CNS  **  2  +  DDDIST  **  2) 
1  **  0.5 

CTIME  =  CDIST  /  18. 

COMPUTE  SAN  FRANCISCO  AIR  STATION  DISTANCE 

IF  SEVERITY  CODE  LESS  THAN  3,  AIRCRAFT 
FROM  SAN  FRANCISCO  WOULD  NOT  RESPOND. 
PROGRAM  BRANCHES  TO  ARCATA  AIR  STATION. 

40  IF  ( ( IPER.LT.3) .AND. ( 1PR0P.LT.3) )  GO  TO  60 
IF  (ISEA.EQ.8)  GO  TO  47 
SFEW  =  DDDIST 

IF  (PLAT.EQ.41.)  SFNS  =  PMINS 
IF  (PLAT. EQ. 40.)  SFNS  =  60.  -  PMINS 

IF  THE  DISTANCE  OFFSHORE  OF  THE  DISTRESS  IS 
LESS  THAN  100  MILES?  ASSUME  THAT  AN  Hri-52 
CAN  BE  USED  IF  IT  IS  ESCORTED  BY  A  FIXED-WING 
AIRCRAFT  FROM  SAN  FRANCISCO.  ONE  HOUR  IS 
ADDED  TO  TRANSIT  TIME  FROM  SAN  FRANCISCO 
FOR  THE  HELICOPTER  TO  REFUEL  IN  ARCATA. 

IF  ( {DDDIST. GT.O. ). AND. (DDDIST. LT. 100.)  )  SFTIME  = 


1  3.7  +  (SQRT( DDDIST 


2  +  SFNS  **  2)  /  85.  ) 


RESPONSE  TIME  FOR  SAN  FRANCISCO  HH-?3F  IS 
CALCULATED.   NO  ESCORT  FOR  THIS  HELICOPTER 
IS  REQUIRED. 


IF  (( DDDIST. GT.O. ) . AND. ( DDD I  ST . LT . 100 , )  ) PTIME 
1(SQKT(44100.  *  DDDIST  **  2)  /  120.) 


43 
45 


IF  ((  DDDIST. GT.O 
IF  (DDD  I  ST.  GT.O.  ) 
IF  (PLMINS.GT.5. ) 
IF  ( PLMINS.LE.5.) 
IF  (SFEW.GT.25.) 
SFDIST  = 
GO  TO  4  5 
SFDIST  = 
SFTIME  = 
PTIME  =  I 
GG  TO  49 


). AND. (DDDIST. LT. 100, 

GO  TO  47 

SFEW  =  ( PLMINS  -  5. ) 

SFEW  =  (5.  --  PLMINS) 

GO  TO  43 

SQRT  ( SFNS  **  2  *  SFEW  : 


) )  GO  TO  4  9 

*  0.7547 

*  0.7547 


2) 


SFNS  *  SFEW 

3.7  +  SFDIST  /  85. 

.75  +  SFDIST  /  120. 


CALCULATE  RESPONSE  TIMES  tQR    LAND  CASES. 

47  IF  (ISEA.EQ.8)  SFTIME  =  3.7  +  DDDIST  /  85. 

IF  (ISEA.EQ.8)  PTIME  =  1.75  -'-  (DDDIST  /  120.) 
IF  (ISEA.EQ.8)  GO  TO  49 

NOTE:  THIS  TIME  BASED  ON  3  HOURS  FOR  HELG  TO  GE1  TO 
SAN  FRANCISCO  FROM  SAN  DIEGO  OR  ASTORIA  PLUS  1 
HOUR  FOR  REFUEL  PLUS  2  HOURS  FROM  SAN  FRANCISCO  TO 
ARCATA.   MISSION  THEN  PROCEEDS  FROM  THERE. 
THIS  DELAY  NOT  NECESSARY  IF  HH-3F  HELO 
STATIONED  AT  SAN  FRANCISCO 


SFTIME  =  (SQRT (DDD I  ST  **  2  + 
PTIME  --  3.73  f  (SQRT  I  DDDIST  ■ 
49  CONTINUE 


SFNS  **  2)     /  120.)  +  6. 
*  2  +  SFNS  **  2)  /  12  0.  ) 
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c 
c 

C  COMPUTE  DISTANCE  FROM  ARCATA 

C 

C 

C  CHECK  FOR  LAND  CASE.   IF  SEVERITY  CODE 

C  LESS  THAN  2,  BYPASS  THIS  SEGMENT. 

C 

C 

60  IF  (ISEA.EQ.8)  GO  TO  65 

IF  ( ( IPER.LT.2) .AND. ( IPR0P.LT.2) )  GO  TO  00 

AEW  =  DDDIST 

IF  (PLAT. EQ. 41 . )  ANS  =  PMINS 

IF  ( PLAT. EQ.40.)  ANS  =  60.  -  PMINS 
C 

C  IF  THE  DISTANCE  OFFSHORE  IS  LESS  THAN  100 

C  MILES,  IT  IS  ASSUMED  THAT  A  HH-52  HELICOPTER 

C  FROM  ARCATA  COULD  BE  USED  IF  IT  IS  ESCORTED 

C  BY  A  FIXED-WING  AIRCRAFT  FROM  SAN  FRANCISCO. 

C  ONE  AND  ONE-HALF  HOURS  DELAY  TIME  IS  ADDED 

C  TO  ACCOUNT  FOR  THE  TIME  TO  GET  THE  FIXED- 

C  WING  AIRCRAFT  TO  ARCATA. 
C 

IF  ( (DDDI ST. GT.O. ) .AND. (DDDIST. LT.  100.) )  ATIME  =  1.5  + 
1  (SQRT( DDDIST  **  2  +  ANS  **  2)  /  85.) 

IF  ( (DDDI ST. GT.O.). AND. (DDDI ST. LT. 100.) )  GO  TO  6  9 

IF  (DDDI ST. GT.O. )  GO  TO  6  7 

IF  (PLMINS.GT.5. )  AEW  =  (PLMINS  -  5.)  *  0=7547 

IF  (PLMINS.LE.5.  )  AEW  =  (5.  -  PLMINS)  *  0.7547 
C 

C  IF  THE  DISTANCE  OFFSHORE  IS  GREATER  THAN  25 

C  MILES*  THE  COUNTER  IS  INCREASED  BY  ONE. 
C 

IF  (AEW.GT.25.)  NRFAR  =  NRFAR  *  1 

IF  (AEW.GT.25.)  GO  TO  63 

ADIST  =  SQRT(ANS  **  2  +    AEW  **  2) 

GO  TO  65 
63  ADIST  =  ANS  +  AEW 
65  ATIME  =  0. 

IF  (AEW.LT.2  5.)  ATIME  =  ADIST  /  85. 

IF  ((AEW.GT.25.) .AND. (AEW. LT. 100.))  ATIME  = 
1  1.5  +  ADIST  /  35. 

IF  (ISEA.EQ.8)  ATIME  =  DDDIST  /  85. 

GO  TO  69 
67  ATIME  =  6.  +  (SQRT  (DDDIST  **  2  +  ANS  **  2)  /  120.) 

IF  (DDDIST.GT.100.  )  NRHH3  =  NRHH3  +  1 
69  CONTINUE 
C 
C 

C  THE  VALUES  OF  TIME  ARE  MADE  VERY  LARGE  IF  THEY  ARE 

C  EQUAL  TO  ZERO  IN  ORDER  THAT  THE  COMMAND  TO  SELECT 

C  THE  MINIMUM  VALUE  FOR  RESPONSE  TIME  WILL  NOT 

C  SELECT  A  ZERO  VALUE 
C 

c 

80  HTYME  =  HTIME 
CTYME  =  CT1ME 
ATYME  =  ATIME 
SFTYME  =  SFTIME 
PTYME  =  PTIME 

IF  (HTIME.EQ.O.)  HTYME  =  999. 
IF  (CTIME.EQ.O.)  CTYME  =  999. 
IF  (AT  1ME.EQ.0. )  ATYME  =  999. 
IF  (SFTIME.EQ.O. )  SFTYME  =  999. 
IF  (PTIME. EQ. O.J  PTYME  =  999. 
IF  (HDIST.EQ.O.)  HOIST  =  999. 
IF  (GDI  ST. EQ. O.J  CDIST  =  999. 
IF  ( SFDIST.EQ.O. J  SFDIST  =  999. 
IF  (ADIST. EQ. 0. )  ADIST  -  999. 


c 

C  MINIMUM  DISTANCE  OFFSHORE  IS  CALCULATED. 

C 

DIST  =  AM  INI  <HEW,AEW) 

IF  (DIST. LE. 150. )  COAST  =  COAST  +  1 
C 

C         MINIMUM  TIME  TO  RESPONO  WITH  ARCATA  AIR 
C         STATION  IS  CALCULATED. 
C 

STIME  =  AMIN1(HTYME,CTYME,ATYME,SFTYME) 
C 

C  MINIMUM  TIME  TO  RESPOND  WITHOUT  ARCATA  AIR 

C         STATION  IS  CALCULATED. 
C 

TT1ME  =  AMIN1(HTYME,CTYME,SFTYME) 
C 

C  MINIMUM  TIME  TO  RESPOND  WITH  HH-3F  AT 

C         SAN  FRANCISCO  IS  CALCULATED. 
C 

UTIME  =  AM  INI  (HTYMEfCTYMEi PTYME ) 
C 

C         RELATIVE  EFFECT  OF  THREE  ALTERNATIVES  IS 
C         TALLIED.   OVER5  INDICATES  RESPONSE  TIME 
C  WAS  GREATER  THAN  ONE-HALF  HOUR.   OVER1 

C         INDICATES  THAT  RESPONSE  TIME  IS  GREATER 
C         THAN  ONE  HOUR  IN  THE  COASTAL  ZONE. 
C 

IF  ((DIST. LE. 150.). AND. (TTIME.GT. 0.5))  0VER5  = 
i  OVER  5  +  1 

IF  ((DIST. LE. 150.). AND. (TTIME. GT. 1. ) )  OVER!  = 
1  0VER1  +  1 

IF  ((DIST. LE. 150.) .AND. (STIME. GT. 1.) )  GREAT1  = 
1  GREAT1  +  1. 

IF  ((DIST. LE. 150.). AND. (STIME. GT. 0.5))  GREAT5  = 
1  GREAT5  +  1 

IF  ( (DIST.LE.150o) .AND.  (UTIME .GT. 0.5) )  MORE  5  - 
1  MORE  5  +  1 

IF  ((DIST. LE. 150.) .AND. (UTIME. GT.l.) )  M0RE1  = 
1  M0RE1  +  1 

IF(  ISEA.EQ.8)  LAND  =  1 

IF  ( STIME. GT. 100.)  STIME  =  0. 

IF  (TTIME.GT. 100.)  TTIME  =  0. 

IF  (UTIME.GT.100.)  UTIME  =  0. 
C 

C         SEVERITY  TO  INCIDENTS  IS  TALLIED. 
C 

IF  ( i IPER.EU.3i .0R.( I  PROP. EQ. 3)  )  GO  TO  8  3 

IF  ( ( IPER.EQ.2) .OR. ( IPR0P.EG.2) )  XMCD  =  XMOD  t  1. 

IF  (( IPER.EQ.2) .OR. ( IPROP.EQ. 2) )  T0T2  =  T0T2  +  TTIME 

IF  ( ( IPER.EQ.2) .0R.( IPR0P.EQ.2) )  T0TNU2  = 
1  T0TNU2  +  STIME 
83  IF  t ( IPFR.EQ.3) cOR. ( IPR0P.EQc3) )  XSEV  =  XSEV  +  1. 

IF  ((IPER.EQ.3)  .0R.(IPR0P.EQ.3))  T0T3  =  T0T3  <-  TTIME 

IF  ((IPER.EQ.3) .OR. (IPROP.EQ. 3) )  T0TNU3  =  T0TNU3  + 
1  STIME 

IF  (( IPER.EQ.3) .OR. (IPROP.EQ. 3) )  T0TH3  =  T0TM3  +  UTIME 


C 

c 


c 

c 


WRITE  (6,200)  J, HTIME, CTIMEyAT IMEjSFTIME, STIME, TTIME, 
1  LAND? IPERflPRUP 
200  FORMAT  ( 2X  ,  I  3  ,  i OX , F5 . 2 , 10X , F5 . 2 , 1  OX , F5. 2 , 10X , F5 . 2 , 1 OX, 
1  F5o278X,F5.2,8X. I  1 ,10X,211 ) 
WRITE  (6,220)  PTIME, UTIME 
2  20  FORMAT  ( 5  9X, F5. 2 , 9X, F6 . 2) 

WRITE  (7,222)  ST  I  ME , TT I  ME, UTIME , LAND 
222  FORMAT  (3F10. 2 , 1  OX , 1 1 ) 

TOT  =  TOT  +  TTIME 
TOTNEW  =  TOTNEW  +  STIME 

TCTHH3  =  TOTHHJ  f  UTIME 
990  CCN1 INUE 


13* 


c 

C  AVERAGE  response   times   are  calculated 


c 


c 
c 


AVOLO  =-•  TOT  /  NUMB 
AVNEW  =  TOTNEW  /  NUMB 
AVHH3  ---  TOTHH3  /  NUMB 

WRITE  (6,400) 
WRITE  (6,500)  AVOLD 
500  FORMAT  (//'  WITHOUT  ARCATA  CGAS  THE  AVERAGE  RESPONSE1, 
1'  TIME  IS  «,F5.2,'  HOURS.') 
WRITE  (6,600)  AVNEW 
600  FORMAT <//»  WITH  COAST  GUARD  AIR  STATION  AT  ARCATA', 
1'  AVERAGE  RESPONSE  TIME  IS  «,F5.2,'  HOURS.') 
WRITE  (6,550)  AVHH3 
550  FORMAT  (//'  WITH  HH3F  AT  SFRAN,  AVERAGE  RESPONSE  TIME 
Is  IS:  «  ,F5.2, '  HOURS. « ) 
WRITE  (6,700)  NRLAND 
700  FORMAT  {//•  THE  NUMBER  OF  CASES  ON  LAND  IS:  ',12) 

WRITE  (6,800)  NRFAR 
800  FORMAT  (//'  THE  NUMBER  OF  CASES  OFFSHORE  REQUIRING8, 
1'  ESCORT  IS:  • , 14) 
WRITE  (6,900)  NRHH3 
900  FCRMAT  (//•  THE  NUMBER  OF  OFFSHORE  CASES  GREATER  THAN' 
1'  100  MILES  IS:  • , 14) 

RES2  =  T0T2  /  XMOD 
RESNU2  =  T0TNU2  /  XMOD 
RES3  =  T0T3  /  XSEV 
RESNU3  =  T0TNU3  /  XSEV 
RES3F  =  T0TH3  /  XSEV 

RES23    =     (T0T2    +    T0T3)/(XM0D    +    XSEV) 
RENU23    =    (T0TNU2    +    T0TNU3 )  /  ( XMOD    *■    XSEV) 
RES23F    =    IT0T2    +    T0TH3)    /     (XMOD    +    XSEV) 
WRITE     (6*1400)     RES2,RESNU2 
1400    FORMAT     (//•     AVERAGE    RESPuNSE    TIME    TO    MODERATE-DANGER1, 
1«     CASES    WITHOUT    ARCATA    AIR:     «,F5.2,»       WITH    ARCATA     c, 
2'AIR    STATION:     «,F5.2,«     HOURS.') 
WRITE    (6,1100)    RES3,RESNU3,RES3F 
1100    FORMAT     (//•     AVERAGE    TIME    TO    RESPOND    TO    S EVE  RE-DANGER ! , 
1«     CASES    WITHOUT    ARCATA    AIR:     ' ,F5*2,6     WITH    ARCATA     8? 
2'AIR:     ',F5.2,'        WITH    SFRAN    H3F:     ',F5.2,»     HOURS.') 
WRITE    (6,1200)    RES23,RENU23,RES23F 
1200    FORMAT    (//'     AVERAGE    RESPONSE    TI ME-MGDER ATE/SEVERE » , 
1'     CASES     WITHOUT    ARCATA    AIR:     «,F562,!        WITH    ARCATA     S 
2'AIR    STATION:     ',F5„2,«     WITH    SFKAN    H3F:     «,F5.2, 
3'     HOURS.') 
WRITE    (6,1300)    XMOD, XSEV 
1300    FORMAT     (//'     NUMBER    OF    CASES         MODERATE:     «,F4.0, 
1«  SEVERE  :       '  ,F4U0) 

WRITE     (6, 1500)COAST 
1500    FORMAT     (//'THE    NUMBER    OF    CASES    IN    THE    COASTAL', 
1'     ZONE    IS:     «,F5.0) 
WRITE    (6,1600) 
1600    FORMAT     (//«     THE    NUMBER    OF    CASES    IN    THE    COASTAL    ZONE', 
1«     IN    WHICH    RESPONSE    TIME    WAS    GREATER    THAN', 
2'     ONE-HALF    HOUR    IS:     » ) 
WRITE    (6,1660)     0VER5,GREAT5,M0RE5 
1660    FORMAT     ('     WITHOUT     ARCATA    AIR:     ' yF4.0,, 

1'       WITH    ARCATA    AIR:     «,F4.0,'        WITH    SFRAN    H3F:     J,I4) 
WRITE     (6,1700) 
1700    FORMAT     (//'     THE     NUMBER    OF    COASTAL    CASES     IN    WHICH1, 
1»     RESPONSE    TIME    WAS    GREATER    THAN    ONE    HOUR    IS:     ') 
WRITE    (6,1700)     0VER1,GREAT1 ,M0RE1 
1760    FORMAT     ('     WITHOUT    ARCATA    AIR:        ',F4.0, 

1'        WITH    ARCATA    AIR:     ',F4.0,'        WITH    SFRAN    H3F:     «,I4) 
STOP 
END 
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